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Beitrag zur schiefen Kreiszylinderschale* 
Von E. Schatz, Miinchen 
Mit 7 Textabbildungen 


Zusammenfassung. Die mit Hilfe der Tensorrechnung in allgemeiner Form aufgestellten 
Beziehungen lassen sich fiir das vorgelegte Problem spezialisieren. Dabei fiihrt man die geometri- 
schen Gréfen in die Gleichgewichtsbedingungen ein und driickt die SchnittgréBen durch die 
Verschiebungen aus. Als Ma8 fiir die Schiefwinkeligkeit wird ein Steigungsparameter f angesetzt 
und die Verschiebungs- und Schnittgré8en in Taylor-Reihen nach Potenzen von f entwickelt. 
Fir f=0 geht die schiefe Schale in die bekannte orthogonale iiber. Eine Fourier-Doppelreihe 
liefert ein partikulares Integral, das die Belastung enthalt. Die Hinzunahme homogener Losungen 
gestattet die Erfillung der Randbedingungen. 


I. Allgemeines 


Die vorliegende Arbeit befaBt sich mit der Aufstellung der Gleichungen fiir die 
Verformungen und die SchnittgréBen einer schiefen Kreiszylinderschale. Diese Auf- 
gabe tritt an den Bauingenieur heran, wenn beispielsweise ein Verkehrsweg iiber einen 
zweiten mittels Gewolbe schrage iiberfiihrt werden soll. Vereinfachende Symmetrie- 
bedingungen sind dann nicht mehr vorhanden. 

Da die Schnittfithrung im Grundri& nicht unter 90° verlauft, andern sich nach 
Abwicklung der Schale zwei Berandungen nach einem einer Graden — der Schnitt- 
graden im gewoélbten Zustand — iiberlagerten Sinusgesetz. Das bedeutet, daB sich das 
Schalenelement laufend andert. Die iibliche Art, die Gleichgewichts- und Formande- 
rungsbedingungen am Element aufzustellen, erfahrt dadurch erhebliche Schwierig- 
keiten. Es wird daher von der Tensorrechnung Gebrauch gemacht, die es erlaubt, 
einen Zusammenhang in beliebige Raume zu transformieren. 

Das behandelte Problem la8t sich mit Hilfe der Ansitze von H. Neuber?? in 
iibersichtlicher Form darstellen. Dabei treten obere Zeiger auf, die nicht mit Potenz- 
exponenten verwechselt werden diirfen. Potenzen von Tensoren wird man in empfeh- 
lenswerter Weise in Klammern setzen. Zur Einarbeit seien an dieser Stelle die drei 
Bande von A. Duschek*und jener von M. Lagally* genannt. In englischer Sprache 
seien die Werke von Green-Zerna und Shouten erwahnt. 

Fiir die Bedingungsgleichung gewinnt man ein partikulares Integral mit Hilfe 
von Fourier-Doppelreihen. Sie erfassen die Lastglieder. Zusammen mit dem allge- 


* Auszug aus der von der Fakultat fiir Bauwesen an der Technischen Hochschule Minchen 
genehmigten Dissertation. 1. Referent Prof. Dr.-Ing. H. Neuber, 2. Referent Prof. Dr.-Ing. 
habil. G. Worch. 

1 Neuber, H.: Allgemeine Schalentheorie. ZaMM 29, 97—108 und 142—146 (1949). 

2 Neuber, H.: Aufstellung und Integration der Grundgleichungen der biegesteifen Kreis- 
zylinderschale mit Hilfe der allgemeinen Schalentheorie. Wissenschaftliche Zeitschrift der TH 
Dresden 2, 511 (1952/53). 

3 Duschek, A., und A. Hochrainer: Grundziige der Tensorrechnung in analytischer Dar- 
stellung. 3 Bande, Wien: Springer. 1954, 1950, 1955. 

4 Lagally, M. und W. Franz: Vorlesungen tber Vektorrechnung, 5. Aufl. Leipzig: Aka- 
demische Verlagsgesellschaft Geest u. Portig. 1956. 
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meinen Integral, das die ,,offenen‘‘ Randbedingungen befriedigt, gelingt die voll- 
standige Lésung. Man bendtigt zwei komplizierte Funktionen F und Fz. Die erste 
gestattet die Erfiillung der Verschiebungsbedingungen langs des waagrechten Randes 
(Auflager), die zweite bringt die SchnittgroBen am gewdlbten (freien) Rand zum Ver- 
schwinden. 

Oben wurde erwihnt, daB das Schalenelement sich im Fortschreiten von einem Auf- 
lager zum anderen andert. Das findet mathematisch seinen Niederschlag in den ver- 
anderlichen Koeffizienten bei den Bestimmungsgleichungen, die eine strenge Losung 
sehr erschweren. Mit Hilfe der Stérungsrechnung®® kann man die Steigung f als Sto- 
rungsparameter auffassen und die gesuchten Groen in Taylorreihen nach Potenzen 
von f aufbauen. 

Die Absicht, gebrauchsfertige Formeln vorzulegen, scheiterte daran, dab die 
Funktionen Fr nur iiber die kubische Resolvente einer Gleichung vierten Grades zu 
ermitteln sind. Die Lésungen haben Zwischenkriterien zu erfiillen. Numerisch ist der 
aufgezeigte Rechnungsgang grundsatzlich gangbar, wenngleich er erhebliche Arbeit 
macht. Mit einer elektronischen Rechenmaschine diirfte die Lésung fiir alle verainder- 
baren GrdBen ohne besondere Miihe aufzufinden sein. 


II. Geometrie 


Wir legen das System nach Abb. 1 und 2 vor. 


x=m+E, 
y=rcosy, {| (1) 
2=7 sin'@, | 


Nach m, gy und r aufgelést: 


m= a —2t ; 


—arctg — 
Be teat (2) 
r = const., 

l=2rsiny. 


z und / in (1) eingesetzt, liefert 


sin 
C= mM eas 


2 sin y ’ 
y=? COs @, (3) 
of sI0 Y, | 
Z=2g= a8, (3a) 


° Collatz, L.: Numerische Behandlung von Differentialgleichungen, 8S. 80. Berlin: Springer. 
1951. 


* Courant, R. und D. Hilbert: Methoden der mathematischen Physik, 2. Aufl. S. 296. 
Berlin: Springer. 1931. 
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dann folgt allgemein 


a, = x = wt (a*) = % (2*); k= 1, 2,3 (3b) 
Co lhe 


«* = a, erfafit das kartesische (orthogonale), «* das beliebige System. Da in der Aus- 
rechnung bei beiden Systemen gleiche Zablenzeiger auftreten, unterscheiden wir, wo 
es notwendig wird, die kartesischen durch Uberstreichen. In der allgemeinen Form 
wollen wir lateinische Zeiger (k,1,m...) im kartesischen System und griechische 
(x, B, y, 0...) im beliebigen verwenden. 
Zur Transformation bendtigt man die Ableitungen 

0 EE " ce eR 

eae oe a= 1,2 
In Tabellenform 


ae aE eg adi (5) 
am ow ke op ee 
9 On  Ff.008P  9Y az 


=—rsin 9 ap 1 OOS: 


dp 2siny dg F 


Wir wollen noch den Hinweis machen, daf die kartesischen Zeiger k gemaB (3b) 
oben und unten stehen kénnen, ohne den Tensor zu verandern. Das gilt fiir die Zeiger 
a, B,... nicht mehr. Sie miissen zur Darstellung von kontravarianten Tensoren oben, 
bei kovarianten unten stehen. Wir wollen die ersten bevorzugen. 

Der metrische kovariante (auch Fundamental-) Tensor der Schalenflache ist defi- 
niert durch 

CH CH= Gap: (6) 


Uber doppelt auftretende Zeiger (hier iiber /) wird summiert. 
Ausgerechnet in Matrizenform 


l j cosg 
2siny 
tap = (7) 
feosp , , ficos®y 
2siny se 4sin? y 
Die Determinante der a@,, ist 
Det da,=7". (8) 
Aus 
= lira y 
qe8 t..= 6% { 9 
By Y | =0 fiir a+ y (9) 
erhalt man iiber den schiefsymmetrischen ¢e-Tensor 
1 
\a — fira=1, p=2 
Va 
Eup = O- 7. eth 0 fir «— (10) 
—|a ow Ae tn ga on 1 


den kontravarianten metrischen Tensor 
Gok eee a. (11) 


15* 
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Ausgerechnet 
fr cos? gy feosp ° 
ei 4+ —— — 
2D v2 
pe 4r2 sin? y 7? sin y (12) 
—fcos@p at 
272 sin y r 
Aus Gl. (2) findet man die reziproken c¢ und kann damit 
c* cha? 
bestatigen. 
Wir errichten auf der Schalenmittelflache die Normale 
Ny= NF Sh th im EtPof OF: (13) 
Man errechnet 
N,= (0, —cos y, —sin ¢). (14) 


N*‘ = N; ist demnach ein Einheitsvektor. 
Der kovariante Kriimmungstensor der Flache (Schalenmittelflache) ist gegeben 


durch 
k 


<= 15 
Da B ax - he) 
In Tabellenform folgt 
siamese ] 2 3 
a 
ack 
AE 0 0) ) | 
; (16) 
act a cf ‘ 6 ; { 
0 x be 0 x 
a 3 fsin y Res 
ye aay rcosg —rsing. 


Im pao nonnaclen Raum (Schalenmittelflache) andert sich die Basis der Tangentenvek- 
toren cf, von Ort zu Ort. Da er im kartesischen Raum eingelagert ist, bezieht man 
die Anderung auf diesen. Die dazu notwendigen Christoffel-Symbole zweiter Art sind 
wie folgt definiert: 


ck 
a 9&8 


= | (17) 


aa 


Die gemischten Differentiationen verschwinden, so da alle Christoffel-Symbole auf 
der Schale Null werden. 
Der kovariante Kriimmungstensor lautet ausfiihrlich: 


Uton. 
baa=(j a . (18): 


Den gemischten Kriimmungstensor erhalt man aus 


bg =a", , (19) 
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se ZO) 
ei 1)! OF 70 (20) 
: f cos 1 

2) Nema 


a 
B 


7at= At 2 
Es ist: 63 + bj, wogegen b¥2= 621, b,, =o, dy.=d, USW. 


Aus 
ip sha adele es (21) 
folgt 
f? cos? feos 
hob — 473 sin? y 273 sin y (22) 
—f cosy 1 
273 sin y 7 


Die mittlere Kriimmung folgt aus 


ee ae Se a ; 
Ha 1 x)= 2 De ae (29) 
und die Gau&sche aus 
] 1] a 
K=-,-— =Detbe=0. (24) 
Cr Qe 


Hier soll nur der EinfluB des Eigengewichts untersucht werden. Der Kraftvektor 
hat dann die Richtung der negativen y-Achse: 


p= (0, —7, 0), v=], 2,3. (25) 


Die Zerlegung in einem beliebigen Punkt 
mit den Grundvektoren 


Mm; = ( i 0, 0), 
Qi= (Snr , —rsing, 7 cos 9), (26) 
fo) — cos g,—sin ¢) 
fordert Abb. 3 
p=AmM+ UG trYT: (27) 
und liefert* die skalaren Werte 
fcos psing 
he De rainy ” 
ak a (28) 
y= pcosg. 


* Vgl. FuBnote 3 auf 8. 193. 
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Die Betrage der Vektoren (26) sind 
Imi|=1, 


loi] =/r2+ Proos? p (29) 


4sin?y 


|7;| =) cos? y + sin? p= 1. 


Die Grundvektoren (26) enthalten neben den Tangentenvektoren m; und g; den Nor- 
malenvektor 7;. Nach (4) haben wir fiir die Zeiger «, f,... nur die Zahlen 1 und 2 
zugelassen. Wir sind daher gezwungen, fiir die Normale die Zahl 3 hinzuzunehmen 
und dann die allgemeinen Zeiger /, uw, » statt «, B, y zu beniitzen, die die Zahlen 1, 2 
und 3 annehmen kénnen. Auf diese Weise erweitern wir den Flachenvektor der Schale 
zum Raumvektor. Dann erscheint die Belastung in folgender Darstellung: 


COS y Sin ~ 2 cos? p 
| p eee : psin g J/1+ faa ‘ p cos 9} (30) 


Der Gleichgewichtsbetrachtung wird ein kontravarianter Spannungstensor zu- 
grunde gelegt. Wir wollen ihn voriibergehend w*® nennen — sein kovarianter Gegen- 
spieler sei w,3; — und seine Beziehungen zu o,; im kartesischen System angeben. 


zh 
ee 


Abb. 5 lehrt aus der Anschauung daB ui? =u?! sein muB, da uw! und u?? durch den 
Schwerpunkt des schiefwinkeligen Elementes gehen und fiir das Torsionsmoment 
um den Mittelpunkt § keinen Betrag liefern. Das la8t sich bei schrittweiser Aus- 
rechnung bestiatigen. Man erhalt: 


u Ox; Ox 
Opp = UB CMe eek 
ea ane aab (31) 
Oxy 82 
uy =a ik OR; OXk (32) 
u LL Gan are 
af q  “y “Bd (33) 
1 Or, OX, 
U = — Peer es a 
gen Pa aah (34) 
=a ax% aah 
Yo Oa; Oa; (35) 
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6 ae 
Die u, und w® sind reziproke Systeme; ebenso die dazugehérigen Tensoren wu, baw. w°. 


Vel. aN 6 und 7. 


Abb. 7 
Uy, Lut 
(36) 
Us, Lut; Us. Lut 
Fiir die Randnormalen folgt langs m= const. 
ee ry ae feos ysin yp f cos? p 
mi = 830 9; Ne=(r, 2siny ° aes | 
i 2 eos? y (37) 
rn, | =]/ 9 at emis 
und langs gy = const. 
nn; = &;;, m; N;, = (0, sing, —cos 9), 
jk ™; ( x. ?) | (38) 


| ®n;| = //sin? y + cos? gy =1. j 

Die Spannungen w*®, iiber den Querschnitt von —h/2 bis + h/2 integriert, ergeben 

die Langskrafte L**. Sie stehen nicht senkrecht auf ihren Flachen und werden daher 

auch nicht als Normalkrafte bezeichnet. Vgl. Abb. 5. Die Langskrafte L*® haben mit 
beliebigen Skalaren a, b, c und d die Form 


DP=a( 1, 0, am 
fcosy. ae 
rao(ee, rsing, frcos@), 
(39) 
LXx=c( 1, 0, OF 1). ae 
\ feos aves 
rm aa(bee, rsing, 17 COSq@j. 
Zerlegt man sie in Normal- und SON SEY so folgt 
NG 40 
N=. ] c ee er pr (40) 
\/" a5 4sin?y 
mN, =a a +b, (41) 
2 sin / +f id 
? 4sin?y 
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°N,=—ad mil) (42) 
/ 1+ f? cos? p 
47° sin? y 
feos y 
ON, =c+d 2rsin y (43) 
[ feostp 
y2 + 47 sin? y sin? y 
Die Transformation des Laplace-Operators liefert * 
ae@ , 3@ 
A®=a ae anB foal (44) 
Ausgerechnet und vereinfacht 
ap oP 
Ag= 0 xo 7 (a 35): (45) 
Mit ¥ d— cs ai): oe folgen die Varianten 
: OG ( am 4 
ag @@ 
a ea ae =AGb—YV@, (46) 
oz 
szaaue (A D)=AD+VS. (47) 


III. Gleichgewichts- und Formainderungsbedingungen 


H. Neuber** geht von einem kontravarianten Spannungstensor aus, der auf der 
Schalenmittelflache wirkt. Durch Entwicklung in eine Taylor-Reihe in Richtung der 
Normalen gelingt es, die Umgebung der Schalenmittelflache tiber die Dicke hf zu er- 
fassen. Die Integration von — h/2 bis + h/2 liefert die Langskrafte L** und die Momente 
M*®, Nach den Gl. (23) und (25) bei Neuber*** hat man 


a0), — bg MY, — 68 Ml, 4 6; MeO), + p= 0, (48) 
M*|,,+6,, °° KM? + p?=0. (49) 


L*°|, bedeutet die kovariante Ableitung nach x*. Wegen (16) und (17) geht die kovariante 
Ableitung in die partielle iiber. Man erhalt: 


aI) ol 10M" | -feosp oM™ 


& yd 
am op r om 2Pin yea g. we p> (50) 
OL el pane Me Py ae Peon Ta Mass ge 51 
am ay r om r op 2rsiny am r? (51) 
ez Miu oz M22 oz M22 
eee 8 
omit smedp i tog Te oe (52) 


Die technischen Schnittgr6Ben (iiberstrichen) hangen mit den tefsoriellen wie 
folgt zusammen 


N,=14, N,=rL", N=? L, (53) 
M,,=M4, M,,=rM, M,,=r? M%. (54) 


* §. FuBnote 3 auf S. 193. 
** S. FuBnote 1 auf S. 193. 
*** S. FuBnote 2 auf S. 193. 
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_ Die Verschiebungen sollen sein * 


in Richtung m; ‘i oe ee awe 
qu / f? cos? p zi 
yp ¥ 47° sin? y 
in Richtung ¢; ee oat —rV2 ee 
Vaz 2 
in Richtung —7; Ww = V3; 
Dann wird 
_ FI Poste 
i bl | 1 47 sin?y U, 
ye V (56) 
V2=W. ) 
Fiir den Zusammenhang von L*? und V* erhalt man** 
apis a a} > ha : a 
L*6 —Gh last Ve, + at® Vol, — 26° V8 4 —2— ash (Vy, —BxV8)], (57) 
ausfithrlich 
PSG 2a os +2 oe a ae et P|, (58a) 
I= Gh la 2 4 an 2F + anso + amto _opew4 -4 ao], (58b) 


aV aV Qu 
22 1 A SrAO Neti oa Gy) WY) 
L Gh |2a re es ee ae 2b eg) 


am | (58c) 


bu bedeutet die Querdehnzahl. Weiter wurde zur Abkiirzung 


oU aoV 1 
ame Yr Q iGee (59) 


gesetzt. 
Fiir die Momente folgt 


he 
M8 = —G= \at? aP® V3.4 at? V4, 4) . (60) 
Ausfiihrlich mit V?= W 
h a ug ew E > 02? W be 
Mu=—@5|@ 1p 2 a 2 + ose toy 4W—PW)), (618) 
h3 ew ew EE oo? W 
Mi2— Ge & qt ail + a}? On a5 pe qil q22 am og + ql2 G22 ap a . 


tats H (4 ean 


ew ew 
Me=—@F|(a aye OW 420 atl ai oo + (a2)? a2 (AWW). (61) 


Die Querkrafte ermittelt man aus 
Q=M"|,, (62) 


* S. FuBnote 2 auf S. 193. 
** Aus 1, 2 auf 8. 193. 
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so dak 
= CRU PES TUES 
Oi = =e arte (63a) 
— aM 0 M22 
Q,='C-=r a9 karen (63b) 
wird. 


IV. Aufstellung der Verschiebungsgleichungen 


Die Gl. (58) und (61) werden in (50) bis (52) eingesetzt. Nach Ausfiihrung der 
Differentiationen und Einfiihrung der Werte in (50) bis (52) folgt bei Vernachlassigung 
der Momentenanteile in den beiden ersten Beziehungen, wenn der Laplace-Operator 
nach (45) verwendet wird, 


Agee 2 (un 22) 4 a (a2 22) 4 2 Iwo (222) | 


l1—pom am} ' am ap l—yp op am 
la ( na) lgeuae ou ee eee us 
" rap am) ' l—p dy ~ oe) — sals a4 (64a) 
a 
: (20), Les 7 0 {au 0 (a2 | 
Pg Cale dean ke we | “|s(4 w)+3.(Cw)|| 
(aes 
Tayany =0, 
F) aU F) aU 2S aU r) aU 
Pr ays aca Ce esr dacs a ueey ad Toles 
t) av 2u {9 av \\ l-+p @ aV 
SAM n gp (ra?) . feplawl?t ae ae lt ae) (64b) 
tg F) F) F) 
— 2} (r 612 W) + aa (rb W) + ar E (BAW )i-E ae (a2 w)|t a: 
fae 
2 on 22 
oV 2 aoV 
age ule (640) 
he : 
+ |AAW+ PEW Ot, (A+) (AW— YW)| + 
‘ 2 fl p* 
ae Oo Ty at Tone ay Oe 
Das System (64) 14Bt sich kurz schreiben 
a Ras 
Dy VE + = , Samy 
Vit p;=0 ‘ten (65) 
Ks enthalt die Abkirzungen 
) a ah a2 
ww (.)= (Far). 2 (...), ; (66) 
F) 
O (2. )=aP se (a) 55 (++) $a Ze OM) 55 (...) + 
F) a2 


Guages (67) 
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V. Ermittlung eines partikuliren Integrals mit Hilfe einer Fourierdoppelreihe 
fiir die vorgegebene Belastung 


Das partikulare Integral enthalt die Belastung p’, befriedigt aber nicht die Rand- 
bedingungen. Es wird die Addition eines entsprechenden allgemeinen Integrals not- 


wendig. Zunachst wollen wir die Gl. (64) verwenden und die folgende Doppelreihe 
vorgeben: . 


Peep ty ; 
O meh = £00 + x >) a;,cosajm+ 4 >” c;, sin am + 
j 7 


& 
at on COS Bn p + _¥’ S’ajn Cos am cos Brg + 
n j n 


+ >») S' ¢insin am cos Bap + 4S” bon sin Bag + 
J n n 

+ >) Dd) bin cos ajm sin Bap + > DG, sin am Br G. 
J n Jj n 


In Richtung » wahlt man die Grenzen von —y bis + y, in Richtung m von — 2d bis 
+ 2d. (Mit —dbis +d werden alle trigonometrischen Glieder Null, so daf folgt p = p.) 
— Ks 14Bt sich jede Belastung durch (68) ausdriicken. — Durch die Grenzen ist iiber «; 
und f, verfiigt wegen 

sin oj 2d=sinjx wird a =j5q 
und snf,y =sinnzx wird £,=n : : 
Fiir f=0 treten keine Koppelglieder auf, so daB wir im Orthogonalen iiber m und » 
getrennt integrieren konnen. In Richtung m allein folgt nach? wegen j= 2d 


1 
cag) 


ee eee COs aM, g=1, 3,5... (69) 


Die sin-Glieder verschwinden wegen der Symmetrie. 
Fiir die Integration tiber gy formt man die Lastglieder aus (30) um 


: feospsing _ f sin 2 
Pome. Drei 4rsiny’ 
: 7? cos? p 
2— sin yi SE a ae 
P p ? oF 472 sinty (70) 


; 2 Cos? f4 cost p {8 cos’ p 
=psin (1+ sorcm 12874sinty © 16-647r6sin®y ease 


p>? =p cosy. 
Ersetzt man in p? die Potenzen des cos durch Ausdriicke von mehrfachem Argu- 


ment, so erhalt man 
p?=psing (ky +k, cos2y9+k, cos4y+...) 


mit 
f? 3 f4 
ky=14 16r2sinzy 1024r4sinty ’ 
73 ee el ie 
k= 16 72 sin? y l 16 72 sin? y} ’ (71) 
Wee 
1024 74 sin4 y 


7 Worch, G.: Betonkalender 1953, II. Teil, S. 92, Tafel 13. Berlin: Verlag Ernst & Sohn. 1953. 
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Nun einige Integrale fiir die weitere Ausrechnung 


‘ i sin (a+b) a sin (a — b) x 
J sin ae sin be de = raw ys) 
heer cos (a + b) x cos (a — b) & 
J sin aw cos ba da = — “% ae 8 ee ies 
; ‘ cos (a + b) x cos (a — b) & a 
| cos aa sin ba da — ile) Saas 
—b 
{cosax cos ba dx = 4 ey ra + i ats oe 
My | ese 
| cos? ax dx = 5 + 4 sin 2ax | 
ve} 
{sin axax = eagle 2ax 
[x cos dx =cos x4 «sin « 
(74) 
Jesin xdx=sin «—2 cos x | 
¥ ; 
ay sin? Bgl, 
eae ee ihe wae 
¥, 
f > Bn P _ wine? fe I (75) 
feos er jooneg ef arene pn es 


f-\l a oe ee | 
+2 0 —2Br o+2Bn 


+ 
1 [ sin y cos op sin B, p dy “= 
-y 
__ 9 (Pr? — Bn — Bng) sin y cosey cos Buy BNE 
Pi 2pe(l+e)+U—e8 ~ pe— PPA +o) +U—epe (76a) 


any, 


| I fsing cos ep cos B, yp dy =0. 
=i 


(76b) 


Fir V und W setzt man Reihen gleichen Aufbaues, nur mit anderen Koeffizienten an, 
symbolisch 


) (77) 
Ww > ee 


Man bildet von U, V, W die Ableitungen und setzt sie in (64) ein. Man erhalt dann 
mit (12) 


4 


——————— 


oe 
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2 a (1+ me) 3— mu f cos au i aU 1+ yn fsng 0U 
47? sin? y |] a m2 1— yp 2rsiny amragp 7 Ag nN ag? 2rsiny Sepia 


+724 (1 si ram av ltyu feosp eV 2fcospsing a@V 
4rsinty/amrag 1l—yp 2rsiny Pag? 4r2sin2y ram a 


by hoes fsng oaV | 
l—yp2rsny reg | Oo) 
eehe2 P? cos? p OW... 26 OW, 4 2 feosp oW 2 fsng W 
l—yp 4rsinty ram l—pn ram l1—p2rsiny P89 1l—4p 2rsiny 72 = 
_ pfsin2y 
~ 4Ghrsiny’ J 


l+yp feosg @2U- ltu #87 
l—yw2rsiny dm? 1— uw Omrdg 


+(1+ een pa 3—" foosg @8V 1 2 aV Easy, 1leV 
4r-simzyjam? 1—yp 2rsiny amrog ' 1—p rag? | Irsiny Pom 
2 feosp OW 2 aw 


: | : ao = 


l—yp2rsny ram l—uw rap 


= — yz Bin. (hy + he cos2¢9+k, cos4g+...), 


—2u 9U 
1 ram 

—feosg aV 2 dV 

2rsiny rom Il1—yp rd ! 

78¢ 

Coe abe LOL e, 9 aw ogee pares 2s ee) 
' 6 1l—wlamt ' “amirag? ' rag mae l—p r 
__ peosp 
PAG h 


re 
Fe 

~J 

@ 

og 

— 


Hier sind die Ableitungen von (68) fiir U und nach (77) die anderen Ableitungen 
V und W einzufiigen. Die Integration tiber gy von —y bis +y und iiber m von — 2d 


bis + 2d liefert nach et Zwischenrechnungen : 


Tag Cin + Ge BE Gin t TE Ft Ba fin — Te | 85 Din = 0, 
kurz 
041 Cin t+ O12 fin + 13 Pin = M1a= 9, (79a) 
analog fiir die Gleichungen (78b) und (78c) 
Oey Cin + Mae fin + %23 Din = M20; (79b) 
O31 Cin + O38 fin + X33 Din = Ta; (79c) 


‘wo 


hee tine od bas 
=p + ae Pn 


1+u 1 
— —, & Pns 
fale # HP (80) 
Dit rmeed 
13 — mm Ph J 9 


Ta 0 ? 
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‘lp 
Stee oH Hi Bns 


2 
Xen = O + aly 


(81) 
2 
Cr eae payne. 
2p 1 ‘ 1 {__sin(1+fn)y , sin(1— Br) | 
St Le pera ty; aaa 
Ha = Gh a i yd | 1 + Bn r ee ee 
od 
Su Ea ge eke 
2 
92 ae Bn 
h2 ] 2 \2 2 1 Sat 
, 2 \9 
%33 = § ale | ) iia pe? 
2p (—1)24-) | sin(1+Bhn) y | sin (1 — Pn) 7 


F fon 


a Oh aj yd ae be ey 
Die Gleichung (78c) liefert, mit cos £, py multipliziert und integriert, noch ein von Null 
verschiedenes Integral mit zwei neuen Konstanten. Eine zweite Gleichung dazu liefert 
die erste Gl. (78a), so da man ein weiteres System von zwei Gleichungen erhalt. 
Die Anteile der Gl. (78b) verschwinden. 


P11 Go tP13 Dip = M15 = 0 


(83) 
Yai Cig + O33 Dio = %sb5 
wo 
bu 
Oy, = 4, As = me 
a9 rps (That gS: 4 1 1 
tan Vag, re Os witout l—-anr’ (84) 
oY _p 2siny % G1) 

Tih a= Or 7136 = Gh ajpyd ( 1) ? 


Aus (79) und (83) lassen sich die Vorzahlen ausrechnen und in (68) einsetzen. Es folgt 


liyo , = 
Or gia Dy Co Sie so 2: > Cin SiN a M COS Br—, 
Jj 


Ves = ss D> fin COS a m sin B, —, (85) 


J n 


ee Re 
W (n,Q) = ae Dig COS a ™m ae D) Pin CO8 a; M COS BrP. 
ij eae 


Dieses Ergebnis hat die Form der Gl. (92) aus’, wenn man beriicksichtigt, daB dort 
vom Rand, hier von der Mitte aus in Richtung m (dort x) gezaihlt wird. Das bedeutet 
ein Vertauschen von sin «%;m mit cos am und umgekehrt. AuBerdem ist dort auf 
die trigonometrischen Anteile mit 8, verzichtet und Aa/a fiir om gesetzt. 


8 Worch, G.: Betonkalender 1958, II. Teil, 8. 38. Berlin: Verlag Ernst & Sohn. 1958. 
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Verschiebungsgleichungen 


nach Potenzen von f vor. 


ist zunachst f=0 zu setzen und anzusetzen 


ee oh Ua 20, PU, i Se PU, 
Vea ot Vy et Ve i ee = DY, 
W=W,+fW,+/? WET? Wastes. t DFT Wr. 


Kine solche Entwicklung wird als méglich vorausgesetzt. 
Mit f=0 folgt aus (64) 


= 244-00 Doser -Vaig Lh? 


_ VI. Aufstellung weiterer partikulirer Integrale fiir die verkiirzten 


AU 4 tare? U ave Taw 02 Vi epLen O Mees 0 
l1—pom '1l—pamreg l—pr am ~ 
l+p #U ! ‘pigeon av 2ET AW i paing 
l—y dmrag l—yp rodg l—prrep Gh 


, 
| 
l—p rom eerer eee 6 |l—wu SN le 


auftreten, und liefert 


h2 r2 a4 F 
paw 4444Fh + am = 0. 
Die Lésung lautet 
F= ¥V Yc, ev? cos am ig isha: 
erate : dg Be 
mit 
/ <1 ; oY hz 
alae eee yg LENZ? k = 12 (1 — p2) 12? 
ausfiihrlich 
F= N’cosajm [ o% {A, cos »,p + Asin  p} + 
4 +e4?{A, cos », p+ Aa sin »p}+ 
‘ 
+ elle? {As cos v2 9 + Ag SiN ¥9 Y} + 
fete? {Az cos 7, 9 + Ag sin vePt |; 
wo 


es /2(V/e +75) + (es) t#+ a) 


9 Girkmann, K.: Flachentragwerke, 3. Aufl. Wien: Springer. 1954. 
* 8. FuBnoten 5, 6 auf 8. 194. 
** S FuBnote 2 auf S. 193, FuBnote 8 auf S. 206 und Fu8Bnote 9. 
*** S FuGBnote 2 auf S. 193. 
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Hay ee 
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Das homogene System (64) wird mit Hilfe der Storungsrechnung* gelést®. Wir 
fassen die Steigung f (Abb. 2) als St6rungsparameter auf und gehen schrittweise 


Fir f=0 ist die Losung bekannt**. In der vorausgegangenen Fourierentwicklung 


(86) 


(87) 


Die Rechnung lauft wie bei Neuber*** weiter, da die Gl. (87) in gleicher Art wie dort . 


(88) 


(89) 


(90) 


(91) 


|e ia) Gal) 
V0 +a Ga) "a 
- (Van) * a) a) 


Von Neuber* iibernehmen wir 


be 
Vy 
Vo 


— 2 03 

Vy = Uh = — omg? ’ 

it. 29 ¢@ oF 3 
Vy =Von =a oe (4F + FS), (92) 


Der zweite Zeiger h soll auf die homogene Losung der Gl. (78) hinweisen. Zusammen 
mit (85), die wegen der enthaltenen Lastglieder den Zeiger / erhalten sollen, folgt 
als Losung 


U., a Osi ae Os, | 
ss — Vo: ok Vaks (93) 
Welw Se was | 


Jetzt sind die Randbedingungen zu befriedigen. Mit Uo; aus (85) und U >, aus (92) 
wird aus (93) fiir das Glied jn: 


: 2 : 3? 
U, =¢j, Sin 4 m Cos Bap + > y % sin Oj 5 oa [Pp], 


: 2 7) : es 
Vo =fjncos a m sin B, p+ ={- 2 a? cos a m oo [py] + cos a m-— Ip [v}} 
(94) 
e2 
W = Din Cos 4m cos Brgy + 2 {o4 cos am [p] — 2 «? cos am Roe [o]+ 
of 

== CORIO; ™ aga [e]}. 

[p] bedeutet die eckige Klammer aus (91). 


Als vierte Gleichung kommt bei gelenkiger Lagerung an den Kampfern M22 
aus (61c) hinzu (f= 0). 


Gh f uw #®@Ww 1 #ew 
22 
M™ — 6 r2 (Tp pas aa 


Fir y= + y fordern die Randbedingungen 
Oo Vy We Ma: (95) 


Damit hat man acht Gleichungen fiir die GréBen A,... As. Sie erscheinen fiir jedes n 
als Ausdriicke von Cjn, fin Und pjn, da sich die sin «m und die cos am herausheben, 
die sin B, p verschwinden und die cos £,m --1 oder —1 werden. 


* S. FuBnote 2 auf S. 193. 
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Nun sind noch die Bedingungen lings der Rander m= + d zu erfiillen. Wir gehen von 
einem Ansatz der Form 


Fr= PIPE i eev? ebrm (96) 
aus und fithren ihn in (88) ein. Dann entsteht 
here fo : 
mae +e} +e = : (97) 
Die @, sind nach (90) bekannt, wahrend die ¢, zu ermitteln sind. Mit 


C¢,.r=e und 02+ ¢«2=z 
wird daraus 
h2 1 ; . 2 i. 
7? 12(1—p2) * +z == 22 0, +0, = 0. (98) 
z laBt sich nicht allgemein anschreiben, da es iiber die kubische Resolvente von Zwi- 
schenkriterien abhangt (vgl.1°). Es 148t sich aber in jedem Fall numerisch ermitteln. 
Der Form nach gleicht die Lésung dem Ausdruck (90), so daB wir ansetzen kénnen 
Q...4=+4,+7 7, 
1 i1 (99) 
Cs...8= + Ye+7 Ho- 
Man erhalt als Lésung fiir #z, wenn man zunichst aus (91) nur das Glied 4, heraus- 
greift, das in seinen Faktoren die Gré8en », und », enthalt, 


Fr=e%® cos 9 [ e&” {Iy, cos ny m+ I,, sin n,m} + 


{ 
--e% {T5, cos n,m + Ig, sin yy m} + 
+ et2™ {151 cos n,m -+ Ie: sin 7, m} + 
+e {Tz cos n,m -+ Ig: sin 4 m} |. 


(100) 


Fir A, wire beispielsweise der zweite Zeiger bei den J durch 3 zu ersetzen. 
Nach Ausfiihrung der Operationen (92) addiert man (100) zu (94) und erhalt 
Uy =Vo1 + Von + Vor, 
Vo = Vor + Von + Vor, (101) 
Wo= Wot Won+ Wor. 
Mit Hilfe von (58) erhalt man die Langskrafte, mit (61) die Momente. Zur Bestim- 
mung der acht J,: stehen die acht Bedingungen 
Me Le = MVS MY =0 fiir m=+d (102) 


zur Verfiigung. Der Rechnungsgang soll fiir Z"' im Prinzip erlautert werden. Bei 
f=0 wird nach (58a) fir L4=0 


aU av 
am | ¥ rag iF 


Schreibt man [m;] fiir die eckige Klammer von (100) — fiir die verschiedenen J entspre- 


10 Schleicher, F.: Taschenbuch fiir Bauingenieure, 2. Aufl., 1. Bd., S. 47. Berlin: Springer. 
1955. 
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chende Zeiger — so erhalt man, wenn man sich auf A, und U beschrankt, aus (94) und 
(100) iiber (92) 


: a 
0 = Cin oj COS a m COS Bry + —, %} COS a; ap {et? A, cos »yp} — 


2, &02 


7 dp 


{er A, cos 7,9} ne [m,] firm=+d. 


Hier sind in [m,] nur die J,, (e=1, 2,...8) unbekannt. Mit den Bedingungen (102) 
lassen sich aus A, die J,,, I,, bis Jg1, aus A, die 1,2, I. bis Jg2 usw. errechnen. 
Setzt man fiir U, aus (93) Fm, ,), So entsteht unter Hinzunahme von Fr symbolisch 


= | pean cee Pom 
U = Ping) 2 oS eer 3 (103) 


Die Randlésung Fz ist eine gediampfte Schwingung in Richtung m, die die Randbedin- 
gungen von (95) nicht stért und jene von (102) erfiillt. Man hat so die Méglichkeit, 
die Schnittgr6Ben auch in Richtung m von einem Rand zum anderen zu verfolgen. 

Wenn die Ausrechnung auch einigen Aufwand erfordert, so bietet sie nach den 
obigen Darlegungen keine grundsatzlichen Schwierigkeiten. 


VII. Aufbau einer vollstandigen Lésung unter Beriicksichtigung des Steigungs- 
parameters f mit Hilfe der Stérungsrechnung 


Setzt man (86) in (78) ein, so folgt aus dem Koeffizientenvergleich fiir /*: 


2 (02 Un 
a(S er 

1 2 (OF: 1 + lb sin iI (e] Cay 
+5 : 

r2 ag? l—w2rsiny r om 
I+ yp Vy a0 
1—p\omrog 
cosgsing 1 A@Vn-2 , l+mu1_ sing @Vn-1 
'l—-p r 2rsiny ray 


cos? gm 02 Un—2 3—pu cosm 62 Un-1 
4r2sin2y am? l1—yp 2rsiny dmrog 


— 


coszg 02 Vn_2 1+ cosp Vy_1 
472 sin? y Omrag l—yp 2rsiny rdg? 


2 : 
4r2sin?y r om 


(104a) 


74 cos?@ 1 A@Wn-2 2 ASW, , 2 cosm 1 @Wn-1 
l—yw4rsint?y ro am l—wr dm l—p2rsny r rd 
2 sing Whe 
; = 0. 
I—w2rsiny rr 
l+mu cosgp 8 Un-1 poe a2Un a ) 
l—yw 2rsiny am? l1—p Omrdg 
an Cos g d2Vn_2 3—" cCOsp a See 
‘am | 47r2sin2y am? l—yp 2rsiny dmrdag 
2  @V, , Qing 1 @Ve-1 
1—pwredag? ' 2rsiny r om E (104b) 
‘nurse cosg 1 A@Wn-1 2 ols OW yy 
'l—p2rsiny ram l—wr rep 


ey nein +, feos? @ 
= 1 13 2s SRC OE ie 


472 sin? y * 


| 
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» 7 
2u 1 aUp 
l1—yw r om a | 
4 2 cos p UO Vir oan 2 Ow, 
' 2rsny r om l—yw r rép 
We | 1 (7 Wn | 2 cos? p 04 W—2 cost~ 04 Wry_a 
© 6 [l—p\am " 4r2sin2y ams 16rtsinty amt 
4cosp Wn-1 cost 4 Wn_s & EW 
2rsiny am rag 2siny @meréeg Am? 72 dG 
coste 04 Wn=2 4cosqm 64 Wn-1 aa Ww 
| 4 B | E 
' 472 sinty am@rag? 2rsiny dmrag? | rags ‘ Nz ©) 
, 2sing 1 Wri , 6cosmsing 1 Ways 
2rsiny r ame ' 8rsiney r ams 
S5cospmsing 1 63 Whi , 6 sing 0 Wr-1 
2rsinry r amerdp | r 2rsny Omredy? 
2— yp costg—sint?g 1 02Wn-2 , 2 cosp 1 0? Wn-1 ! 
l—yw 4r2sin?y 72 r2am2 ' 1—yw Qrsiny r dmrag 
Die Wi, 
ip Pie ee 0. 


Die ersten beiden Gl. (104a) und (104) enthalten neben dem n-ten Glied die beiden 
vorhergehenden, die dritte sogar vier vorhergehende. 


Wir gehen nun so vor, da wir ansetzen 


ul 


CAUS Eel One Dre Ueie >, fs a Ae: Pf: Ue aos 


1 : 
fe an(G) 1 aG fp a 
= DEF = Gy. He roy @ eae ee (105) 


1 1 


1 
— Go {Ga PG 


1 
‘Danach haben wir aus (78) Pain zu bilden und erhalten, wenn wir die Glieder 
U,, V. und Wo, die aus dem vorigen Abschnitt bekannt sind, auf die rechte Seite als 
Lastglied bringen: 


2 a2 Uy, La Ga l+p a2Vy Dh, AL CONIC 
1—p om? | P ag? '1l—pomrdp 1l—-pr am 


3—u cosy 082 Uo alk yp sing 1 0Uo , 
=T_—pw2rsiny amrdp 1—p2rsiny r am 


(106a) 
Peet Nee 
ae 1—p 2rsiny COS P 72 ag? Py rag 
2, 1 ve AW» 5 Wo sin 2 p 
ST AG ay (cos ray Spat ane TP Gharsiny’ 


16* 
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LeU 0 Van 2 1 0M 

l1—p emrop am? l—p 7 dg? 1—pw rom 

l+p cosy 0Uo 3—p @Vo 2sing 1 dVo 

~T—p 2rsiny am fi l1—pamrog 2rsiny r am 
2 cosp 1 @Wo 


l—p2rsny r om‘ 


= >=> 2 = 


(106b) 


ou teu 2 120K P| 1 (ae eras ) 


l—wr am l—yw r rag 6 |1—p\ am 720g 
ihn 
Bi PS Mi beriney pss ise 
2cosy 1 dVo h2 1 4cosy { 04Wo a4+Wo 
— —s-« Qrginy + om ub 6 esd eee amirag | amrag ( (106c) 
2sing 1 Wo 6sing 1 03 Wo 
2rsiny r ams 2rsiny r amrag? 


2 cospm 1 maar 
l—p2rsiny r aomrdg)* 


In Bezug auf die Lastglieder der auBeren Belastung mu8 noch erwahnt werden, daf 
in den Gln. (104a) und (106a) nur ein solches fiir » = 1, also fiir /*=f vorhanden ist. 
In den Gln. (b) haben nur die geraden Potenzen Anteile, daher ist die Wurzel in (104b) 
aufgefiihrt. Man entwickelt sie gema8 (70). Die Gln. (c) haben nur fiir f= 0 einen Last- 
teil. 

Die Gln. (106) gleichen in ihren homogenen Bestandteilen den Gln. (87), wenn 
man fiir 

fe eB ew l+p eu 
l—p am am? 1—p dm 


setzt. Man gewinnt aus ihnen U,;, Vin, Win. 

Mit den Lastgliedern 2,;, 7,, und a3; von (106) hat man wie oben bei (78) eine 
Fourierentwicklung anzusetzen und erhalt unter Hinzunahme der homogenen Lésung 
analog (93) bzw. (101) 


U, =O 7 + Ui, + Up, 
Vi =Vii + Vin + Vir, (107) 
W,= Wir + Wirt Mie. 


In dieser Weise kann man den Rechnungsgang in Abhangigkeit von der Steigung f 
bis zur gewiinschten Genauigkeit fortsetzen. 


Zur Fourierentwicklung formt man Ausdriicke wie (91) 
A,e? cos 7,9, <A, e2? sin », 9, 
A,ze"? cos», 9, Age? gin », 

durch Addition und Subtraktion von 


3 A, e1* cos »,y und } A, e cos », p 
zu 


2 (A, + As) ©03 uw, pcos py + 4 (A,— As) Gin M1 P COS ¥1 
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um, da 


2 (eM? + e-?) = 603 wy y 
und 


} (eh? 0-9) = Gin, y. 
Es treten dann folgende Integrale auf 


“Ginax| . 
ee sin Bx cos yuda = 
ee es (L052) 
as a + (B+ y)? at + (B — 7)? 
[=o 
uae cos Bx cosya da = | 
Cos ax oe te 
oe 108b 
ae (108b) 
nit a Gina acos(B+y) x 4 a Sin « wz cos (B—y) 
a2 + (B+ y)? a2+t (B —y)? 
micawe e : 1 
oes sin Bx sn ya dx = 
aa 108 
0 (108¢) 
“|= _ «Gina x cos (6+y) x 4. % Gina cos (8 —y) x 
i a+ (B+)? a2 + (B— 7)? 


In gleicher Weise wie wir nach (86) U, V und W aufgebaut haben, geschieht dies 
mit den GréBen (58) und (61). Dann folgt beispielsweise fiir L144 nach (105) 


Dil rial ae | Cae i TAGE) a5 


Fiir den linearen Stérungsanteil hat man (58a) nach f zu differenzieren und dann 
f=0 zu setzen. Die hdheren Std6rpotenzen erhaélt man gemaf der angegebenen Ent- 
wicklung. 

Dieser Rechnungsgang ist fiir die anderen Schnittgr68en analog zu fiihren. 


(Hingegangen am 1. Juli. 1959.) 
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Compressibility Effects in Circumferential Inlet Distortion 
in Axial Compressors* 


By M. Z. v. Krzywoblocki, Urbana, Illinois 
Part I 
With 4 Figures 


Summary. There exist several approaches to various involved phenomena in the flows in axial 
compressors, like circumferential inlet distortion, rotating stall, etc. But these papers are concerned 
mainly with an incompressible fluid. In the present work the author proposes some similarity 
laws which allow one to extend the results obtained in the domain of an incompressible fluid to 
the domain of a compressible fluid. Due to the nature of the similarity laws, the results so obtained 
are approximate and refer to the subsonic flows only. In particular, in the present paper, which 
is the first one of the series, similarity laws are applied to Ehrich’s solution valid in the incom- 
pressible domain. 

In the Appendix (Chapter 5) the author outlines a procedure of obtaining a first, rough approxi- 
mation to the coefficient ds, eq. (4.2.5). 


Introduction 


The recent development of the dynamics of flow in axial compressors refers to such 
complicated phenomena like circumferential inlet distortion, rotating stall, ete. But 
the work done up today, refers mostly to an incompressible fluid. 

In most cases, the used equations are linearized and the results are approximate. 
Due to the linearization process, there seems to be quite natural the attempt to construct 
the similarity laws which may allow one to generalize the results so obtained to the 
domain of a compressible fluid. Such an attempt is presented in this paper. 

For illustrative purposes, the construction of the similarity laws is demonstrated 
on some particular methods derived in the domain of an incompressible fluid. 

In particular, in the present paper, the first in the series, the approach by Ehrich 
is chosen to be generalized to the compressible flow regime. The subsequent papers, 
to be published soon, refer to the methods of Yeh, Bragg and Hawthorne, and 
Ronnie and Marble. 

The last paper of the series will contain some experimental informations on the 
subject in question. 


1. Solutions in the Domain of the Incompressible Fluid 
1.1. Ehbrich’s Solution 


Ehrich! proposes the following solution of the circumferential inlet distortions in 
axial flow rotating machines: The circumferential inlet axial velocity component 
distortion for upstream of the rotor, u,, consists of the mean axial velocity, uw, and 
of the superimposed axial velocity with a uniform swirl angle characterized by a slope a,: 


U,=u-+ e, sin [n (y—a, x) r-], (Vaasy 


where ¢, denotes the amplitude of velocity distortion, 7 mean radius of passage and n 
the order of the harmonic component of the distortion. Since the actuator disc type of 


* The approach presented in the present paper was developed during the author’s association 
with the Compressor Aerodynamic Design and Development Unit, General Electric Company, 
Cincinnati, Ohio, U.S.A., with Mr. R. A. Novak. as its supervisor. The author wishes to express 
his deep appreciation to Mr. R. A. Novak for the numerous discussions on the subject in question. 

1 Ehrich, F.: Circumferential Inlet Distortion in Axial Flow Turbomachinery. J. Aero. 
Sci. 24, 413—417 (1957). 
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analysis is applied with the variations of parameters in the radial direction neglected, 
the consideration is reduced to a two-dimensional plane flow (cascade theory approach). 
The tangential velocity component is given by the formula: 


v,=a, {u+ e sin [n (y—a, x) r-]}. ; (1.1.2) 


The velocity far downstream of the rotor possesses a sinusoidal variation of a different 


amplitude and is shifted in circumferential direction with respect to the mean flow. 
Hence: 


U,=u-+t e, sin [nm (y—a, x) r3tq,l, (1.1.3) 
Vg= 4, {U+e,8in [n (y—a, x) r++q,]}, (1.1.4) 


where ~, denotes the circumferential orientation of the phase of distortion. The other 
equations, considered by Ehrich are: Two-dimensional continuity equation and 
the radial component of vorticity: 


Moot Uy We Yee tb, =e Gdn ln (1.1.5) 


To simplify the problem, Ehrich uses the linearization and the assumption that the 
vorticity is transportet along the lines of mean swirl; in the region upstream of the 
rotor (called region [1—2]), Ehrich assumes: 


V,x—Uyy=C (y—G, &); (1.1.6) 
and downstream (region 3—4): 


(1.1.7) 


Fig. 1. Nomenclature 


The general expressions for the velocity components satisfying all the foregoing equa- 
tions are chosen by Ehrich in the form: 
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In the region [1—2] (— co <a<0) (with 2 and 3 denoting the points immediately 
in front and after the blade row): 
u=u+e, sin [n (y—a, x) r-] +e, sin (nyr4+,) exp (nxr-); (1.1.8) 
v=a, {u+e,sin [n (y—a, x) r-]} + €, cos (ny r+.) exp (nxr-); (1.1.9) 
in the region [3—4] (0 <x <+ oo): 
u=u-+e, sin [n (y—a, x) r-1+,]+¢3,8in (nyr-1+ 3) exp(—nxr-*); (1.1.10) 
v=a, {u+e,sin [n (y—4ay, 2) r-1+9,]} —¢€3 cos (ny r-1+ 3) exp (—n xr-1). (1.1.11) 
The nomenclature is explained in Fig. 1. 
2. Similarity Laws for Ehrich’s Approach 
2.1. Fundamental Formulas 
In the case of a compressible flow, the continuity equation has the form: 
(0 &),x+(@%),»=0; or: [(@+ Qa) U],«+ [(o + ga) ¥],»=9, (2.1.1) 


where 0 denotes the mean density and g,the disturbed density. To find a similarity law 
in the present case, we introduce the streamfunction, y (x, y), satisfying identically 
the continuity equation: 


w= 2 Pry = OF Yaa Peo 
Equations (1.1.6) and (1.1.7) give: 
07" Osx Pra 057 Yi xe 0 Oey OS ig SG (2.1.3) 
From the equation of motion taken along a streamline, 
dp = — 0 d(4V*)+ o dP, with P,, = P, etc., (2.1.4) 


where P, denotes the x-component of the external force acting on the system. Using 
the expressions: 


V2=0-* [(y,s)? + (y,)7]; 0,2= C7? B23  c? =(eplédo)s; -(S= entropy); 


calculating the derivatives p,,, etc., from equation (2.1.4) and using them in equation 
(2.1.3) furnishes the result: 


(Cis 8) Prax t+ (L—v? -*) p57 — 24064 y,.,+ 9c YP, —UP,) = 
=—oC [1—(u?+ v?) c-?]. (2.1.5) 


Assume a running coordinate, s, along a streamline; the interval of s [— oo, 0] cor- 
responds to the interval [1—2]. At the position s=0 (= point 2) there is some heat 
addition Q (s), to the particle moving along the streamline s, where Q (s) denotes the 
total heat amount added to a particle from s = 0 to s. Similarly, the work of the external 
force along the streamline on the distance [0, s] is: 


Nee \de = P'G)n PP (2.1.6) 


Let the bars over the symbols denote the properties of the medium in the mean, un- 
disturbed flow; the symbols without the bars refer to the regions in front of and after 
the blade row. Then the energy equation has the form: 


& (wu? +o) + (y — 1) c? = 3 (wu? + 0) + (y— 1)? —Q (s)— P(s), (2.1.7) 
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with “c” denoting the velocity of sound. Application of the linearization procedure, 
similar to the classical one, gives an approximate result: 


c# = eq? ~ 6? + (y—1) (Qm+ Pn); (U2 +02 —u2— v2) ~ 0, (2.1.8) 


where Qn and Pn denote some constant average values of Q(s) and P(s), referring 
to the interval in question; these values of @Q,, (s) and P,, () must be taken from some 


_. experimental data and are considered to be parameters of the system in question. 


Then we may put down: 


uc? w uc,” = M2; v2 c-2 w 6, = M,?. (2.1.9), (2.1.10) 


As the next step, we apply the classical linearization procedure to equation (2.1.5). 
The justification of this procedure applied to the stream function equation is the fol- 
lowing: Fundamentally, the stream function y can be decomposed into two parts: The 
stream function y referring to the undisturbed uniform incoming flow of the velocity 
components u and ¥; the stream function wa referring to the disturbed flow and disturbed 
velocity components. Similarly, we can assume: u=u+uz; v=v+v, ete. 

Assume that in the incoming undisturbed flow there are no forces acting and the 
vorticity is equal to zero. Then equation (2.1.5) may be decomposed into two parts: 
One referring only to the incoming undisturbed flow: 


(1 —@? c~*) Prax + (1—0? —2) p,yy—2 UVC Pyzy = 0. (2.1.11) 


The remaining part of equation (2.1.5), which is left when the expression (2.1.11) 
(equal to zero) is inserted into equation (2.1.5) will refer to the “disturbed” stream 
function, yz. During the process of linearization we may neglect the terms containing 
the components of the force if the force is small. The influence of the external force 
and heat is taken into account in the expression for the velocity of sound. In the case 
the external force is large, it should be taken into account in both the incompressible 
and compressible flows. The result of the linearization is: 


(1 —M,?) pase + (1 — M,?) Pasyy = — Om f [(1—(M? + M,*)]. (2.1.12) 


Apply a transiormation of coordinates from the (2, y)-plane into the (a’, y’)-plane, 
with 2,, A), Ay, 4c, denoting the constants: 


g=Aas Y=HAy; CHKCS vd @TY)HA ve (2, 9)- (2.1.13) 
With: 

1—M,=£.2; 1—M=8.?, (2.1.14) 
equation (2.1.12) furnishes: the formula: 


Bu? Ae? Ag hye er + Bo? Ay? At Vay y = — Om Ae!’ (Bk + BA 1). (21.18) 


The value of the density, 9, is calculated from the equation of energy with: 


p p= (0 4)" exp [(S—S) ¢,']; (2.1.16) 

po-t=po- (9e@-1)"-» exp [(S—S) ¢,"]; (2.1.17) 
e=olv—)y7k @+e—wv—w)+ 

+@Q (s) +P (s)]ep4t1} exp [(S—8) ]]"™. (2.1.18) 


The mean value 9,, in the interval [1—2] is assumed to be equal to: 


Om (1—2) = Om1- (2.1.19) 
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The mean value in the interval [3—4] will be calculated in section 2.2. 


We require that all the three coefficients of yy. es Wy; and C’ in eq. (2.1.15) 
be equal one to each other, thus reducing this equation to the form in the incompres- 
sible flow: 


But A= 223 Bo® On! Bu? Ax? = Ay Ae! (Bi? + Bo? — 1); (2.1.20) 
On’ Bo® Ay? = Ay Ag! (Bu? + Bo? — 1). (2.1.21) | 

With /,=f, we have: 
Ay = Bus (Bu? Bo®) Ag= Om Ay (Bu? + Bo? — 1). (2.1.22) 


Assume that the plane (2’, y’) refers to the domain of the incompressible fluid flow, 
and the plane (2, y) to the compressible fluid flow. Then from the definition: 


Wire = — 0a! = Ay Ae pare = — Ay By! 004% — Bu Me. (2.1.23) 
This seems to justify the assumption 


Ag 07" Br Prete eis fix 


Thus the correlation formulas are: 
pa’ (2", ¥')= On Bu Bo pa (@, Y); Vd =— Yaw © On BuQ Va © Buta; (2.1.24) 
Ud =Yary = Ay Ay Yay © Om Bo Q Ud © Bota (2.1.25) 


Assume that the rules analogous to the correlation formulae (2.1.24) and (2.1.25) are 
valid for the incoming flows, with y, uw, etc., referring to the uniform incoming flow of 
the compressible fluid and wn, Un, etc., referring to the uniform incoming flow of 
the incompressible fluid. Thus, we propose the relations analogous to relations (2.1.24) 
and (2:41.25): 


Yin (x' y')= Ay P(x, y)s Y (2, Y) & 6 (Bu Bo) Vin (2", Y’); (2.1.26) 


hy = 07 Bu Bos te = Bri = Ba. (2.1.30) 
Using these values in equation (2.1.11) and keeping in mind that in an incompressible 


flow, the velocity of sound approaches the infinity and consequently the last term in 
equation (2.1.11) tends to zero, one obtains: 


Bi2 a2 Wins x! x + Bi? Ay? Vinsy’y’ ae 0; Vy? Yin = 0, (2.1.31) 


which is exactly the stream-function equation (Laplace equation) in an incompressible 
irrotational flow. Hence the correlation formulas between the incoming velocity com- 
ponents in compressible and incompressible fluid flows, respectively, are given by 
equations (2.1.27) and (2.1.28). 

Thus we have: Compressible fluid: 


U=uUutug; v= vty; (2.1.32) 
Incompressible fluid: 


UW =Uin + Uj; V' = Vin + 04; (2.1.33) 
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Correlation formulas: 


v = fz Wins Ua © 0 Om Br? ua’ © Bol ud’; (2.1.34) 
=e he Vins Va © O-) Om Bri va’ & |e dped OF ae (2.1.35) 
U=U + Ug © Bz) Uin+ O-! Om Bo? Ua’ Bs? Win + By? ud’ ; (2.1.36) 
V= V+ © Be Vint O-} Om Bet Vd © BO Vin + Be vd. (2.1.37) 


To avoid lengthy notations, we shall use the following combined symbols: 


| 


u= Bo) Unt BF! ud’ = pO w’; (2.1.38) 

V= Bz) Vin t+ Bo va = BO V'. (2.1.39) 

This enables us to calculate the first appreuhnntion to the density from eq. (2.1.18): 
0=0 [{—1) 73 B40 flow? ploy’) + 

WAS (8) e po 1} expy = Oc, 1) | es), (2.1.40) 


with S=0. The pressure can be calculated from eq. (2.1.16), temperature from the 
equation of state. We may find the relation between the tangents of absolute flows 
angles in both domains: 


Oye qty = poe Vy’ =a, ee Uy’; (2.1.41) 
or es 
0, = 0% (BT Un p, 4 ey to2) =a," (By tn + Beye); (2.1.42) 


which furnishes the approximate value: 
apap ay (2.1.42a) 
A similar reasoning applies to all a,’ ((=1, 2, 3, 4). 
Equations (1.1.8) to (1.1.11) will be used in the forms: 
os <2 =< 0: 


w= BO) {tin + €1' sin [n (y’ —ay' 2’) r'-1]} +655 eo sin (n y’ 71+ y') exp-(n a’ r’-1) ; 


(2.1.43) 
8 = By Ay’ {Un + By" sin [nm (y’ — ay! x’) roast +B Eo" cos (n ap eae + 2’) exp (n a’ Pos 
; (2.1.44) 
0<%<-+ co: 
w= BO in + €,' sin [n (y’ —a,' 2’) rato, }}.+ 
+ Bo} 64! sin (ny! 7-1 +95’) exp (— nz’ r'-1); (2.1.45) 
y= PO a," {Uint eq sin [2 (y' —a,' x’) r'4+9,'}}— 
we: Bo oy cos (n-y’ y-l +3") exp (—nz’ ys). (2.1.46) 


$29 LHS os : / oe 
One more item,may be of some interest, namely the relation between the: densities 
in both domains. With the symbols {Uin, Vin, Pin, Qin} denoting the fluid flow properties 
in the incompressible, incoming undisturbed flow and {u’, v', p’} denoting the total 
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velocity components and pressure in the domain of the incompressible disturbed flow, 
one obtains from the Bernoulli equation: 

4b (u? +02) + Din 0; =F (wu? +07) +P’ 0; (2.1.47) 


This may serve to calculate the pressure p’ if all the other characteristic properties 
are known. If one assumes that the “Bernoulli constant” is the same for both com- 
pressible and incompressible flows, one obtains: 


4 (U2 +05) + Pin Oy =E(W+0) +7 (y—1)*pet=r(y—1)* pei, (2.1.48) 


with the subscript “‘s’’ denoting the stagnation conditions. Hence, with 


uU — B Uin' v = fis Vin} 7 oo a (o ge Ds Re: 


the relation between the density in incoming compressible fluid flow, e, and the 
density in the incoming incompressible fluid flow, 9;,, is given by the relation 


0= 0. {(y—1) 74 FE (LB) 0, +4 (1B?) ©, + Bin Oi!) D,* Qf. (2.1.49) 
2.2. Matching Conditions for the Axial Velocity and Rotor Exit Angle 


Assuming that the annulus area is constant, we have the equality of the entering 
and leaving axial fluxes: 


Oo U2 = 02 U (—9, Y) = 03 Us = 034 (0+, y). (2.2.1) 
Using eq. (2.1.38) one obtains: 
Om BO" ue! (— 0, Y) = Oma BO-* ug’ (0+, 9). (2.2.2) 


Use of eqs. (2.1.43) and (2.1.45) furnishes the formula: 
Omy {B{~* [tin + er’ sin (ny’ r’-1)] + B53 e9' sin (ny! r’-? + ,')} = 
= Oma {BQ-"[Uin+ €4' sin (ny’ 1-1 +4')] + By} és’ sin (ny’ r’-1+¢,')}. (2.2.3) 


Decomposing this equation into the matching phase components of cos (n y’ r’-), 
sin (ny’ r’-1) and the terms free from trigonometric functions, we obtain the value 
for Qn, and two matching conditions: 


Oma = Br ns Omi (2.2.4) 
[B53 + Boe (€2’ ef?) cos y,"] = 
= Bos Bri [Bos (ea €{*) + €08 ya’ + B55 (€9' £17) Cos G']; (2.2.5) 


Bo (€2' €¢*) sin @o’= Bot Bos [B,} (Ea’ €7*) Sin ya’ + Bos (e9' 6(1) sin Ps |. (2.2.6) 
It seems justified to assume: 
MM? w M?, ~ ur c-*; M*, » M?, wwe”, (2.2.7) 


Bor wx Boa; Bos ms But: (2.2.8) 


Then the conditions (2.2.5) and (2.2.6) reduce to those given by Ehrich. 

With the assumption that the relative exit angle of the medium leaving the rotor 
coincides with the geometric exit blade angle regardless of the relative angle of incidence, 
one obtains: 


U3 6,=v3s+U; U=wheel speed, (2.2.9), 
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with 6,=tan B, denoting the tangent of the relative angle of flow. Using equations 


(2.1.45) and (2.1.46) and the technique explained above, furnishes two matching 
conditions and a condition for Uy-: 


bs [By (4 enn) COS M4’ + Bo; (€3' €{ *) cos Ps | = 


= Bua Ga’ (64' (1) cos p,! + B73 (es’ ef) sin py’; (2.2.10) 

bs [Bri (ea ef *) sim py’ + Bo} (e5’ ef) sin gy'] = 
= By, %4" (e4' ef) sin py’ — Bo (eq! &/-) cos Gy’; (2:34 1) 
U ui) = Bo Os —Byi Ma - (2.2.12) 


2.3. Energy Condition 


Assume that in the interval [1—2] there are no losses, no heat addition and no 
work done along any streamline: 


2 (us? + 017) + (y— 1)? py OF =} (ue® + 022) + (y—1) py 031; (2.8.1) 


inside the rotor there takes place the heat addition and the work of external forces: 
hence: 


2 [Wo?+ (ve+ U)?]+y (y — ie Po 03! = 
=2 [Us + (Vs + UP] + (y 1) Ds 0st —Qz2-s (8) — Po-a (s); (2.3.2) 
no losses, no heat addition and no work of the external force in the downstream region: 
2 (Us +s) +y (y—1) Ds Ost =F (uP? +00) +7 (y—1) pay. (2.3.3) 
Combining these three equations furnishes the formula: 
uy? (1+ a,?) —u,? (L442) + 2 y (y—1)*? (py 14 — Da Oy") = 
= 2 U (v3—v.) — 2 [Q (s) + P (s)]o-s. (2.3.4) 


The following values are used: 


Uy = By Win + Boy Uys (2.3.5) 
U4 = Boa Min Boy Wu; (2.3.6) 
Vi=Byj Vin + Bi; %; (0 = 2,3). (2.3.7) 


In using equations (2.1.43) and (2.1.45) and the values of uw, and u,, the following 
assumptions are followed: 


The symbol w, denotes the magnitude of the velocity very far upstream of the blade 
row; because the exponential function is the fastest convergent function, the term 
containing the exponential function in equation (2.1.43) is neglected ; for the same reason 
the analogous term in equation (2.1.45) is neglected. 
One term retained in eq. (2.1.43) and one term retained in eq. (2.1.45) contain the 
coordinate x’, considered to be a parameter. Practically, two magnitudes of this para- 
meter must be chosen in each case: The terms containing the coefficient (1-+ a,*) refer 
to uw, and to the negative value of x’ upstream (— oo <x <0); the terms containing 
the coefficient (1+ a4?) refer to u, and to the positive value of x’ downstream (0 <x < 
pase 

Inserting the similarity laws into eq. (2.3.4), using the general expressions for the 
velocity components in the domain of the incompressible fluid flow (u;’, v,’), given 
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by eremons (2.1.43) to (2.1.46), matching the phase components of sin (n y’1r’) 
and cos (n y’ r’-!) and ihe mean value terms, and neglecting the terms of order (é’*), 
we obtain the set of three conditions*: 


ne os (1+ a,?) cos (n a,’ x’ r’-*) = BA Be ay se: ) [cos 4’ cos (na,’ #’ r’-+) + 

+sin (na,’ x r’-1) sin y,’] (1+ a4?) = U uj) [aa' Bis (eq' €¢*) COS Yq’ + 

+ Bz} (ea! ef) sim ga! + B;3 (e9’ €(-1) sin ga’ — Bt ay’); (2.3.8) 
eRe Bot (1+ a,?) sin (n a,’ x’ 7’) Soap BS (ex f7) [sin py,’ cos (na,’ x’ r’-1)— 


—sin (n a,’ x' r’-1) cos p,'] (1+4,7)=U u;" [ay’ BF) (e4' €7*) sin v4’ — 


— B73 (é3' €; ) COS ps’ — B 5 (E2' €, |) COS po’); (2.3.9) 
where: | 
ay Oy” BL Boys Gal Oy" Bos Beas (2.3.9) 


EL (Bit + But ar) — (Bos +Bus 4°) ] Mint p 
9 iss Lo Eps O1 —Pr01 )=U Un (But a,’ —But ay’) — 
~[Q (s)+ P (8) 2-2. (2.3.10) 
2.4 Set of Matching Conditions 
The six matching conditions (2.2.5), (2.2.6), (2.2.10), (2.2.11), (2.3.8) and (2.3.9) 
allow one to calculate the six unknown variables: 
(25':6, *) Sin pes , (€q' 67 7) COB Me’; (eg) 24 1). sin.ps 5 
(254, )' COS Gas (6, 2.57) -singy. ((es ec, | Cong, 
Following Ehrich, one can use a matrix representation for the six equations in ques- 


tion. In this, each magnitude of independent variable x is considered to be a fixed 
parameter. One obtains: 


Ay On = O43 90 yaeGyam © 0 

0 aj; Oa geo aaa =—Ba 

0 0 G33 —b33 0 —Dgy Meera 0 (2.4.1) 
0 0 bss gg—— Ogg OO 0 

0 —D,. 0 —Dgg Giza —Dig W 

bse 0 63; 9 Bsa Usa x 


where the column matrix ¢ is: 


Sat fs 0 cos 2! 


* There is a misprint in Ehrich’s paper in eq. (17a). It should be: 
(1 + a?) — (1 + aj?) (e1/€1) cos gy = 
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and. the coefficients have the values: 


@y2=B,2; G13=B)1 Boa Bos; G14=f)) Boa Bois (2.4.3) 
@s3=Bi3; bs3=b38,3; bs4=b B;1— a,’ B}; (2.4.4) 
Asa = Br} Boe (1+-447); sg =e By; (1+a,°); (2.4.5) 


Asa =(By4 Ds Brida’) Bri Gy’ +B,4 By i Cos (nay! a’ r’-1) (lt az); (2.4.6) 
bs2=(Br1 bs —By 4 Gy’) Bua; bss =(Br4 b3—Bui 4’) Bis: 
bsa=fy1 Bi sin (ma,’ ar’) - (1 4+-a,2); 

W = — B71 Be! (1+.4,2) sin (n ay’ x’ 1’); (2.4.8) 
X =By1 Bos (1 +442) cos (may! x 1’) 4+ Bot ay’ (Boibs—Bui Ga’). (2.4.8 a) 


After some elementary operations, the above system reduces to the system of two 
equations: 


Ay; (Eq &'-") sin pg’ + Aye (€4' €1*) CoS 4’ = (Way, — Boi yg) O ge D595 (2.4.9) 
oe (og &. SI gy = As (ey & * ).CO8 4" =X Oy, Oss Usa (2. 
A, = (33 + O33) (— G2 Usa) — (Are O53 + 213 O59) Aas Oa; (2. 
A yg = Digg Dag (= Gyo Ggy) — Oy 4 Beg 054 O55 — Gre Ugg Ogg Ona; (2.4.1 : 
Ao; = — (a3, + 35) (Aya D5o + Gye O54) + (Gyo O53 + O13 O59) B35 O34; (2, 
Aaa = (Aga Ucg 1719 Osa) ag 09a + (Ae O55 + O15 Ose) De, +01 Gea Osa Ogg { 
This furnishes: 
tan yg,’ = [(— W ay4.+ Boi Giga )eAns Ug aes Ags Oia-Cag Onn 
> [(W Gy. — Bot G33) Aas O54 O59 —X Ay Gye G33 Osa]; (2.4.15) 
(ex R= = Walaa Bat G35) Aes Daa Osa ct X Ans Gia Ose Osa)? + 
+ [— Way.+ >) 43s) Ao: Osa O52 +X Ary G12 Mga bsal}?} (Aye Asx — Ary Age)@?. (2.4.16) 


Assume the flow through a stator (stationary row of blades); then the wheel speed 
U =0,,or Re bs— oe Lae Oy 


and, as one can easily verify: 


b= 05,0; is, = By! Bo, (1 + a4") cos (n ay’ x’ r'-1); (2.4.17) 
X =f} B;} (1+.4,%) cos (m ay! 2’ 7-1). (2.4.18) 
This gives with B3,=(b3 6);—44' Bi 3) © 9: 
tan ~4/\y-0=9, (2.4.19) 
in agreement with Ehrich’s result. Similarly: 
Ay, = — (a2, +535) 124543 Arz= — (G54 433 + M53 O54) Gre O34; (2.4.20) 
Ag = — (2, + 625) Gyo 0543 Asz—= — (054 O33 — 54 V5) G12 O34; (2.4.21) 


= Bot Bot (1 +42) cos (ray! 2’ 7’), (2.4.22) 
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and after division by 4, b,: 


(e4' ef?) = [(X Aye @g3)? + (X Ay Gg5)?] * [(A54 O33 +438 954) (433 + B33) 12 O54 — 
— (b54 533 — G33 Asa) (a3 a b:) 4 A54]*. (2.4.23) 


For: Bor Bra ~ Bos © Boa © lL, a ow,9, 


this result agrees fully with that of Ehrich. 

Equations (2.4.15) and (2.4.16) give the magnitude of the phase shift and of the 
attenuation of a distortion through an isolated rotor. Similarly, equations (2.4.19) 
and (2.4.23) give these magnitudes for a flow through a stator. The equations and 
the procedure, presented above, is directly applicable to a cylindrical lamina of a 
radius “‘r’’, under the assumption that there is no radial velocity component, no vari- 
ation in the radial direction and that the radius itself is considered to be a parameter 
independent of the differentiation and integration operations. In this case, the equa- 
tions of continuity and vorticity have-the forms: 


t*(0 Cy)ig + (6. Cea 0: SO Cgc: (2.4.24) 


These equations are reducable directly to equations (1.1.5) if one uses the substitutions: 
C,=v; r6=y. Hence, in this case, the method, presented above, is applicable 
directly. 


3. Second Approach 
3.1. Fundamental Assumptions 


One can propose another approach to the solution of the problem in question. 
Assume steady state flow conditions and the validity of Ehrlich’s solution for the 
velocity components, eqs. (1.1.8) to (1.1.11) and write the continuity equation in the 
form: 

(u Z2+-¥ 2,5) +(a,.-+2,,)= 0; *z2—l 0. (i 
The second paranthesis is equal identically to zero, since Ehrlich’s solution satisfies 


identically the equality div 7=0. The first parenthesis in eq. (3.1.1) is identically 
equal to zero, if the following relation is satisfied: 


%x (%y) 1 =—vUu-, (3.1.2) 
This is true if: 


2=U (y—a, «)—rn- e, cos [n (y—a, x) r-1]— 7 n— £4 cos (r-1 ny + qe) exp (n x 1r-). 


(3.1.3) 
Since @=exp Zz, one can write for — oo <a <0: 
0 = 0. exp {A [u (y—a, x) —1r n-1 e, cos (nm (y—a, 2) rt) — 
—1r n~ &, COs (7-1 n y+ Hz) exp (nr! 2)]}, (3.1.4) 
and for 0 <2 <4 oo: 
0 = Ooo CXP {A [u (y—a, x) —rn e, cos (n (y—% 3) Pz) 1.4) ti 
—rn- €, cos (n1-1 y+ M3) exp (—nr- a)]}, (3.1.5) 


where the coefficient A takes care of the proper dimensions of the expression (it may 
be equal to unity). The exponential functions are the fastest convergent functions. 
Moreover, they converge for all finite values of x or y. Hence, we shall restrict ourselves 
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to finite values of a and y, and apply the expansion exp x=1+ a+ sre +5.., 


preserving only the first terms (in the first term containing the y-coordinate we 
preserve only the term 1): 


for —co<a# <0: 
0 = Qo exp (— A wa, x) {1—Arn- ze, cos [n (y—a, x) r] — 
—Arn- e, cos (r+ ny+ qs) exp (nr- a)}; (3.1.6) 
For. 0 <7 <4 co: 
C= 0. exp (— A ua, x) {1—Arn-'e, cos [n (y—a, 2) r+¢,]— 
—Arn- ce, cos (r-- ny +3) exp (—nr-} a)}. (3.1.7) 


For an incompressible fluid flow one has to put down 4=0. 


3.2. Matching Conditions 


Using the above values of the density in the first matching condition 9, w,= 05 Us, 
preserving only the terms of order (¢), and matching phase components of sin (n r-! y) 
and cos (nr-1 y) furnish the first two conditions: 


1 + (€, e{') cosy, +A rn u (é, ef) sin py = (€4 E11) COS gp, + 


+ (és e1!) cosys+ Arn u [(€5 71) sin pst (Ee, E72) SiN Ya]; (3,221) 
(ep €;1) sing, — A rn u [1+ (eq €,~1) Cos Hy] = (€4 €72) Sin gy + 
+ (e3 €71) sng,;— Arn! u [(€5 71) cos pg (€4 €12) COS yg], (322.2) 


with for A =O reduce to Ehrich’s conditions. 
The two other conditions are purely geometrical and are given by the rotor exit 
angle: 
u; b,=0,+ U. (eq. 2.2.9) 


Since in the present case, the velocity components are identical to those in Ehrich’s 
approach, these conditions remain the same as in the incompressible fluid flow: 
bs (€4 €11) COS My +)z (Es e171) COS Pg = Ay (Eq E11) COS Wa + (Es €7') SIN YQ; (322-3) 
bs (€4 €12) SIN Pa +)z (Eg €1') SIN Y3 = Aq (Eq €1') SIM Py — (Ee E71") COS Ps. (3.2.4) 
Inserting the general expressions for the velocity components into the energy equation 
(2.3.4), eliminating terms of order (e?), and matching phase components of sin (n 7 y) 
and cos (nr-! y), one derives the last set of conditions: 
[cos (r-1 a, x) + (e, e11) exp (n x 7-1) Cos ye] (1+ a”) — 
— [cos (r-1 a4 X) (€4 e711) cos pg sin (7-1 a, x) (€, €11) Sin py + 
+ (es €71) exp (—n 1-1 a) cos ps] (1-42) =U U4 [(€, e7?) sin yy + 
+ (€g é;1) Sin Pg + Gq (Eq €1') COS Py — 44] ; (3.2.5) 
[(eo 12) exp (n x r-1) sin y,—sin (r-* a, x)] (1 + 4,7) — 
— [(e3 611) exp (—nr- a) sin p+ (€, €1’) Sin Py COS (ra, %)— 
— (€4 €51) 008 gy sin (7-1 ay 2)] (1-042) =U w [— (e2 13) cos p— 
— (€3 &11) COS py t (Eq €1") Gq SiN Ya]; (3.2.6) 
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uw? [ (L-+a,%) — (1+44?)] +2 (y—1)* y (Pi @1*— Pa 2) = 
=2U (a,—a,) u— 2 [Q (s) + P (s)]e-s. (3.2.7) 
As above, one has to choose two values for x: for — co <a <0 [terms containing 


(1-++a,2)] one has to choose the negative value of x, and for 0 <% <~+ oo [terms 
containing (1-+a,”)] one has to choose the positive value of 2. 


3.3. Matrix 
The matrix has the form: 
li Ais 1 —,, —1 — Ay. —l —l 
1 es Ofatet sexe Cig peasy) Ay2 G12 
pir Ess aes, 0 
0 1 OB 0 (bs — a4) i Fee - (3.3.1) 
0 0 bs 1 (bs — a4) 0 0 
A52 bsz 52 bss Mas bsa x 
— bse 52 — bs3 O52 — bsa 54 ba 
with: 
Qyg=Arn U; dse=(b,—4,4); ~Gsa—=(1+2,7) sin (ra, @); 
bs2 = — (1+a,?) exp (nar); bss =(1+4,?) exp (—n xr); (3.3.2) 
bsa = (1+ @4”) cos (7 Gq &) + (03 — Gy) 4; 
X = [eos (r- ay x)] (1 +442) + (03-44) (3.3.3) 
Y = (sin ir a, 2) (1) + @,*). (3.3.4) 


Multiply the second row by (— 4a.) and add the first row to the second row; similarly, 
multiply the first row by a,, and add the second row to the first row; then the first 
two rows divided (1-+ a?,,) are identical to the first two rows in the Ehrich’s case: 


Tet eer On 0 
ee 0d Seen teeese Ieee es 


Since the factors a,, or A (which take the compressibility effects into account) do not - 
appear in any other row, this implies that the compressibility phenomena do not 
have any influence upon the attenuation, if the Ehrich approach is used. 

The fina] result is: 


Ay, (&4 1") Sin Py + Aye (Eq 71) COS Py = (X + bra) (1 +857); (3.3.6) 
Ag, (€4 71) Sin Mg + Age (E4 €71) COS Mg =(Y + ase) (1+5,2), (3.3.7) 
where : 
Ay = (1 +55?) (as2-+ sa) — (2 d52 b3 + 52+ dss) (bg — 4); (3.3.8) 
Ay. = [2 ds2 (bs — 4) +052 + Dsa] (1 O51) (2 452 bs + bs2 + bs3) (bs — 4) bg; (3.3.9) 
Ay = [2 As2— bg (bs2 + bsz) (63 — @4)] — (bs2 + bs) (1 +052); (3.3.10) 
A y9= [2 As2— bs—(b52 + bss) ] (bg — G4) bg + [52+ Osa — (bg a4) (b52 +053) ] (1 bia ee ) 
3.3.11 


The consecutive procedure is analoguos to that given in Chapter 2. 
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4. Third Approach 
4.1. Fundamental Assumptions 


Assume that for — oo < x =a ()e 


eS eae sees 
uU=uU+ = é1nSin [nm (y—m, 2) r-1] + »’ eo, sin [n yr + on |exp(nxr-1); (4.1.1) 
n= n=1 
ae a : 
v= 4,4 ote ae Wy Ein SIN [n (Y¥—Qn x) r We 

n=1 
oa au an COS [n yr-4++qo,] exp (n xr-); (4.1.2) 

n=l 


and for 0 <a%<-+o: 


“=uy+ 2 fap BIN [n (Y¥ —Qan x) tani 


+ >’ e3, sin [n y r-1+ 93,] exp (—n x r-3); (4.1.3) 


n=1 


V=G,U+ 2L an Ean 8M [N (Y — Aan X%) 71+ Yan] — 
ad 


— >’ és, cos [n y r-1-o3,] exp (—nxr-). (4.1.4) 


The expression for the density is obtained in the manner analagous to that one used 
_ in Sectiom 3: 


fOD —cot< a= 0: 
0=0..[exp(— A ua, &)] 4 r >) n-} Ein COS [M (Y— Qin %) 7-4] — 


n=1 


co 
eH! 


EWG he 


n= 


n— &, Cos [rn ¥+ Gon} exp (n Zale (4.1.5) 


for 0<a%<-+oco: 


v = 0,. [exp (— A ua, @)] f —A rom €4n COS [1% (Y—Gan X) 11+ Yan] — 


—Ar >) &3nn— cos [r-1 2 y+ Pan] exp (—nx ro} : (4.1.6) 
n=1 


4.2. Matching Conditions 


Comparing the. terms containing the factors sin (n yr~*) and cos(n yr) in the 
equation of continuity furnishes the infinite set of the following pairs of matching 


conditions (one pair for each n=1, 2, 3, ...): 
. 5m = 2 =4 sth te 
(€2n €,1) SIN Yan —A TM —) U (Ean €,1) COS Yan— (€3n €j,) SIN Pan + 
+A ru n- (én 1) COS Pan— (Ean €j,) SID Pan + 


+Arun-! (an &,) COS Pan=Arun; (4.2.1) 
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Aru n- (én €1) SiN Pant (Ean &],) COS Pan—A UT N™ (Ean 2;)) sit Gan 
— (€3n €;1) COS Pan—A UPN (E4n €y,) SIN Pan — 
— (€4n €,1) COS Panz=—1, (n=1, 2, 3, ...): (4.2.2) 
The next condition, used to furnish the two subsequent matching conditions, the 
rotor exit angle condition, will be expressed in a different form than that given in 


previous sections. To allow for the variations of the ratio (v3+ U) uz!, we shall use 
instead of eq. (2.2.9) the formula: 


us tan (Bs + 63) = U3 (63 +43) (1 —bsd3)-*=05+ U; (4.2.3) 
tan’ B,=0,;; tan-o,=—d;: (4.2.4) 
where: 
dg = >) An éan [sin (ny7r-!)] cos Ysn+ >” Bn Egn [Cos (n yr-1)] SiN Pan + 
n= n=1 
+ >) On Egn [sin (n yr )] sin Yen D> Dry &gn [cos (ny 7r-1)] COS Yan; (4.2.5) 
n=1 n=1 


where the coefficients A,, B,, C,, D,, unknown for the time being, are constants or 
parameters depending possibly on (x, y). A comparison of the free terms and terms 
containing sin(n yr-1) and cos(nyr-!) in eq. (4.2.3) furnishes equation ub,= 
=a,u+U, and an infinite set of the following pairs of matching conditions (one 
pair for each w= 1, 2, /3,"...): 


(€3n €;.) (SiN Yan) (UC, +a, ub; C,+ U bs C,—1)+ 
+ (€3n €,,) (COS G3n) (4 Anta,ub, A,+U bs A,+63)+ 
+ (€4n €,) (COS Pan) (bs — Gan) = 0; (4.2.6) 


(Ean €,) (SiN Pgn) (U B, +a, ub; B,+U b, B,+b3)+ 
+ (én €;,) (COS Yan) (& D, +a, ub, D,+U bs D, +1) + 
t+ (Giq 2) (Bint ay) (0g 0a) 0, Cv Ske oe (4.2.7) 


The energy equation furnishes and infinite set of the following pairs of matching 
conditions: 
U u- (én €;,) SiN Yan — [exp (n w7-)] (1 +a, Gin) (€2n €1) + COS Pan + 
+ U u-! (ean €;,) SiN Yan+ [exp (—n x r-1)] (1+ G4 Gan) * (€3n €1) COS Pan + 
+ (1+ ay Gan) (Ean €[,) SiN (M Gan % 7-1) SiN Yan + 
+ [U U4 Qan+ (1+ 4 Gan) COS (0 Gan £ 1-1)] (Ean €;') COS Pan= 
=Uutant (l+4, din) cos (nai, x 7-1); (4.2.8) 


lexp (n «7)] (1+ a1») (Ean €,) SIN Pan + U U (Ean €,) COS Pan— 
— [exp (—m wr*)] (1+ aq Gan) (€3n €;)) SiN Yan + U U- (én €;1) COS Pan — 
— [(U U- dan) + (1+ G4 Gan) COS (0 Gan © T-1)] (Ean E71) SiN Pant 
+ (1+ @q Gan) (€4n €],) SiN (M Gan ZN) COS Pan = 
= (1+ 4 a) sin (1,212), (n=1,2, 8, .. 2). (4.2.9) 
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4.3. Matrix 
The matrix equations have the form (foreach n= 112 301).A): 
Aas 1 = has —] — Ars al w/a oan Pan 
1 4, a B19 | Ay. Ay2 
0 0 (P67, C,)) (05 — 627A.) 0 —(b3—dan) on 0 
0 0 (63 +535 B,) (1 +633 D,) (53 — Gan) 0 | ek Oak 
52 bse 52 bs3 sa Dsa x 
— bso 52 — bss A52 — Dsa O54 ¥ 
(4.2.10) 


where the used symbols denote: 


A,=Arun; ds=UU; asa= (1+a4a4,) sin (n a4, x7-1);  (4.2.10a) 


bs3= (U+a,Ub3+U bs);  bsa= — (1+ a, ain) exp (n wr); (4.2.10b) 
bsg (14-4 Qin) exp(—n ar); bsa=U U-) dant (1 +44 Gan) COS (n Ga, x 1-1); 
(4.2.10c) 
X=UO0 utarnt (1+, ain) cos (n 1, x7); (4.2.10d) 
Y == (1-1 g7ay,) sin: (way re) (4.2.10e) 


Notice the similarity between the matrices in the second and third approaches. 
After some manipulations eq. (4.2.10) takes the form: 


ae (63 — an) (1 —bg33 Cn) + Ko (As2 + Asa) (1 — 33 C,)} (€4n Cig) SIN Pan + 
+{ — Ky (bs — dan) (bs +33 Bu) + Ke [2 (63 — Gan) As2 + (b52 + 54) (1 — B53 Cn) ]} ° 
* (€4n 1) COS Pan= (X +52) (1—b35 C,) Ko; (4.2.11) 


HS (bs — an) (1 —bs35 ©.) awe K, (bs2 + bsa) (1 ee bss C,)} (€4n Sy) sin Pan+ 
41K (b3—4an) (bs +033 Bn) + Ky [—(63— Gan) (52 + 053) + (G52 + M54) (1— Bg Cn) ]} 
: (€4n é,,) COS Pan = (Y + as2) (1 — bss C.) is: (4.2.12) 


where the used symbols denote: 
k= 2 (b,+b3 A») 52+ (1 — bs5 CG) (O52 + bss); (4.2.13 a) 
K,= (63; +533 B,) (63 +-b33 Anas +633 Dn) (1 — bss Ca) (4.2.13b) 
Ky= (bg +33 An) (bs2+ 058) + 2 G52 (1 — bss Cr). (4.2.13¢) 


The remaining procedure is the same as in the previous sections; it refers to an 
infinite system of sets of equations and gives a system of formulas for: 


(€4n es B tan Pan> (n= ie 2, 3, see ), (4.2.14) 


separately for each » equal to 1, 2, 3, ... The parameters A,, B,, Cn, D,, should be 
given. There may be possible to propose some procedure for calculating these para- 
meters on the base of energy considerations, but this will not be discussed here. 
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5. Appendix 


A Procedure of Obtaining a First, Rough Approximation to the 
Coefficient d, [equation (4.2.5) ] 


Using the technique of complex variable functions, Brong? considers the influence 
of the inlet distortion upon the axial compressor surge. The considerations refer to an 
incompressible flow domain. In spite of the fact that the results of Brong were not 
yet generalized to the compressible flow regime, in a very rough, first approximation 
one may apply his results to the compressible flow. Anyhow, since in the final step 
of the calculations one has to apply an iteration 
procedure, the initial step of these calculations 
may begin with the results obtained in the 

Fig. 2. Resultant Relative Inflow. incompressible flow. The symbols used in 

Velocity Pattern Brong’s report are re-adjusted to the symbols 
used in the present paper. 

Brong assumes that the variation of the geometrical angle of attack in front of 
the rotor blade row, and the variation of the resultant relative flow velocity in the 
region of inflow velocity are in the form of a stepwise function (Fig. 2). The variation 
of the effective angle of attack, when the influence of the induced velocity due to the 
shed vortices is taken into account, follows more or less the dotted line. 

The following symbols are used (Fig. 3): 


Piers i Son effective angle of attack in front of a blade row; 

ES ae ae minimum geometrical angle of attack; 

Ps ieee maximum geometrical angle of attack; 

Vien see resultant relative flow velocity in the region of high inflow velocity; 
W ga nice resultant relative flow velocity in the region of low inflow velocity; 
| Seen resultant relative flow velocity behind the rotor. 


Brong derived the following formula: 


sin Bez 25 (V 02/V 12) sin Boz oY Vs tll - 
sin Piz — (Voo/Viz) sin Boo 1+ Vsti/l 


hy, hol Fe , J 


A= Const 


(5.1) 


Knowing the quantities: Bo.; B12; Voo; Vie; 
3 V,; or rather sin Bjo) sin. Bisd Vos? Veep ee 
from the approach discussed in this paper, 
we can calculate sin Bf... We can express 
the magnitudes sin Bo, sin B,., in terms of 
Sere : trigonometric series. The velocities V 93, V 45, 
V, are expressed in terms of trigonometric 
series. The magnitude sin f,2 (or tan £2) 

i will obviously be expressable in the form 
dry, | of trigonometric series. 

Assume a flow with a circumferentially 
uniform inlet velocity distribution. Then 
(see Fig. 3) the angles 6, and £,, are fixed for some fixed velocity pattern in the axial 
compressor. Then also the angle: 


Be2+ (2/2 — B,) =const.= K, (5.2) 


Fig. 3. Nomenclature 


2 Brong, E.: The Effect of Circumferential Inlet Distortion on Axial Compressor Surge. 
General Electric Co., Aircraft Gas Turbine Division, Jet Engine Department, Technical Information 
Series, No. R 57 AGT 645, 16, Sept., 1957, Evendale, Cincinnati, Ohio. 
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is constant. This angle will be constant in the compressor rotor regardless of whether 
the inlet velocity distribution is or is not circumferentially distorted. 

Knowing the magnitude K, and using the formula (5.1), one can calculate the 
magnitude sin £,2 (expressable in terms of trigonometric series) from formula (5.1) and 
the magnitudes of the angle £, and of sin 6. in the case of a circumferential inlet dis- 
tortion from formula (5.2). Obviously, sin £, will be again expressable in terms of 
trigonometric series. 

Next, we can pass formally to the absolute velocity components and apply the 
formula (VI — 9)* p. 36 in Reference? (see also Fig. II — b, p. 3 in Ref.#) to calculate 
the angle £,, corresponding to the angle f, in Fig. 3. Finally, we can return to the 
relative velocity components, thus obtaining the angle f, and finally sin f, (or tan Bs) 
in Fig. 3 in terms of trigonometric series. 

The outlined method allows one to obtain the first approximation to the angle f; 
and sin £; (or tan B;) in the case of the circumferential inlet distortion. 

In the method, outlined above, the approach used is rather stationary or quasi- 
stationary, since practically the interference phenomena due to the nonstationary 
character of the flow are neglected. 

The question how to express the variation of the geometrical angle of attack and 
of the relative flow velocity in the inlet cross-section will be discussed next. 

To preserve a continuity with the proposal discussed in Section 4 of this paper, 
it seems that the application of a series analogous to that one in eq. (4.2.5) is the most 
proper way of doing that. 

The used symbols will denote: 


WESC 8 resultant relative inflow velocity (having Vy, as the maximum and V,, 
as the minimum); 
Pawn cain geometrical angle of attack in front of the blade row. 


Instead of dealing with the angles, we shall deal with the tangents of the angles. 
Thus, we may propose: 


Ve (yy=thyot+ a [hn cos ("") + ], sin (““)| : (5.3) 
tan B. (y) = 2 M+ ) [me cos ("") + pn sin ("") |; (5.4) 
where the coefficients ky, In, Mn, Yn, are given by the formulas: 

| ee = fe. “COs (=2) Vonly) dy vate Oedint das (5.5) 
eee ={" “sin | Via (yy) Gy. == 0, b, 250 ae (5.6) 
i +. [ "cos (*") tan Bs (y)dy; m=0,1,2,...! (5.7) 
=H locecitad MI dy; n=0,1,2 (5.8) 

Pr==)_, SIn\| an 65 (yay, w=O0n1, 2%.-. , 


Using these formulas, we may apply the procedure, described above, and to calculate 
B, and tan B;. But comparing this proposal with the formulas (4.2.3) and (4.2.4), we 
see that what we actually calculate is tan (6;-+ 63), where ; really refers to the mean 


* See eq. (5.46) at the end. 
3 Novak, R.A.: Some notes on the Fundamentals of the Design of Multistage Axial Com- 


pressors. General Electric, AGT, Jet Engine Department, Mimeographed Notes, 1955, Evendale, 
Cincinnati, Ohio. 
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flow; this value is supposed to be known. Hence, this procedure enables us to find 
the value of tan 6,=d, in terms of the trigonometric series with the coefficients k,,, 
ln, Mn» Pn. Comparing the result, so obtained, with eq. (4.2.5) one can obtain the 
magnitudes of each of the coefficients A,, B,, C,, D,, for fixed values of ¢3, and @5n 
[i.e., compare the terms containing sin (n yr-1) and cos (n yr-")]. 

Another way of approaching the solution of this problem is to consider the downwash 
behind each blade. Let us at first consider the results of the elementary theory of the 
downwash as given by the airplane wing airfoil theory. 

The first, very crude approach to the downwash phenomenon is the following one: 

Assume that the action of the airfoil on the air is confined to that part of the air 
which passes through a certain area S perpendicular to the direction of the motion, 
depending in its demensions upon those of the airfoil. In passing through this area S 
the air experiences the reaction of the lift ZL. In consequence of this there is a downward 
momentum communicated to it in a unit of time equal to L. In the steady state 
conditions the mass of the air which passes through S in a unit of time is equal to 9 V S. 
Let the symbol w, denote the downward velocity communicated to it, thus: 


d= 0S Vidi Geils LO ada): (5.9) 


The downward motion implies a kinetic energy. The amount of this energy communi- 
cated to the air in a unit of time is: 


B= st 0.8 ty. (5.10) 


In generating this energy, certain amount of work must be done. Calling this amount 
in unit time D; V, and calling D; the induced resistance, we have: 


D;=E[V =} oS w2=}t Lo V?8. (5.11) 


By some reasoning not given here, we may get the area S in the form: 


S= 71>", shalt span, (5.12) 

hence, 
D;= L?|(2 0 V2.2 b2), (5.13) 
Wo = Lo V x b2). ; (5.14) 


This downward motion cannot be generated at once; it is usually assumed that half 
of it is imparted to the air before it reaches the airfoil, the other half being generated 
after the air has passed the airfoil. Hence, at the airfoil itself there exists a downward 
velocity of the magnitude: . 


w= wi2=—Li(2 0 wb). (5.15) 


When this latter velocity is combined with the original velocity, V, of the air we obtain 
the resultant velocity V (V?)+ w,?/4 sloping downward under an angle y, given by: 


tan p=wy/(2V), or: p=wy/(2 V). (5.16) 
We deduce the relation: 
D;=Lw|(2V)=Lo; p=D,IL. (5.17) 


Thus, the effective angle of incidence is: 


%y = A—. (5.18) 


Compressibility Effects in Circumferential Inlet Distortion in Axial Compressors 233 


Let us express in the formula: 


D;= L?|(2 9 V? x 62), (5.19) 
both Z and D; by means of the corresponding coefficients, then: 


C= oa S . AR (2 b)2 MS 
a bp > ot aes <a (5.20) 
C? 
Cp be 2] 
a) eA ie Ce 
Thus: : 
D; Cp: CL ee: 
ea I CL aAR?’ (5.22) 
IL 


ae: (5.23) 


Concerning the application of the formulas given above, to the axial compressor 
blades, the following changes may be introduced: 


Since one considers the conditions at the trailing edge of a blade, one may assume 
in place of eq. (5.15) the formula: 


Wi, = Wy = L|(o Va 2), (5.24) 


with the following derived formulas: 


tan Yr Py = W, V-; (5.25) 

Dew Ne) \=4 Los: oy =—2 Die (5.26) 
Cor =4— Pr; (5.27) 

Pr =2 DJL = 2 Cyi|Cr=2 Ci/(a AR); (5.28) 
Lo = %—Py =a4—2C_/(x AR). (5.29) 


Moreover, the velocity V should be taken to be the local velocity which varies along 
the blade. Also, the coefficients C; and Cp; should be taken as the local aerodynamic 
coefficients which may vary along the 
blade. In case the incoming velocity is a 
function of time ft, this dependance should 
be included, although the final result, 
eq. (5.29), does not depend upon V = V (i). 


Similarly, if the angle of attack « is a func- Z 
tion of time, i. e., «=« (t), this fact should N Bite ie 
be taken into account. The above reasoning : . / 
does not include the interference pheno- S 

mena, cascade effects, etc., which may Fig. 4. Induced Velocity 


partly change the final conclusions. 
There exist more refined formulas for the calculation of the downwash behind 


a wing (see Durand?, Vol. 2, p. 151). Let us quote below a result from reference* 
(Durand, Vol. 3, p. 80). Using the Biot-Savart Law the following result may be 
derived for the induced velocity (see Fig. 4): 


ie ae 1 eee 
dw,= Fy 4a (008 ¢ + Cos DUT sere aa) (1+ sin g) + 
1 x 
d 5.30 
By <a ee ee ere 


4 Durand, W.F.: Aerodynamic Theory, Vol. 1 to 6, Berlin: Springer. 1935. 
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where: 
: 1 : l 
sin 9 = —————— 5 sn YP = 5 (5.31) 
ig Vai + yz)? +P 4 Rn = og)? 4a 
(Yr + Yo) : bee (Y1 — Yo) ¢ 
cos 9 = ——34 es cos Y= (5.32) 
4 Vir + yo)? +P id V1 — Yo)? +P 


When considering the conditions at the trailing edge, the value of / may be assumed 
to be equal to approximately 50 percent of the cord length. In the case of an axial 
compressor blade one may assume that y,=lenght of the blade and that there is only 
one tip vortex attached to the blade, the other end of the blade being attached to 
an infinitely large wall (the hub); this gives: 


we AR Wf dl 1 : 
E cos e+ 4a, i+ snp] dy. (5.33) 


The value of the velocity w, at the point P is given by the formula: 


1 ] 


BS Al : 
ae, Oy, 4x lz COs = (ith) (1 + sin y)| dy. (5.34) 


The circulation J’ is related to the lift coefficient of an infinitely long span wing by 
means of the well-known formula: 


Pat Or V ¢ (5.35) 


where c is the length of the chord. In the case of a finite span wing one can use the 
formula given in Ref.>, p. 69: : 


Pine (Y)= 4 U Cine tine (tine — “F*), (5.36) 
where the subscript “ine”? denotes the incompressible flow conditions, and the profile 
constant, in, is the slope of the lift coefficient curve for the profile; to the first 
approximation it depends only on the Reynolds number of the flow and the character 
of the boundary layer on the wing. One usually assumes xn. to be known. 

Concerning the influence of the compressibility (p.70 in Ref.*) the following 
formulas may be recommended: 


Thi a (5.37) 


where /’(y) denotes the circulation in the compressible fluid flow system ; 


cin (y= ely) Bs pt=(1— 2), (5.38) 

dine (Y) = & (Y); ; (5.39) 
: Line (y) = 

Wine (y) = 7 | (y aS dy, (5.40) 


where the quantities without the subscript refer to the compressible fluid flow and 
the symbol ‘8 indicates that the Cauchy principal value is to be taken; 


dey ai I 6 IY) 4- 
P(W= 95 U 6 rine [0 mr ® |G dy|. (5.41) 


Sears, W. R., Editor: General Theory of High Speed Aerodynamics, Vol. VI of High Speed 
Aerodynamics and Jet Propulsion, Princeton University Press, 1954. 
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For a wing of elliptic circulation distribution one gets the formulas: 


ae Chine = 
Wing = COnst. =U, ay Wea (5.42) 

Ola etes “ine : . ~ 
Line ig Zine (at A Reine) Ones (5.43) 
w=U O|(a AR): (5.44) 
Cp Bt wml AR) & (5.45) 


Eq. (VI-9) on p. 36 in Ref.3 has the form: 
By = tan [2 (1+) tan By, + (1—g) (1+) tan B,], (5.46) 


where: 
m = Stagger angle; 
6, = entering angle (between the axial direction and the relative entering flow 
direction); 
B, = leaving angle (between the axial direction and the relative leaving flow 
direction); 
q=K [2 ¢(l2.cos B,)=1]-; 
1 = length of the chord of a blade; 
= spacing distance; 
K =f (t/l, Bm) = lattice effect coefficient, known from the design characteristic 
properties ; 
C= 2 K-Si0_0; 
6 =the angle between the f,, and the f, directions. 


(Received, July 10, 1959) 


Die Stabilitat der axial gedriickten Kreiszylinderschale mit 
verdinderlicher Wandstarke 


Von H. Wagner, Karlsrnhe 
Mit 5 Teatabbildungen 


Zusammenfassung. Es wird die elastische axialsymmetrische Ausbeulung des Kreiszylinders 
mit verinderlicher Wandstirke behandelt. Dabei wird auf die Reduktionsmethode (Matrizenkalkul) 
zurickgegriffen. Durch Einfithren des Integralmittelwerts fir die Wandstirke lassen sich einerseits 
bekannte Beulwerte fiir die ‘einfachsten Wandstarkegesetze nachweisen, andererseits wird damit 
auch die Berechnung von Beulwerten bei komplizierterem Wandstarkenverlauf erméglicht. 


I. Einleitende Bemerkungen 


Wenn auf eine kreiszylindrische Schale (Rohr) an den beiden Randern gleichmaBig 
verteilte Druckkrafte in Richtung der Erzeugenden wirken, so wird man mit einer 
zweifachen ,,Knickgefahr‘‘ rechnen miissen. Einmal kann dieser Zylinder, sofern er 
lang genug ist im Verhaltnis zu seinen Querschnittsabmessungen und sofern die auf- 
gebrachten Druckkrafte einen gewissen , kritischen‘‘ Wert iibersteigen, wie ein Stab 
in einem Bogen ausknicken (Euler-Knickung). Dabei hat man vorauszusetzen, dab 
die zugehérige ,,Knickspannung“ in den elastischen Bereich fallt. Zum anderen aber 
ist es auch méglich, insbesondere wenn Wandstarke, Radius und Rohrlange zusammen 
mit den Druckkraften in bestimmter Gré8e auftreten, da sich bei gradlinig bleibender 
Zylinderachse in der Zylinderwandung ringsum Falten bilden, welche bei einer fort- 
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dauernden Belastung zum Bruch fiihren (Abb. 1). Diese Art von Knickung, welche man 
dem allgemeineren Fall der nichtaxialsymmetrischen Ausknickung (wobei die Quer- 
schnitte nicht Kreise bleiben, sondern sich zu Vielecken verformen) subsumieren kann, 
nennt man anschaulich ,,Beulung‘‘. Fiir kreiszylindrische Rohre mit konstanter Wand- 
starke unternahmen fiir diese axialsymmetrische Beulungsart Lorenz! und Timo- 
schenko? erstmals theoretische Untersuchungen. 

Nur von der zuletzt genannten (elastischen) axialsymmetrischen Ausbeulung des. 
Kreiszylinders handelt die vorgelegte Arbeit”. 

Die mathematische Behandlung der Eigenwertprobleme der Elastizitatstheorie 
geht davon aus, daB man die Gleichgewichtsbedingungen fiir den verformten 
Zustand des betrachteten Tragwerks aufstellt. Danach ergeben sich die zugehérigen 
kritischen Werte als Eigenwerte eines Randwertproblems. An Stelle der Gleichgewichts- 
methode zur Aufstellung der das Problem beherrschenden Differentialgleichung kann 
man sich auch der Energiemethode bedienen, und dies mit zweifachem Vorteil: 

Erstens hat diese Methode den Vorzug, da mit einem verhaltnismafig einfachen 
Ansatz die Differentialgleichung durch formales Rechnen gewonnen wird. Die gesuchte 
Differentialgleichung ist namlich die Eulersche Gleichung, welche dem formulierten 
Variationsproblem vom Minimum der potentiellen Energie zugeordnet ist. Dieser 
Lésungsweg wird oft als der indirekte bezeichnet?. 

Zweitens aber, und dies ist die manchmal wichtigere Ausniitzung der Energie- 
methode, gestattet diese eine direkte Behandlung des Variationsproblems, indem sie 
zur Gewinnung von Naherungslésungen herangezogen werden kann. Fiir diesen Zweck 
leistet auch die Anschauung vorziigliche Dienste (Ritz-Ansatze), wenngleich auch 
manchmal der (elementare) Rechenaufwand betrachtlich ansteigt. 

Fir den Fall der Kreiszylinderschale, wie sie im folgenden betrachtet werden soll, 
resultiert aus den vorstehenden Betrachtungen eine gewohnliche lineare homogene 
Differentialgleichung vierter Ordnung als Randwertproblem. Diese Differential- 
gleichung ist bekanntlich dieselbe wie jene fiir den Balken auf nachgiebiger Unterlage’. 
Jedoch nur bei konstanter Wandstirke ist das Eigenwertproblem einer leichten 
Behandlung zuganglich, weil dann die Differentialgleichung konstante Koeffizienten 
besitzt®. 

Wesentlich erschwert wird die Aufgabe bei Vorgabe eines Wandstirkegesetzes. 
So haben erstmals Biezeno und Koch‘ die Stabilitat eines zylindrischen Fliissigkeits- 
behalters mit bereichsweise konstanter Wandstarke untersucht. Der Einflu8 einer 
stetigen Veranderlichkeit der Wandstirke bei einem Kreiszylinder unter axialer 
Druckbelastung wurde in neuerer Zeit von Federhofer’ dadurch in Betracht gezogen, 
daf er nur geringe Querschnittsanderungen zulieB, wobei er dann mit Hilfe der Stérungs- 
theorie fiir das Eigenwertproblem geschlossene Naherungsformeln fiir die Beullasten 
bei lmear und quadratisch zunehmender Wandstarke entwickeln konnte. Dabei — 
um es zu wiederholen — beschrinken sich alle Aussagen allein auf elastische, nur 
axialsymmetrische Ausbeulungen des Kreiszylinders. 

In der vorgelegten Arbeit wird nun die Forderung einer nur geringen Wandstirke- 
anderung nicht mehr so strikt gehandhabt werden miissen. Es wird sich zeigen, da 
erstens auf einfache Weise im wesentlichen dieselben Naherungen fiir die Beullasten 
erzielt werden wie bei Federhofer, wenn man sich auf geringe Anderungen der 
Wandstarke beschrankt, und dafi zweitens die Behandlung des allgemeineren Falles 
der beliebigen** Wandstarkeinderung méglich ist. Dies geschieht folgendermaBen: 

Die Kreiszylinderschale wird durch senkrecht zur Zylinderachse liegende Schnitte 


* Man kann diesen Fall den ,,klassischen‘‘ nennen. 
** | beliebig*‘ in dem Sinne, da® der Charakter als Schale gewahrt bleibt. 
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in solche Teilzylinder zerlegt, in denen sich die Wandstarke nur gering andert. Dabei 
kann das Wandstirkegesetz eine beliebige* integrierbare Funktion der Langenvariablen 
der Schale sein. Fiir einen solchen Teilzylinder wird die dem Problem zugeordnete 
Differentialgleichung gelost, indem fiir die Koeffizienten der Differentialgleichung die 
linearen Mittelwerte fiir diesen Integrationsbereich herangezogen werden. Damit 
erhalt man in diesem Bereich z. B. fiir ein lineares Wandstarkegesetz denselben Kigen- 
es bei quadratischem Wandstarkegesetz fast denselben Kigenwert wie bei Feder- 

ofer. 

Das Aneinanderfiigen dieser Teilzylinder zum gesamten Kreiszylinder geschieht 
dann mit Hilfe des Reduktionsverfahrens, dessen adaquates Ausdrucksmittel die 
Matrizenrechnung ist. Mit diesem Reduktionsverfahren ist in neuerer Zeit eine Reihe 
wichtiger technischer Probleme gelést worden (vgl. 8 und die dort angegebene Literatur, 
ferner® und den Ubersichtsartikel ). 

Das Reduktionsverfahren 148t sich auf solche Probleme anwenden, die durch 
gewohnliche lineare Differentialgleichungen beschrieben werden kénnen. Kann man 
ein Fundamentalsystem der dem Problem zugeordneten linearen Differentialgleichung 
angeben, wie z. B. bei Differentialgleichungen mit konstanten Koeffizienten, so ist 
eine Behandlung nach diesem Verfahren zweckmabig. 

Will man danach komplizierte technische Probleme mathematisch einfach formu- 
lieren, so hat man dem gegebenen System ein Ersatzsystem zuzuordnen. Ein solches 
Ersatzsystem liegt z. B. vor, wenn man eine stiickweise konstante Verteilung der 
einzelnen GroBen (wie etwa der Steifigkeitsbelegung) vornehmen kann. Die genannten 
Teilzylinder stellen in diesem Sinne das Ersatzsystem fiir den gesamten Kreiszylinder 
dar, wobei der Ubergang der zu verfolgenden Gro%en (Krifte, Auslenkungen usw.) 
von einem Teil in den nachsten durch Matrizenmultiplikation geleistet wird. Damit 
hat man es in der Hand, gegebenen Forderungen nach ausreichender Genauigkeit 
der Rechnung durch entsprechend feine Gliederung des Ersatzsystems zu begegnen. 

Dies alles bedeutet fiir die gestellte Aufgabe ein einfaches Transformationsproblem 
gewisser einspaltiger Matrizen (Spaltenvektoren), wobei die Transformation selber 
mittels sogenannter ,, Ubertragungsmatrizen“ und ,, Ubergangsmatrizen‘“‘ vorgenommen 
wird. Die Beulbedingung fiir den gesamten Kreiszylinder stellt sich dann dar durch 
Nullsetzen einer zweireihigen Unterdeterminante einer Multiplikationsmatrix, der 
sogenannten ,,Beulmatrix‘*. Zwar wird diese Bedingung, insbesondere bei mehr als 
zwei Teilzylindern, sich in wenig iibersichtlicher Gestalt prasentieren. Dieser Umstand 
begegnet jedoch hinsichtlich der numerischen Prozedur dann keinen sonderlichen 
Schwierigkeiten, wenn man den Einsatz von digitalen Grofrechenanlagen ermdég- 
lichen kann. 


II. Die Differentialgleichung fiir axialsymmetrische Deformation 


Ks sei 

ie die auf die Liangeneinheit des Schalenrands bezogene, tiber den Rand gleich- 
maBig wirkende axiale Druckkraft, 

Py, die kritische Last (,,Beullast*‘), 

l die Lange der Kreiszylinderschale, 

(p der Radius der Schalenmittelflache, 

a die Langenvariable (ec =0 und x=/ kennzeichnen die Rander der Schale), 

h (x) die veranderliche Wandstirke der Schale, 

w (x) die Radialverschiebung der Schale, positiv gemessen nach auBen von der 


unbelasteten Schalenmittelflache aus, 


* beliebig’’ in dem Sinne, daB der Charakter als Schale gewahrt bleibt. 
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bo 
i) 
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die Umfangsrichtung oder tangentiale Richtung fiir einen Schnitt «= const., 
die Flachennormale der Schalenmittelflache, 

der Elastizitatsmodul, 

die Poissonsche Zahl, 

die Formanderungsarbeit der inneren Krafte, 

die Arbeit der auBeren Krafte, 

das axiale Flachentragheitsmoment fiir einen Stabquerschnitt, 

die Spannung, 

die Dehnung, 

auf die Langeneinheit einer Zylindererzeugenden bezogene Kraft in 
Richtung ¢, 

aT die Komponente von 7 in Richtung 1, 

has ein Mittelwert der Wandstarke. 


Hea Qk ke st 


Abb. 1 Abb. 2 


Fir die Herleitung der Differentialgleichung ist es zweckmafig, sich aus der 
zylindrischen Wandung einen Liangsstreifen der Breite 1 cm ausgeschnitten zu denken 
(Abb. 2) und fiir diesen ,,Stab“ die Differentialgleichung aufzustellen. 

Diese Herleitung soll nun mit Hilfe der Energiemethode geschehen, damit die hierfiir 
zu formulierende Variationsaufgabe spaiter (in VII) fiir die dortigen Zwecke erweitert 
bzw. fiir Naherungslésungen verwendet werden kann. 


Das zu einem Minimum zu machende Gesamtpotential /7 =A; — A, setzt sich aus 
folgenden Anteilen zusammen: 


1. Forminderungsarbeit A; der inneren Krafte. Hierzu zahlt ein Term, welcher 
von der Biegung herriihrt, 


l 


ae 
A) =| Jw! de, (1) 


Q 
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worin wegen der Beschriinkung der Dehnungsméglichkeit des ,,Stabes‘‘ an Stelle Z 
in bekannter Weise zu setzen ist 


E 
iT a (2) 
Der angenommene Querschnitt des Stabes mit rechten Ecken muB8 namlich ein solcher 
bleiben, d.h. der Spannungszustand kann nicht einachsig bleiben, wenn der Stab 
sich mit den anschlieBenden Zylinderteilen in Zusammenhang befinden soll. 

Infolge der tangentialen Dehnungen , 
werden Umfangsspannungen oa, hervorgerufen. 
Deren iiber die Wandstirke genommene Resul- 
tierende T' ist auf die Langeneinheit bezogen: 


Es 


T =oh=e Eha=*w, (3) 


Diese Kraifte T hat man sich an Stelle der 
fehlenden anschlieBenden Zylinderteile langs 
des Stabes angebracht zu denken. Sie haben 


eine Komponente d7' in der Flachennormale n 
vom Betrag (Abb. 3) 


Abb. 3 


Eh 
Oe Me Ww. (4) 


Diese Komponente (4) stellt einen der Ausbiegung w proportionalen Widerstand 
dar, der das Ausbeulen zu hindern versucht. In (4) spielt der Faktor = dieselbe Rolle 


wie die Bettungsziffer in der Theorie des Balkens auf nachgiebiger Unterlage. 
Damit laBt sich dann der zweite Term der Formanderungsarbeit angeben zu 


l 
A® = Ta hwdz. (5) 
0 


2. Arbeit A, der A4uBeren Kriafte. Lediglich die Arbeit der auBeren Krafte P, die 
auf dem Verschiebungsweg geleistet wird, gibt hierfiir ihren Beitrag 
l 
age | w’2?da*. (6) 
0 


Damit ist nun die Forderung (das Variationsproblem) gegeben 


l 
=| eee w''? (v) + BN) (x) — Pw’? (x) dx = Min. (7) 


2 


T= A;—A,= 


0 


Die Eulersche Gleichung fiir das Variationsproblem 
I 
|F (x10, wi, w’) de = Extr. 
0 


lautet 


a2 (37 d (3) ! oF a (8) 


az \aw’’ du \aw’ aw 


* Der Faktor 4% kommt aus der Taylor-Entwicklung fir den Verschiebungsweg (ds — dx) 
mit ds = 1+ w® dz, d.h. die Krafte P sind bei der betrachteten Verschiebung von Anfang an 
voll wirksam. 
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Fiir (7) nimmt sie die Form an 


a2 EB h3 (x) 
da? |12(1—» 


aE (2) = 0. (9) 


2) w+ Pw (x)| + 


Dies ist die bekannte Differentialgleichung fiir die axialsymmetrische Deformation 
des Kreiszylinders!*. Damit die Differentialgleichung (9) dem Reduktionsverfahren 
zweckmakig angepaBt werden kann [das will besagen, damit man fiir (9) ein Funda- 
mentalsystem bequem angeben kann], wird fiir h (w) der lineare Mittelwert hj, tiber 
den in Frage stehenden Integrationsbereich gewahlt 


i) 
hu=> [ih (a) dae. (10) 


Uber das Wandstiarkegesetz h (x) wird zu Beginn des nachsten Abschnittes das 
Notige gesagt. 
Beachtet man (10) und schreibt man zur Abkiirzung noch 


Par ae ad 2b as P, (11a) 
po= rug (11b) 
so erhalt man aus (9) die Differentialgleichung 
wll (x) + 2 ap-w"’ (x) + Baw (x) =0. (12) 
Die Randbedingungen fiir das Problem seien durch 
w(O)=w' Hw Vw’ (y=0 (13) 


gegeben, es handelt sich also um ein beidseitig gelenkig gelagertes Rohr. Der gesuchte 
niedrigste EKigenwert ist durch P= P;, gekennzeichnet; er ist nach (11a) proportional 
ZU Oo. 


III. Die Lésung der Differentialgleichung fiir einen Teilzylinder 


Zuvor wird folgende Einschrankung vereinbart: Die Kreiszylinderschale mit 
veranderlicher Wandstarke (Abb. 4) werde so in m Teilzylinder aufgeteilt, daB innerhalb 
eines solchen 7-ten Teilzylinders (7 =1, 2,..., m) sich die Wandstarke nur gering 
andere, wobei die Aussage ,,geringe Anderung“ noch priazisiert werden wird. 
Betrachtet man nun die i-te Teilschale mit der Linge /; fiir sich, so sei fiir diese Teil- 
schale das folgende Wandstirkegesetz gegeben 


h; (x) =hy-1 [1+ wu fi (x)). (14) 
Es bedeuten f; (2) eine monotone stetige Funktion mit 
fi (wi-1) = 0 
und 
fi (wi) =1, 


ferner h;_, bzw. h; die Wandstarken am Anfang bzw. Ende der i-ten Teilschale sowie 
ui einen Parameter, welcher die relative Anderung der Wandstirke im Intervall 
[w;-1, %] angibt, 


hi 
a= hina Sir (15) 


» 
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Die Abmachung, daf sich innerhalb des i-ten Teilzylinders die Wandstirke nur 
>gering™ andere, soll dann besagen: 


Es sei y; eine kleine Zahl derart, da8 


per<i (16) 


gelte. Dies bedeutet, daf in allen folgenden Rechnungen y;? gegeniiber uw; vernach- 
lassigt werden darf. 

Im Gegensatz zu’ muB also die Forderung (16) nicht fiir die gesamte Kreiszylinder- 
schale gelten, sondern nur innerhalb einer Teilschale. 

Im folgenden werde nun die Integration der Differentialgleichung (12) mit den 


Randbedingungen (13) fiir eine einzelne Teilzylinderschale durchgefiihrt. Wegen der 
Ubersichtlichkeit werde dabei vereinbart, 


dafi zunachst tiberall der Index 7 weg- 
gelassen werden soll. Nach dem Integra- 
tionsproze8 wird er dann wieder — wo 
notig — eingefiihrt. Dies bedeutet dasselbe, 
als ob die Integration iiber die gesamte 
Kreiszylinderschale erstreckt wiirde, d. h. 
die Variable x aus [0,/] genommen wiirde*. 

In der Differentialgleichung (9) war 
fiir die veranderliche Wandstarke h (2) 
der lineare Mittelwert (10) genommen 
worden, woraus Gleichung (12) entstand. 
Bedeutet nunmehr 


Stelle i 
Stelle i‘ 
j-ter Teilzylinder 

V Stelle i-1 4 


Stelle(i-1)' 


SSS: 


é= 


ae 
1 ? 


so hat man wegen (14) fiir den linearen 
Mittelwert (10) 


| 


hu=hof +n f (Olde | (17) 


0 
hu=h, [(l\+uF], Abb. 4 


wenn man setzt 
1 


F=\f (64. (18) 
0 
Hierin ist h, die Wandstarke am Anfang der Schale (Abb. 4). (Spater hat man dafiir h;—1 
zu setzen, wie aus den vorstehenden Bemerkungen ersichtlich ist.) 
Die Differentialgleichung (12) schreibt sich somit, wenn Striche jetzt Ableitungen 
nach & bezeichnen, 


wil¥ (£) +2 aw" (€) +f? w (g)—0. (19) 
Die Randbedingungen sind 
w (0) =w" (0) =w (1) =w"” (1)=0 (19a) 
ferner hat man jetzt 
d 12(1—) 2 5 
S aan Ty eee) 


* Nach dieser Bemerkung wird wbrigens auch verfahren. 


¢ 18 
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12 (1 —»2)T 


9 
re Nas 


pi (20b) 
mit ha nach Gleichung (17). 
Die Integration der Differentialgleichung (19) tiber die charakteristische Gleichung 


fiihrt zur allgemeinen Losung 


4 
w(é)= 24% (6), (21) 
fa 
wenn die Funktionen w, (€) die Fundamentallésungen sind. In (21) stehen die vier 
mathematischen Integrationskonstanten c,, die aber im Sinne der Mechanik nicht 
direkt deutbar sind. Wollte man sie durch 


We == (0); > wh ey ne te 


ausdriicken (diese GréBen haben eine direkte mechanische Bedeutung, sie sind namlich 

proportional den Auslenkungen, Kraften usw.), so wire einige Rechenarbeit erforderlich. 

Es wird also zweckmabig sein, die Differentialgleichung (19) als Anfangswertproblem 

zu behandeln. Der Kalkiil hierfiir ist in der Laplace-Transformation gegeben; diese 

gibt die Lésung der Differentialgleichung (19) direkt in den AnfangsgroBen wy, w%, 

wi, und w,”’ an. Diese Tatsache ist fiir die weitere Rechnung geradezu erwiinscht**. 
Wahlt man noch die Abkiirzungen 


a=a—| «— p? (22a) 
b=at) a2 — B2 (22b) 


und setzt man tiberdies «? = 8? > 0 voraus*, woraus a? > 0 und b? > 0 folgen, so lautet 
die Lésung von (19) mittels Laplace-Transformation 


b2cosaé — a®cosbé , Bsinaé—assinbé ,, cos a& — cosbé 
w (€)= Wy hi ga Wo a b (62 — a?) + Wo b? — a? 7a 


1, bsina&é—asinbé , 
+ Wo a b (b? — a?) (23a) 


und die Ableitungen bestimmen sich danach zu 


bsinb&é—asinaé 


rg a2bsinb&é —ab?sina& , 2 cosa&—arcosbé i 
Ww ( ) = Wp b2 — @2 Wo b2 — q2 Wo b2 — a2 + 
17, C08@E— cosbé 
Slaw eer pe ae (23b) 
w"' (é) a@b2cosbE—arb2ecosaé ,a@bsinbé—ab?sinaé ,, bcos bE— ar cosaég 
lamers 0. b2 — a2 Wo b2 — a2 Wo b2 — a2 = 
bsinbé— asinaé 
| Lage 2 . 
wo — (230) 
wi(E) = w, a b2 sina & — a2 be sin b & ern a2 b2 cosbE — a®b? cosaé l 
b2 — a2 ial b2 — a2 
Lop , @sinaé— bsinb& pe b2cosb&— a? cosaé 
0 b2 — a2 0 b2 — a2 


(23d) 


ei Der Fall a?<#? hat fir das vorliegende Eigenwertproblem keine Lisung, wie man durch 
austihrliche Rechnung nachweisen kann. Daf tiberdies a?>0 ist, folgt (ganz abgesehen von der 


technischen Bedeutung von a) aus der Selbstadjungiertheit und Volldefinitheit des Eigenwert- 
problems, welche Realitat der Eigenwerte verbiirgen. 
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IV. Die Matrizenmethode 


Ks erscheint angebracht, zuvor den wesentlichen Gedanken der Reduktionsmethode 
zu erlautern. 

Durch den in III vorgefiihrten IntegrationsprozeB lassen sich die Kraft- und 
Deformationsgrofen, deren Anderungen man gerade verfolgen will, am Anfang des 
_t-ten Teilzylinders (bei x;_1) durch jene am Ende (bei 2;) ausdriicken, und zwar ist 
diese Beziehung linear. Dieser Umstand legt auch die Anwendung der Matrizen- 
rechnung nahe. Entsprechend lassen sich nun die Kraft- und DeformationsgréBen 
am Ende 2; des -ten Teilzylinders auffassen als die Anfangsgréfen fiir den benachbarten 
(1 + 1)-ten Teilzylinder, jene wiederum fiihren durch die Integration einer Differential- 
gleichung fiir den (¢+1)-ten Bereich zu den EndgréBen dieses Teilzylinders usf. 
Endlich kann man also die Kraft- und Deformationsgré%en am Ende des gesamten 
Kreiszylinders durch die des Anfangs ausdriicken in Form linearer Beziehungen. Die 
Beulbedingung selber ergibt sich dann einfach durch Nuilsetzen einer zweireihigen 
Unterdeterminante einer Produktmatrix. 

Kehrt man wieder zum i-ten Teilzylinder zuriick, so lassen sich in den Glei- 
chungen (23a bis d) fir §=1 (d.h. fiir e=2;, denn im 7-ten Bereich gilt eigentlich 


oe <=)" die GroBen 


Dette (oc 0. 1, 2, o8p eam nt, B, te. 1) 


fiir die Endstelle 1’ des Teilzylinders (Abb. 4) ausdriicken durch die GréBen w'® (0) = 
=w'?, der Anfangsstelle (¢— 1) in folgender Weise 


a i 7 , G vr (i) yer 
w, =A w,_, + Alpw_, +43 wl, + Aw, (24a) 
Wi = ACen Ary erly, Wea Aa (24b) 
Lj i) 7 , :) wr (2) pyy227 ¢ 
wy = Af) w,_, + AS w;_,+ A} wil + AST i”, (24¢) 
wi = Aw, ,+AQw! +49 wi, + AQwl%. (24d) 
Nun definiert man als Zustandsvektor tv; die einspaltige Matrix 
‘Ww. : 
ae 
Wit 
y;-1= j (25) 
w.’ 


(entsprechend !v;) sowie die , Ubertragungsmatrix®‘ {;, welche die Ubertragung des 
Zustandsvektors 1v;-1 in den Vektor vy; vornimmt, 


fae Age A Ap 
x | ah 48 48 48 oe 
wlth ocAlttee dean 

AQ? ae Ag ag 


* Dies impliziert auRerdem, daB in den Ausdriicken (20a, b) baw. (22a, b) an Stelle J die 


GréBe li zu stehen kommt. 


18* 
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Die Elemente dieser Ubertragungsmatrix 2; lauten nach (23) und (24) 


(i) __ D2 cos aj — a? cos B; | A® bi? sin aj — a,* sin bj 
BO bi? — a;* ; EO Oe 
Al) — 008d = 008 bi UNG pete 
13 b? — az? > 14 aj bi (bi? — ai?) 
, aj? bj sin bj — a; b;? sin a; ; j 
AO ani As} =A}} 
) bj sin b; — a; sin a; (i) (i) 
AY == > An =Aly 
23 bi — a;? : | (27) 
; a,” b;? cos bi — ai? bi? cos. ai (i) — 4) | 
AB heat ae 
. bcos b, — aj? cos ai (i) (i) 
AN = b22 — aj? 2 AN a 25 
a; b;2 sin a; — a? 6? sin b; (i) (i) 
Ay} a bi? — aj? ; Ap=4y 
ee az? S1N A, b;? sin b; f (AG) 
Ble — b?2 — a?” : sie ea 2 


Die Matrix Y; ist dem 7-ten Teilzylinder zugeordnet, der hochgestellte, geklammerte 
Index 7 in (24) bzw. (27) zeigt diese Zugehérig- 
keit an. 

Ubrigens mu8 sein 


det 2; = const., 


weil man die Wronskische Determinante zu (21) 
genau in det YY; umformen kann. Diese Konstante 
bestimmt sich sofort zu 1, wenn man in der Deter- 
minante fiir die Gleichungen (23) =0 setzt. 

Die Ubertragung des anfanglichen Zustands 
(1v;-1) in den des Endes (iv;-) 1aBt sich nunmehr 
mittels (25) und (26) als eine lineare homogene 
Transformation auffassen 


vy, = WW; m;-1*. (28) 


Beim Ubergang vom i-ten Teilzylinder in den 
benachbarten (7 -++ 1)-ten hat man zu beachten, daB 
an dieser AnschluBstelle ein stetiger Ubergang der 

Abb. 5 Kraft- und DeformationsgréBen erfolgt. In der Tat 

bedeutet das hier gewahlte Ersatzsystem (bestehend 

aus m Teilzylindern) nichts anderes als eine Aneinanderfiigung von m Teilzylindern 

verschiedener, aber konstanter Wandstirken, weil hy nach Gleichung (17) fiir jeden 
Teilzylinder einen anderen konstanten Wert h;,, hat. 


higM = 


lis] 


Stelle i 
Stelle i' 


Stelle (i-1) 
Stelle (i-1)' 


* Fur den Fall statischer Probleme, welchen Gleichgewichtsbetrachtungen am unverformten 
Element zugrunde liegen, hat Schumpich® (S. 30) ein allgemeines Verfahren zur Aufstellung 
solcher Ubertragungsmatrizen % angegeben, bei dem nicht von einer homogenen linearen Differen- 
tialgleichung n-ter Ordnung mit konstanten Koeffizienten, sondern von einem dieser Differential- 
gleichung entsprechenden System von m homogenen linearen Differentialgleichungen 1. Ordnung 
mit konstanten Koeffizienten ausgegangen wird. Man kann sich dabei schon von vornherein 
zur Lésung dieses Systems des Matrizenkalkils bedienen, mu jedoch dafir noch eine Inter- 
polationsaufgabe lésen. Der hier beschrittene Weg itber die Laplace-Transformation ist aber 
einfacher, weil er véllig schematisch verlauft. 
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_An der AnschluBstelle ¢— 7’ (Abb. 5) 


mussen demnach folgende Stetigkeitsfor- 
derungen erfiillt sein nach folgende Stetigkeitsfor 


Wy = WU; 
l, Wi =I, W; 
hing Ui wy =h; ou G1 Oy ae 
va Wis = hi, 1,M ran 1 W;". 


Diese Stetigkeitsforderungen (29) sprechen den Sachverhalt aus, da® an den beiden 
, schnittufern®‘ der AnschluBstelle 7— 7’ Verschiebung, Tangente, Moment und Quer- 
kraft gleich sein miissen. Die GréBen w’, w’’ und w’” selber sind namlich unstetig an 
den AnschluBstellen. 

Der Zustandsvektor 1; am Anfang des benachbarten (i+ 1)-ten Teilzylinders 
entsteht dann aus ty, durch Multiplikation mit einer sogenannten Ubergangsmatrix”, 


d. h. die Stetigkeitsforderung nimmt fiir die Zustandsvektoren die Form einer Matrizen- 
gleichung an 


ww; = U; 1p;,, (30) 
wobei die Ubergangsmatrix Ul; lautet 
1 0) O 0) 
Li 
0 = 0 
UW= 0 0 eats 1,M 0 : (31) 
i Mia a Wire 1,M 
0 0 0 AS ike : a teal Sr) 
I? bea 


Mit (30) schreibt sich dann die Ubertragung des anfanglichen Zustands (tv; —1) 
der i-ten Teilschale in den Anfangszustand (1v,;) der benachbarten (7+ 1)-ten Teilschale 
unter Verwendung von (28) wie folgt 


= UU; wp, VW Ww; 1. (32) 


In der zu Beginn dieses Abschnittes geschilderten Weise laBt sich fiir jeden Teil- 
zylinder eine lineare homogene Beziehung analog (32) angeben. War der gesamte 
Kreiszylinder in m solcher Teilzylinder zerlegt worden, so ergibt sich rekursiv unter 
Beachtung der Stetigkeitsforderungen an den inneren AnschluSstellen der Teil- 
zylinder* 

Dre = Wn ies 1 Wn —1ls--- UG A, te WM Wo (33) 
oder kurz 
Wm == Vato. (34) 

Diese Produktmatrix S,, stellt nun den linearen homogenen Zusammenhang her 
zwischen den Kraft- und DeformationsgréBen Ww, und i, der beiden Schalenrander. 
Durch Vorgabe der geometrischen und dynamischen Randbedingungen (13), némlich 
hier 

W = Wi, = wv, = w=, (35) 


* Die Ubergangsmatrizen WU waren zu vermeiden, wenn man den Zustandsvektor (25) so 
definieren wiirde, da® er als Elemente die Kraft- und DeformationsgroBen direkt enthalt. Dann 
waren aber die Dimensionen aller auftretenden Matrizenelemente sowohl in tv als auch in Y ver- 
schieden, wihrend so die Elemente der Matrizen 2% und U alle die Dimension 1 haben, wodurch 
als Vorteil bei der spiteren numerischen Rechnung alle Elemente in etwa gleicher GroRenordnung 
erscheinen. Ubrigens haben jetzt die Elemente der Matrizen tv stets die Dimension einer Lange. 
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erhalt man in dem homogenen Gleichungssystem (34) nur dann eine nichttriviale 
Lésung (welche dem instabilen Gleichgewicht entspricht), wenn die zweireihige Unter- 


determinante 
) 
Kp Ke 


Kis ke 


a (36) 


ist. Dies kann man sich an Hand der Gleichung (34) veranschaulichen, wenn man_be- 
achtet, daB die sogenannte Beulmatrix 8,, wie folgt geschrieben wird 
( 
Tp Kgs Hp 
[ae EP RP we waa 
No.2 3 = ° 
Kp ke Re ie 
EP BD Re ee 


V. Behandlung des Sonderfalls m=1 


1. Allgemeines Wandstirkegesetz. In diesem Abschnitt soll ein allgemeines Wand- 
starkegesetz nach Art der Gleichung (17) angenommen werden. Wird dann der gesamte 
Kreiszylinder mit einem einzigen Teilzylinder identifiziert (das ist der Fall fiir m= 1), 
so hat man die Beulbedingung 


et (1) 
KY KY 

1 1) 
KQ KG 


1 1) 
Ay? AY 


1 1 
AY Ag 


=i (38) 


Setzt man die Elemente gemaf (27) ein, so erhalt man, wenn man den nunmehr tiber- 
fliissigen Index 1 weglaBt, fiir (38) 


sina: sinb=0. (39) 
Hieraus folgt 
b baw. a=ka (k ganz) (40) 
und fiir beide findet man unter Beniitzung der Gleichungen (22) 
y 2 + kA a4 
Pai Be — (41 ) 
und wegen (20a) schlieBlich 
Eh, B2 + k4 oct 
Vie 12 (1 —»2) 2 k2 a2 = 9 (k). (42) 


In bekannter Weise1? hat man das Minimum von (42) zu suchen. Es ergibt sich fiir 
piss 
kyo = Ve. (43) 


Voraussetzungsgemals muBte ky ganzzahlig sein, was man jedoch im allgemeinen 
fiir (43) nicht erwarten darf. Man wird deshalb in (42) die zu k, nachstgelegene kleinere 
bzw. gréBere ganze Zahl k,’ bzw. ky’ +1 einsetzen, je nachdem gilt 

Koi Vee ee) (44a) 


bzw. 


ky > V ko’ (ko’ + 1). (44b) 
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Diese beiden Bedingungen (44) erhalt man aus den Annahmen 


p (ko’ + 1) : 
wenn man y nach Gleichung (42) einsetzt. 
Fiir die kritische Axiallast gilt damit sicher 


Pe > 7 (k yi 
Die Abweichungen 
G (ko) —¢ (ko’), 
bzw. ee Gare 
@ (ko) —@ (ky +1) 


sind meist unbedeutend, so daB es geniigt, wenn man q (k,) als kritische Axiallast P,, 
ansieht, mit welchem Wert man ohnehin auf der sicheren Seite liegt. 
Somit erhalt man aus (42) fiir k=k, 


. ER, : fe 
= BI pe 2P (45) 


und mit. der Bedeutung fiir 6 nach (20b) sowie jener fiir ha; nach (17) schlieBlich 


ES = (EEO ey (46) 
— 72) 


Pe 8a 


Fiir «=0 (konstante Wandstarke) wird (46) zu einer bekannten Formel. 


2. Spezielle Wandstarkegeseize. Ist irgendein Wandstirkegesetz vorgegeben und 
ist die Bedingung (16) hierfir als erfillt anzusehen, so hat man entsprechend der 
Gleichung (18) nur den Wert F zu berechnen und erhalt in (46) sofort den kritischen 
Axialdruck. 


a) Die Wandstirke als Potenzgesetz. Dieser naheliegende Fall ist gegeben durch 
das Gesetz 
hey=h (Uwe) Ost al. (47) 


Nach (18) errechnet sich hierzu 
i 
ne +1? 


und der kritische Axialdruck P;, ist dann 


F 


Eh? is uy a it 

ie eae (tage) (48) 
Bei linearer Wandstirkeanderung (n= 1) resultiert aus (48) genau die von Feder- 
hofer?’ angegebene Formel, waihrend im Falle quadratischer Anderung (n= 2) gegen- 
iiber Federhofer eine Abweichung besteht. An Stelle des nun in (48) stehenden 
Faktors 2/3 vor uw erhielt Federhofer den genaueren Wert 


2 5 rho 


peo Ms (Tey ee 


(49) 


Von dem zweiten Term in (49) kann man aber zeigen, daB er in vielen Fallen 
gegeniiber 2/3 verschwindend klein anzunehmen ist: . 
Die errechneten Werte fiir den kritischen Axialdruck gelten ohnehin nur im elasti- 

kr 


= oR ee 
schen Bereich, so daB man mit 7 > unterhalb der Proportionalitatsgrenze des ver- 
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wendeten Stahls liegen mu. Fiir Baustahl St 37 ergibt sich z. B. die Bedingung’ 


ho 
1 = 1OoOM 


(50) 

Ist diese Bedingung nicht erfiillt, so tritt vor der Gefahr eines axialsymmetrischen 
Ausbeulens des Zylinders Uberschreiten der Elastizitatsgrenzen ein, womit sich die 
Verhiltnisse vollkommen andern. Die vorstehenden Ausdriicke fiir den kritischen 
Axialdruck haben dann keine Giiltigkeit mehr. Der Kreiszylinder muB also wegen (50) 


sehr diinnwandig sein*. 
Andererseits darf aber das Verhaltnis 7 nicht beliebig groB werden, weil sonst 


der Faktor (49) kleiner als Null werden kénnte, wodurch in der Gleichung von Feder- 
hofer fiir die Beullast bei quadratisch zunehmender Wandstarke ein kleinerer Wert 
herauskime, als es bei konstanter (Anfangs-) Wandstirke der Fall ware, wo ja die 
Rohrlange nach Gleichung (46) fiir ~=0 keine Rolle spielt. 

Fernerhin zeigt die Erfahrung, daB im Falle /> 2,7 (also fiir lange und tiberdies 
diinne Rohre) die Kreiszylinderschale wie ein Stab in einem Bogen ausknickt. Theo- 
retische Schranken fiir 27/] sind hierfiir angegeben worden’. 

Eine axialsymmetrische wellenférmige Ausbeulung kann demnach nur auftreten, 
wenn / und 27 ungefahr von der gleichen GréBenordnung sind. Beachtet man dies, 
so kann der zweite Term in (49) in der Tat gegen 2/3 verschwindend klein werden, 


5 ics 
2V3(1—») P 


Dy\2 
0,006 (7 . (51) 


Bei Anwendungen wird man also stets zu untersuchen haben, ob (51) tatsachlich — 
gegen 2/3 vernachlassigbar klein ist. In diesen Fallen liefert dann die Gleichung (48) 
fiir Wandstirkegesetze in Potenzform geniigend genaue Ergebnisse. 


b) Die Wandstérke als Sinusgesetz. Fiir ein Wandstarkegesetz 


h()=ho (1+ sin? > §) (52) 


ergibt sich nach (18) 
Pes V2 


und (46) wird zu demselben Wert wie bei linearer Wandstairkeainderung. 

Nachdem es sich gezeigt hat, daB mit den nétigen Einschrankungen innerhalb 
eines Kreiszylinders brauchbare Naherungswerte fiir die kritischen Axiallasten zu 
erzielen sind, wird man nun einen Kreiszylinder mit einer nicht mehr geringen Wand- 
starkednderung (d.h. wenn fiir den gesamten Kreiszylinder die Forderung py? < 1 
nicht mehr als erfiillt angesehen werden kann) in m solche Teilzylinder zerlegen, fiir 
welche einzeln wu? <1 (i=1, 2, ..., m) giiltig ist. Die Gesamtheit dieser m Teilzylinder 
reprasentiert dann das Ersatzsystem, auf welches das Reduktionsverfahren anzu- 
wenden ist. 


VI. Die Falle m>1 


Das Hauptproblem besteht darin, die Beulgleichung (36) zu lésen. Fiir den Sonder- 
fall m=1 ergab sich hierfiir die einfache Beziehung (39). Wesentlich komplizierter 
gestaltet sich die explizite Formulierung der Bedingung (36) bei mehreren Teilzylindern. 


* Das Beulen von Kreiszylinderschalen im plastischen Gebiet hat Geckeler! erstmals be- 
handelt, indem er von den Deformationen der Schale vor Erreichen des kritischen Wertes ausging 
und dann die weiteren Deformationen bei allmahlich zunehmender Belastung verfolgte. Dabei 
wurde an Stelle des Moduls # der sogenannte Knickmodul K eingesetzt. 
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Wie man dann zu verfahren hat, soll nach der Besprechung des Falles m= 2 angedeutet 


werden. 


; Lie Der Fall m2. Um die Beulgleichung (36) explizit anzugeben, berechnet man 
sich die hierfiir bendtigten Elemente der Beulmatrix 


R= A, U, A. 


(53) 


(Wie Ul und % hat auch & in allen Elementen die Dimension 1.) Bezeichnet man allge- 


mein noch die Ubergangsmatrix 


U; wie folgt 


UOT POO” PGA SEG 
0 a 0 1) 
U1, == A ? 
0 0 uy 0 
GPS wgulria re quliky gh 


(54) 


worin die u (9 =1, ..., 4) aus (31) abgelesen werden kénnen, so ist (unter Weglassung 


des oberen Index | fiir 11,) 


1 1 1) 1 
Uy A‘? ab, ASS Uy Att Uy AW) 
U, x Us Aly Us As Us El Us As? (55) 
eS } 
Ub AY Us AY Us Ay Uz Alp 
TAG ug A a A ray AD 
Alsdann ist nacheinander 
\ (2) 4) | 
ea) Ane 2 if 
Ki? _— Ay; Aj2 Uj 
irs 
4 
ns Rhee: 1) 
Kye = 2) As) Ape wy 
Eh (56) 
KR? = YAP Ay 
j=1 , 
1 
2 eal 2 4 
KO = 3'AP 48 wy 
ee | 
sowie 
4 
K@ KD = D' AP AS) AP? Ale wy ux (57) 
ds =1 
und 
4 
2 2) 2) 2 (1) CL) ea 
KG K2 = & Ay; As) Aji Ajs Uj Uk. (58) 
fk= 
Aus (57) und (58) hat man den Ausdruck zu bilden 
P P 2) z(2 
A=KO KY — Ki KR, 
also 
A 
~ < Fy 1 
A= 3) A® AQ [AD AR? — Af? Ap?) wu. (59) 


jek=t 
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Vertauscht man in dieser Doppelsumme j mit &, wodurch ihr Wert nicht geandert 
wird, so folgt 
4 . 
A=— SAP AQ [AW AQ — AP AD) a a. (60) 
jsk=1 
Wegen A=0 erhalt man schlieBlich durch Addition von (59) und (60) die explizite 


Beulbedingung 
4 


: 2 2 2 2 1 1) (1) 4(Q) 
D>) Uj UE ay 7 Aas [Af Au Ee Aj, Ajs = (61) 
jk=l 


Die Summation ist in (61) nur noch iiber k <j zu erstrecken. Die anderen Glieder sind 
teils Null (fiir k= 7), teils treten sie doppelt auf (fiir k > 7). So hat man an Stelle von 
4-4—16 Summanden tatsdchlich nur noch deren sechs in Gleichung (61). Ein aus- 
fiihrliches Anschreiben dieser Summe, wobei fiir die AW die entsprechenden Werte 
nach den Gleichungen (27) zu stehen kommen, mag unterbleiben, weil diese Formel 
sehr umfangreich ist. Es sei nur mitgeteilt, da der Fall 


Oy == == & 


aun 
(meee | 


zur Verifikation nachgerechnet wurde. Dieser Fall bedeutet, da8 die gesamte Kreis- 
zylinderschale in zwei vollig gleiche Teilschalen [also auch hiy=hoy und u,=1 
(o=1, ...,4)] zerlegt wurde und daher somit die Transformation unndtigerweise 
tiber eine Verbindungsstelle in der Mitte durchgefiihrt wird. Dabei reduzierte sich 
notwendig die Gleichung (61) nach einiger Rechnung auf 


(62) 


sina: sinb=0, 


wie das auch nach (39) sein mu, weil sich dieses zweiteilige Ersatzsystem genau so 
gut wie ein ganzer Kreiszylinder hatte behandeln lassen*. 


2. Numerische Behandlung der Aufgabe bei beliebigem m. Es hat sich schon bei 
der Betrachtung des Falles m= 2 erwiesen, daB es gar nicht zweckmafig sein wird, 
explizite Beulbedingungen anzustreben. Vielmehr wird man in praxi so vorgehen: 

Fiir eine vorgegebene und in m Teilzylinder aufgeteilte Kreiszylinderschale wird 
man die Ubertragungsmatrizen Y%; zuniichst einmal berechnen, wobei man sich der 
Rechenkontrolle 

det Yl; ==) (63) 
bedienen kann. Dabei hat man aber in den Matrizenelementen AY fiir P zunachst 
einen festen Wert einzusetzen, welchen man sinnvollerweise nach jenem Naherungs- 
wert P,, orientiert, den man nach Gleichung (46) fiir den gesamten Kreiszylinder 
erhalt. Was hierunter verstanden wird, ist im Sinne der Mechanik anschaulich klar: 
Einmal setzt man in Gleichung (46) fiir die konstante Wandstarke (d. h. ~=0) den 
groBten Wert hax, dann den kleinsten Wert hin der gesamten Kreiszylinderschale 
ein und erhalt so Schranken fiir den gesuchten kritischen Axialdruck**. Die festen 
Werte fiir P wird man also aus dem Bereich innerhalb dieser beiden Schranken wahlen. 


* Eine andere, notwendig richtige Kontrolle fiir (61) wire die Wahl von L=l und = 0 
gewesen. In diesem Fall ist die Rechnung aber so trivial, da8 eigentlich von einer Kontrolle nicht 
mehr die Rede sein kann, weil fast alle Glieder in (61) a priori Null sind. 

** Der mathematische Beweis dieses Sachverhalts ist im Anhang (VII) gegeben. 
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Danach wird man nicht etwa die Produktmatrix ®t, explizit berechnen, weil 
man namlich hierfiir unvorteilhafterweise eine Reihe unnotiger Operationen durch- 
fihren miBte, sondern man wird zuerst den transformierten Vektor Y, wy durch 
Multiplikation der Matrix 2, mit dem Vektor ww, bilden. Damit hat man dann praktisch 
auch den Vektor MU, 2, w, und kann aus diesem durch Multiplikation mit der folgenden 
Matrix %, den Zustandsvektor 1, berechnen usf. Man arbeitet also mit der jeweiligen 
Matrix nur an einem Vektor und spart so erheblich an Operationen (9, S. 19). 

Auf diese Weise erhalt man schlieBlich die Elemente der Beulmatrix &,,, welche 
man zur Bildung der Beulgleichung (36) benotigt. Diese Gleichung hat man als Funktion 
von P aufzufassen. Man triagt A (P) iiber verschiedene P auf und bestimmt durch 
Interpolation den gesuchten kritischen Wert P;,, fiir den zum erstenmal A —0 wird. 
DaB dieses Vorgehen tatsichlich die kleinste Nullstelle P,, fiir A (P) liefert, folgt aus 
der Tatsache, daf die erste Niherung nach Gleichung (46) (hy =Amin und yu = 0 gesetzt) 
zu klein sein wird. 

Die Hauptarbeit bei groBem m ist also (neben dem Aufstellen der Ubertragungs- 
matrizen selber) die Multiplikation dieser Matrizen und die Auflosung der Beul- 
gleichung (36). Derartige Rechenarbeiten kénnen aber von einer GroSrechenanlage 
geleistet werden, wenn viele solcher Rechnungen eine Programmierung und Rechnung 
auf der Rechenanlage wirtschaftlicher erscheinen lassen als Rechnungen an der Tisch- 
rechenmaschine. 


VII. Anhang 


In diesem Abschnitt soll zunaichst der Kinflu8 des Terms y? fiir den kritischen 
Wert P, bei quadratischem Wandstarkegesetz abgeschitzt werden, auBerdem soll 
noch eine Anwendung fiir einen anderen Belastungsfall der Kreiszylinderschale erwahnt 
werden. Beide Male wird von der Energiegleichung ausgegangen. AbschlieBend wird 
ein Vergleichungssatz fiir den ersten Eigenwert formuliert und bewiesen. 


1. Ritz-Verfahren fiir quadratisches Wandstarkegesetz. Hier wird nun die Variations- 
gleichung (7) des elastischen Potentials // fiir eine Naherungslosung bentitzt. Bedeuten . 
jetzt [in Abainderung zu den fritheren Abkiirzungen (20)] 


12 (1 — 92) 2 


Pye B (64) 
ft = _ 2) i is 
so lat sich (7) in folgender Weise formulieren (e= = 
1 
We oft (£) w’”? (£) — a w’? (£) + Bt (é) w? (é)] d = Extr. (66) 


0 
mit © als einer fiir das Folgende unwichtigen Konstanten. 
An Stelle der gesuchten (und genauen) Lisung w (&) wird in (66) eine Vergleichs- 
funktion w* (£) verwendet in Form eines eingliedrigen Ritz-Ansatzes 


w* (E)= A - v (é), (67) 


worin A +0 eine Konstante darstellt. Von diesem v (£) wird lediglich die Erfiillung 
der wesentlichen (oder geometrischen) und restlichen (oder dynamischen) Rand- 
bedingungen (19a) verlangt. Die Naherungsrechnung wird jedoch besser sein, wenn 
vy (€) auch gewisse Eigenschaften der genauen Lésung w (é) besitzt. Im vorliegenden 


Fall setzt man deshalb an 
v (€é)=sink a €. (68) 
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Mit (67) wird damit das Variationsproblem (66) zu einer einfachen Extremumsaufgabe 
; d 
fiir J] (A). Der gesuchte Eigenwert « ist dann wegen aa == () 


if 1 
pe | h (8) v8 (6) dé + | h(E) 02 (6) dE 
Fie a : é (69) 


Fiir die Wandstarkeanderung wird ein quadratisches Gesetz gewahlt 
h (€)=ho (1+ 4p €), (70) 


in allen folgenden Rechnungen wird aber yw? vernachlassigt. Die elementare Rechnung 
fiihrt auf einen langeren rationalen Ausdruck fiir den Eigenwert «, 


=a. (2), (71) 
wenn man substituiert z—k? 2. Wie friiher wird auch hier das Minimum von « (z) 


als kritischer Wert angesehen. Mit os =0 ergibt sich dann die kubische Gleichung 


) 2 3 2 1 9 2 
z (i+ hy— ni) 2 (B2 a ee he) + p—«, (72) 
Zur Losung dieser Gleichung (72) wird ein Stoéransatz 
= ky 2+ Ul? 2p (73) 


in (72) eingefiihrt. Das ungestérte Problem (u = 0) ist sofort losbar, so daB die Anwen- 
dung der Storungsrechnung zweckmabig erscheint. Sieht man die ,,ungestdérte‘ Losung 
als Naherung an, so ist die Losung auch mittels des Newtonschen Iterationsverfahrens 
zu erhalten. 


Somit ergibt sich als Lésung von (72) 
B Tee on 26) oak 15 he % 
a bfimolh +B) ead Bh oe 


Setzt man diesen Wert in (71) ein und entwickelt man nach Potenzen von y, so resultiert 
als Naherung fiir den kritischen Axialdruck 


Bh 2) mn hor 
7 V3 (1 — »2) ( 2/3(1—») PB 


lee 


kr 


| : (75) 


11 hor 17) Br 
te SSS : 
45: 123 (= payee Coen EG 


Fir «=0 liefert der Ritz-Ansatz denselben Wert wie die Differentialgleichung 
(weil naémlich dann v ihre Lésung darstellt), wihrend im Falle eines quadratischen 
Wandstarkegesetzes der Faktor bei w gréBer ist als jener, den Federhofer’ mit Hilfe 
der Sto6rungstheorie fiir das Kigenwertproblem erhielt. Diese letzte Tatsache ist durch 
die Minimaleigenschaft des exakten Eigenwerts gegentiber den Rayleighschen 
Quotienten (69) belegt. 


2. Die in ein elastisches Medium eingebettete, axial gedriickte Kreiszylinderschale. 
Bei der Untersuchung von gerammten Pfaihlen in Form von Kreiszylinderschalen hat 
man es in der Regel mit einem Baugrund zu tun, dessen Bettungsziffer c (w) mit der ~ 
Rammtiefe x zunimmt. 
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Die hierfiir geltende Differentialgleichung 148t sich aus (7) herleiten, wenn man 


den Ausdruck (6) noch ergainzt durch den Arbeitsanteil jener auBeren Kraft, welche 
der radialen Verschiebung w entgegenwirkt 


I 


1 
a ze (a2) w? (a) dee. (76) 


0 
Damit lautet die Eulersche Gleichung des so formulierten Variationsproblems 


d2 Eh? (x) ss b (a 
dae [12d am" (2) +P w(x) + {AO +c (a)hw (a) =0. (77) 


In der Praxis geniigt oft die Annahme einer schwachen linearen Veranderlichkeit 
der Bettungsziffer 


c=e,(1+6£) mit =F, (78) 


wenn mit ¢, die Bettungsziffer am oberen Rande der Schale (= 0) und mit 6 die relative 
Anderung derselben bezeichnet werden. Definiert man analog (17) auch fiir c den 
linearen Mittelwert 


Ei a 6) (79) 


und verwendet man iiberdies fiir h den Wert hy, so fiihrt die Losung der Differential- 
gleichung (77) mit den friitheren Randbedingungen (19a) auf die kritische Axiallast 


Bh* j 2 1 
De Pe jon Sa) a ae ( ) 
r\3 (1 — v2) | A Eh 1+ 3 9). 80) 


Fiir konstantes h =A, ist dies der von Federhofer™ angegebene Wert. Man erkennt 
in (80) die versteifende Wirkung einer elastischen Bettung. 

Mit dem Reduktionsverfahren ware nun prinzipiell die Moglichkeit gegeben, den 
Fall veranderlicher Wandstarke und veranderlicher Bettungsziffer zu behandeln. 


3. Ein Vergleichungssatz. In diesem letzten Abschnitt soll ein Satz bewiesen 
werden, der — anschaulich betrachtet — eigentlich eine Selbstverstandlichkeit ist 
a ALL): 
ae bei einem Stabilitatsproblem an irgendwelchen Stellen des Systems die 
Steifigkeit erhdht und nirgends erniedrigt, so konnen sémtliche Kigenwerte nur gréBer 
oder wenigstens nicht kleiner werden. Das Umgekehrte gilt bei Erniedrigung der 
Steifigkeit.“ , 

Zunaichst muB dieser Satz abstrakt formuliert werden in einer Weise, die dem 
Sachverhalt in dieser Arbeit Rechnung tragt. 

Vorgelegt sei das homogene selbstadjungierte volldetinite EKigenwertproblem der 
Hingliedklasse* 


M [yj=4N [y] (81) 
mit den Randbedingungen 
U,[y]J=9 = (@=1, 2, 3, 4). (82) 
Im einzelnen sei darin 
M [y]= {A (x) yy" («)}" + B (x) - y (2), (83) 


* Die Nomenklatur ist hier dieselbe wie in ¥. 
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N [yJ=—y" (#) j (84) 
und 

U, [y]J=y (0), U; [yJ=y () (85) 

Us lylJ=y" (9), U, [yJ=y" (. (86) 


(85) werden die wesentlichen (oder geometrischen) und (86) die restlichen (oder dynami- 
schen) Randbedingungen genannt. 
Die Bedingungen der Selbstadjungiertheit 


1 
| {u M [v])—v M [u]} da=0 (87) 
0 
und 
i 


| {uN [v]—v N [u]} de=0 (88) 
: 
fiir zwei beliebige Vergleichsfunktionen u, v sowie die Bedingungen der Voll- 
definitheit 
| wM [ujde >0 (89) 
; 


| uN [u]dx>0 (90) 


sind dann erfiillt, wie man durch partielle Integration nachweisen kann. 

Bei entsprechender Deutung der auftretenden Grof8en stellt dann die Gleichung (81) 
genau das Eigenwertproblem (9) dar. 

Im Intervall [0, 7] gilt auBerdem A (x) +0, B(x) +0. Das erstere deshalb, damit 
die Existenz einer Losung der Differentialgleichung (81) gesichert werden kann. 

Dann folgt daraus fiir den Rayleighschen Quotienten 


l 
[uw M [uj dx 


R (u]=*;}———_ > A. (91) 
[uN (ul da 
0 


Verwendet man als uw die erste Eigenfunktion y,, so ist R [y,]=4,. Die erste Eigen- 
funktion y, lost also die Variationsaufgabe, den Rayleighschen Quotienten zum 
Minimum zu machen, wenn u den Bereich aller Vergleichsfunktionen durchlauft. Der 
zugehorige Minimumwert ist der erste Eigenwert /,. 


Nach diesen Vorbemerkungen folgt nun der abstrakte Vergleichungssatz: 


a) Es sei 
A («)= Ag+ A, (2) 
(92) 
B («) = By + B, (2) 
mit 
A,=Min A(x) >0 
x € [9,1] 
: (93) 


By = Min? BB (e)S10 


x € [0,1] 


> 
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also 


A,(z)20, B > 0 in [0, 1]. 
Bezeichnet dann 1 (@) aes pak (94) 


My [y]= (Ao y+ By y, (95) 
so gilt fiir den ersten Eigenwert 2, des Problems M [yJ=AN [y] 
y Digital (96) 


wenn unter 4,* der erste Eigenwert des Problems M, [y]=/ N [y] verstanden wird. 
b) Entsprechend folgt fiir 


A (%)= A,— A, (2) 
3 (97) 
, B (x) = B,— B, (x) 
mit 
Wee Max4 Ai(z)mee 0 | 
x € (9,2) 
Be Max be) | me) 
x € [0,1] 
also 
A, (x) =0, B,(x)=0 in [0,1] 
und 
M, [yJ=(Azy")" + Bry (99) 
fiir das EKigenwertproblem 
M, [y]=2N [y] (100) 
mit dem zugehorigen ersten Eigenwert /1,** die Behauptung: 
Ay x A,**, (101) 
insgesamt also die EinschlieBung 
Ase Apesdlie™ (102) 


Beweis: Es wird nur Behauptung (96) bewiesen; der Beweis fiir (101) verlauft 
dann analog. 
Wegen der Linearitat des Operators M und wegen (83), (92), (95) gilt 


M (y)=M, [y]+ My] (103) 

mit 

M* [y]= {Ay (x) yf" + B, (2) y- (104) 
Damit hat man fiir den Rayleighschen Quotienten (91) 
l l 
Jw Mo (u) dx +fu M* [ul dx 
A, = —— (105) 

fu [u] da 


v0 


und mit uw=4y, 
l I 


[ys Mo Cyd de + fy. M* Lys] de 
pues : (106) 


q 
[aN (yi) de 
0 
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also 
I 


fy M* Cyn] de 
42a (107) 
[ 41 N [ys] dx 
b 


Nun ist aber fiir alle Vergleichsfunktionen wu 


l 


jw M* |u| dee (108) 
0 
ui 
| uN [uj dx> 0, (109) 


ty 
denn fiir (108) folgt durch partielle Integration 
l 


ju {(A, (x) u'’)"" + B, (a) u} da= 
Te: 1 


== UlAray | = [ aw’ (A, w’)’ dx + | B, wu dx 


o 0 0 
5 
l l ei 
=—wu' Ayu" + A,u d+ | B,wdx 
Oo ny 
ennteneananeeemaaanen 
0 


und wegen (94) gilt dann (108). Da® die Gleichung (109) zu Recht besteht, ist eine 
Folge der Volldefinitheit des vorliegenden Problems [s. Gleichung (90)]. Damit ist 
der Teil 1,* <A, bewiesen. 
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Die Verfahren von Ritz und Trefftz in der Reissnerschen 
Plattentheorie 


Von D. Riidiger, Freiberg/Sachsen 


Zusammenfassung. Es wird ein Variationsprinzip aufyestellt, das der dritten Randwertaufgabe 
in der Reissnerschen Plattentheorie aquivalent ist. Randbedingungen treten dabei nicht auf. 
Deshalb kénnen zur Lésung die direkten Methoden der Variationsrechnung ohne Schwierig- 
keiten angewendet werden. 


I. Einleitung 


Die Gedanken, die zur Verallgemeinerung des Prinzips vom Minimum der poten- 
tiellen Energie elastischer Korper fiithren!, werden auf die Plattentheorie unter Beriick- 
sichtigung der Querkraftverformungen von E. Reissner? tibertragen. Randbedin- 
gungen des Variationsproblems treten nicht auf, so dafi an die bei den Verfahren von 
W. Ritzund E. Trefftz verwendeten Funktionen hinsichtlich der Randbedingungen 
keine Forderungen gestellt zu werden brauchen. 

Zur Ableitung des Variationsprinzips werden die Grundgleichungen der Reissner- 
schen Plattentheorie in allgemeinen ebenen Koordinaten® 


Q|.+ p=0, M4, QF=0, 
3 1 5 B 
M = aw 5 HOSP tae eee 8° (Bylo Day) 1s (1) 
ER  5( 


a 1 =) (ows WV \% 
C= RI) frail 


verwendet. 


Der Querkraftvektor Q* und der Momententensor M* 8 legen die Randquerkraft Q 
und die beiden Randmomente 7’ und B fest*. Sind e* die kontravarianten Koordinaten 


1 iidiger: Ing.-Arch. 29, (1960). 

: a eee J. Math. pies. 184 (1944); J. Appl. Mech. 12, A 68 (1945); Quart. Appl. 
Math. 5, 55 (1947). Vgl. hiezu auch: A. E. Green: Quart. Appl. Math. 7, 223 (1949); M. Schafer: 
Z. angew. Math. Mech. 32, 161 (1952); A. Kromm: Ing.-Arch. 21, 267 (1953). : 

3 —D. Rudiger: Ing.-Arch. 28, Festschrift Richard Grammel, 281 (1959). Die Grundgleichungen 
der Reissnerschen Plattentheorie in allgemeinen Koordinaten folgen aus den Gln. (31) und (34) 
fiir verschwindende Koeffizienten der zweiten Fundamentalform der Flachentheorie. 

4 H. Neuber: Z. angew. Math. Mech. 29, 144 (1949). 
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des Tangentenvektors der Randkurve C, dann gilt mit dem Normalenvektor n= ¢q5 eb 
fiir die Randquerkraft 
Q = n,Q" (21) 


und mit dem Randbiegemoment B=n, n, M* ® und dem Randtorsionsmoment 7’ = 
= —n, es M** fiir den Momentvektor des Plattenrandes 


Me =e mn, MVS BP +T nh. (22) 


II. Das Variationsprinzip 
Eine elastische Platte mit der Flache / und der sib tices acct in p ist beliebigen 
Randbedingungen so unterworfen, dafs an den Randstiicken c, ¢, ¢ die Pande 


das Randtorsionsmoment 7’, das Randbiegemoment B und an den Randstiicken c, c, ¢ 
die Bandyerecnohuns ue und die Bepdgercuhinern e* },, n* 0, vorgegeben sind. Es 


ist dann C=c +c, C=c uae c und C=c-+c. Zur Formulierung des Variationsprinzips 
wird die Formanderungsarbeit —6 // der im Gleichgewicht befindlichen Platte bei dem 
Verschiebungszustand 6 W, 60% gebildet. Sie ist mit dem Flachenelement d/ 


— 61 ={ [(Q%\at+p) 6 Waf +| | (A""|, —Q°) 69, df =0. (3,) 
F F 


Fiihrt man die Identitaten 
M8), 694 = (M*? 68,)|, — M8 dds, 
= (Q* 6W)|, —Q* OW, 
ein, dann ergibt sich 
—di=—f [Me d9,),d¢ — f [Q°0(8, + W]e) df + 
F F 
(32) 
+] (Me? 69)\.df + [ [(Q* OW). af + | fpsWaf =o. ; 
F F 'F 


Fiir die beiden ersten Flachenintegrale in (3,) kann mit dem Elastizitaitsgesetz (1) 


| fare? oa) af— [ [Q28, + Wh.)af = > } | 5(M=® dp),) df + 


+ 3 | [ole (0.4 Wada 
E 


geschrieben werden, so dai (3,) zu 


il 


— df= — 5 | | date? 4.) df — 3 [00 (+ Wle)df-+ 
F F 


+ | [ (ate? 69,)/a¢ +f [Qt oW)|a af +] [pow df=o 
ir PF rf 


folgt. Beachtet man noch die wegen (1) giiltigen Beziehungen 
yer Dela == (Mee Os) la = Qe UF > 
Q* Wie (Q* Wie oe Pp W; 


ENS 
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dann erhilt man als potentielle Energie 


d= | | OM? 9,) af +5] fe 5(Q° W) df — 


F 
— {fe (Me? 59,)\, df — tele (Q* 6W)|,.df + 


F FP 
ire : ‘ : 
+5 |p Wyaf — | |powaf=o. 

“le . m. 
Die Anwendung des GauBschen Integralsatzes 


tye Ja df= o ( )n, ds 


F (0; 


auf die ersten vier Integrale in (4,) ergibt 


Alm 2 
dM => $ d(M** d,) nds +5 § 0(Q* W) nq ds — 
C C 
= Mee 60, n, ds — bQ ow N, as + 


Cc 
+3 | [ow W) df — | {pow df= 0, 
Sh 


F 


und nach eis der Randquerkraft und der Randmomente (2) wird 
61 => $ 6(Q W) ds — zo d(T e 9.) ds +5 : (Bn 6.) ds — 
C 
— QW ds + Cae o,)ds — Pais ) ds 
C 
toe J fo W) df — { | pow df =0 
F 

erhalten. Werden die Randintegrale in (43) durch die Linienintegrale lang 3 


dp ee 3 
6, c, & ersetzt, so folgt 


te 
sIT = 700 W) ds — | (Pe 9,) ds + 7 | (Bn 9.) ds + 
> 3 
+5[o@W) ds — = (Te°d,)ds +5 [o(B ns 0,) de — 
1 2 : 


= ie bWds + | 1 6(c% 94) ds — | Bon 8.) is — 


— (QoWds + [rae ,)ds — 


! 


& 
Se 
S 
2 © || 
S 
ce 
a 
iva) 
ee 
—————=—_E_O———e a 


2 
SC, C, 
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Fiihrt man in (4,) die Identitaten 
—QdW=—<d(Q W) +60 W, 
T d(e* 3,)= 6(T e* 0,) — OT e* 8, 
— B d(n* 0,) = — 6(B n* 8.) + OB n* 8, 
ein, dann wird zunichst 


dH =+[8QW)ds — 5 for e*0,)ds-+ +] (Bn° 9.) ds+ 


| mee. 


oe 
_——, ad 


& 
|o(QW) ds — +f 0 Te 9,) ds + >| (Bn 9.) ds— 
2 3 


¢ 


6 |[ € 


Bo(n* 0,.)ds— 


. 
2 


— [Qowds + | Ta(e" 8.) ds— 
_ 


o [ot 


c 


a | ae W) ds 4 [ee 8.) ds — jie n° 9.) ds 


¢€ C 


+ foow ds — fore d, det Ba 9,ds + 


c ¢ 


+3 | fow wa [ [pow a=o 
F 


F 


und nach Zusammenfassung 


oll =5Ja@ W) ds — 5 fares 8) ds + > [o(B ns 0.) ds — 
. 3 ay 
1 a a ! 
—+] 3(Q W)ds + x | (Le 8.) ds — yf (8 n* 3.) ds — 
2 2 = 
— |QoWds + [roe 9.) ds — | Bd(n* 8.) ds + (5) 
i = 


| 30 Was — | OT & 8.ds + lea 9,,.ds + 


+- 


c 


+] fom may — [ (pow df=o. 
F 


5 |[ een, & | 1 


Die Beziehung (5) ist die modifizierte Energiehauptformel, aus der sich sofort das 
Prinzip vom Minimum der potentiellen ate in der Soe Plattentheorie 


herleiten 1a8t. Da namlich in den Pa a . c wae ¢ die Randquerkraft Q und die 


Randmomente 7, B, in den Randstiicken ¢ = t F und = die Randverschiebung Ww und die 


ee 


: 
; 
; 
4 
; 
3 
4 
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Randverdrehungen e* 8. n* é, und auf der Plattenmittelflache die Flachenkraft p fest 
vorgegeben sind, gilt 


2 


c 


1 ( aon aee 1 = : om 
6M =5|5 | (Q—2Q) Wds— 5 (7-27) ce 0,ds4+1{(B—2 B)n* Igds — 
1 : 3 


2 
c 


1 od l a z 
—5[o(w—2I )ds+ x | Tet (8,2 8,)ds— [ Bn (8 —286,) ds + 


1 


1 
+3) {@—-2n Wal =o. 
F 


Hieraus ergibt sich, daB unter allen moglichen Verschiebungen einer Platte, die in 


allen Randpunkten beliebige Werte annehmen k6énnen, diejenigen ein- 
treten, die den Ausdruck 


zu einem Extremum werden lassen. Das Extremalprinzip (6) hat gegeniiber der klas- 
sischen Formulierung den Vorteil, daB an die infinitesimal benachbarten Systeme von 
Verschiebungen hinsichtlich der Randbedingungen keine Forderungen gestellt zu 
werden brauchen. Man kann also zur numerischen Darstellung der Lésung W= 
W (a1, x), 0* = 0* (x1, x) eines vorgegebenen Randwertproblems der Reissnerschen 
Plattentheorie im Sinne von W. Ritz einen Ansatz 


oe— S'e ge, W= Yeo (A) 


' v v A 2 2 : é 
mit willkiirlichen Funktionen g*, m machen. Die Funktionen sind dabei nur so allgemein 


auszuwahlen, daB bei hinreichend groBem n und passend gewahiten Koeffizienten’c. ¢ 
die Lésung der vorgegebenen Randwertaufgabe mit beliebiger Genauigkeit dargestelit 
werden kann. 

Das System von Verformungen (7,) liefert fiir den Momententensor und den Quer- 


kraftvektor 


n E h2 es 2 ¥ al vB Vv el x viv _ vo v | Vl 
MP =m DOO +e OP + vert h(C Gy + C Hohl, 
v= 
(73) 
sf Bie o(l—») a ae reed 
C=piI-*) = \° pt +60) 
und fiir die RandgréBen und Flachenbelastung 
n n n n- n n n n 
=n, Q, T——n, ¢, M*®, B=n, ng U**, p= —Q"*\,.- (73) 
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Die Koeffizienten ‘C. c: y=l, 2, ..., n,a=1, 2 des Ansatzes (7,) werden nun so be- 
stimmt, daB die mit (7) gebildete potentielle Energie 


T= 5| ig ae 0) W ds = Ca 2 T) oo B,de4+ 4 (B 2B) n*d,ds — 
is Ey = 
= ak (W— 2W?) ds + ie e«(3, — 29,) de — 5 Bn(9, — 280) ds (8) 
3 = us 


dem Extremum méglichst nahe kommt. Die Bedingungsgleichungen fiir das Eintreten 
des Extremwertes lauten 


n n 
OID sag CFE. 9 peed Seed Aan aS (9) 

v,@ Vv 

0 ¢ 0¢ 


und haben die Galerkinsche Form. 


Ill. Das Verfahren von Trefftz 


Beim Trefftzschen Verfahren wird der Ansatz 


2 oy N; x ny Vv 
OO er Ps Wi 25m (10) 
v= v— 


verwendet, wobei die Funktionen ¢*, @ Partikularlésungen der homogenen Platten- 
grundgleichungen (1) sind und entsprechend dem Variationsprinzip (6) in den Rand- 
punkten beliebige Werte annehmen kénnen. Bildet man die dem Verschiebungszu- 


n n 


stand (10) entsprechenden RandgréBen Q, 7’ und B, so ergeben sich mit 


Vv 
m 


me Eh WT Seals Vole a iv Vv 
= Baw) cal P OP taame ¥ £65(5,14 + Gor) 


Vs Eh3 5 (l= 9) (ge 
SAYIN OFT p Ney aS 


Vv 
Vv Vv vv ah 
b=ngngm?, t=—n em G=n.g 


die Ausdriicke 


(115) 
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Die n linearen Gleichungen fiir die Konstanten lauten 


— 


| 
to) 
 |J AL s 
Re 
Ss 
-+ 
Os 
Se 
bo 
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Qu 
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+ 
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1 ifs? n n v Vv = 
+ [lied + Te Gq —216°D,) ds — 


Bei Anwendung des Bettischen Theorems 


[Oads—[Peguds + [ Bre h0 as + [Qads— [ ee pads + [ Bue pads = 
1 2 3 at 2) 3 


= [a Was — te" 8, ds + foneGeds + iW ds—[tedads + [bn dads 
: 3 * 3 


a] 
o || 


c c ¢ c 


a! : 5 5 as - Vv 
auf (12) kénnen die n linearen Gleichungen fiir die Konstanten c; »=1, 2,..., ” in 


der Form 
[(e — 0) w ds =e — T) e% Ca ds + {(B-3) n Gads = 
us Es 
_ [lw —W) ds + [ie (i, ag ds — [ons pees) ds — (13) 
Es Z us | 
| 


— | {roa=o: pss lh, Bye cog 
angegeben werden. La8t man in (10) die Anzahl der zur Approximation verwendeten 
Funktionen iiber alle Grenzen wachsen, dann liefert (13) die Koeffizienten des gegen 


die exakte Loésung konvergierenden Ansatzes. 
(Bingegangen am 30. Juli 1959) 


264 Kleine Mitteilungen 


Kleine Mitteilungen 


Beitrag zur Gitterkorrektur nach Betz* 
Von H. Gallant, Graz 
Mit 9 Textabbildungen 


Im Mittelpunkt der Betrachtungen tiber axialdurchstrémte Maschinen steht das Schaufelgitter. 
Man gewinnt es durch Abwickeln von koaxialen Zylinderschnitten durch das Laut- bzw. Leitrad. 
Die Potentialstrémung in einem derartigen ebenen Gitter kann man mit Hilfe der konformen Ab- 
bildung oder nach dem Singularititenverfahren berechnen. Beide Wege erfordern jedoch einen 
gréBeren Rechenaufwand. Bei der konformen Abbildung wird die Strémung durch eine komplexe 
Abbildungsfunktion auf die bekannte Kreiszylinderstromung zuriickgefithrt. Mit weniger Aufwand 
verbunden ist das von H. Sehlichting und N. Scholz®-® entwickelte Singularitatenverfahren, 
das bei maBigen Schaufelwélbungen eine gute Niherung der tatsdchlichen Stromung ergibt. 

In vielen Fallen ist jedoch eine genaue Kenntnis des Strémungsbildes nicht erforderlich, es 
interessiert nur ein einfacher Weg zur Ermittlung eines Gitters, das eine bestimmte Ablenkung 
hervorruft. Hier ist die Anwendung des Verfahrens von Betz! zweckmafig. Dieses Verfahren 
gestattet es, einen Fliigel, dessen Kennwerte nach der Tragfliigeltheorie bei Parallelstromung 
bestimmt worden sind, in die durch den Einfluf8 der Nachbarfligel im Gitter abgelenkte Stromung, 
im folgenden Grundstr6mung genannt, konform einzufigen. 

Dem Verfahren liegt die Annahme zugrunde, da die Wirkung eines korrigierten Profils im 
Gitterverband dann gleich der Wirkung des nichtkorrigierten Profils in der Parallelstromung 
ist, wenn die Abweichungen des korrigierten Profils von der Grundstrémung gleich den Ab- 
weichungen des nichtkorrigierten Profils von der Parallelstr6mung sind. Diese Annahme 1labt 
freilich bei beschleunigter und verzégerter Stromung sowie bei Gittern mit starker Richtungs- 
ainderung Abweichungen erwarten?. 

Um die Grundstroémung zu finden, ersetzt man jeden Fligel entweder durch einen Einzel- 
wirbel von gleicher Zirkulation, der im Druckpunkt des Profils liegt (nach Betz 2. Naherung 
genannt), oder durch eine Zirkulationsverteilung entlang der Skelettlinie, die durch die Druck- 
verteilung gegeben ist (nach Betz 3. Naherung genannt). 

Ein Nachteil der Korrektur nach Betz liegt auch darin, da der Einflu8 der Dicke der Profile 
auf die Stromung im Gitter nicht bericksichtigt wird. 

Es ist nun der Zweck der folgenden Untersuchungen: 

1. den Einflu8 der Profildicke eines nach Betz korrigierten Profils auf die Gitterstromung 
festzustellen und 

2. die Ergebnisse der notwendigen Gitterkorrektur des 2. und 3. Naherungsverfahrens nach 
Betz zu vergleichen. 


I. Einflu8 der Profildicke 


Bekanntlich geben symmetrische Profile im Gitter bei stoSfreier Anstro6mung infolge ihrer 
Dicke immer eine Ablenkung zur Gitterfront. Dieser Dickeneinflu8, der von der Dicke der Profile, 
dem Staffelungswinkel 2 und dem Gitterverhiltnis ¢/l abhangt, ist unabhingig von der Wélbung 
der Profilmittellinie zu errechnen®. Diese Rechnung wurde von Eckert fir Joukowsky-Profile 
durchgefiihrt. Die Ergebnisse lassen sich jedoch, wie Messungen von Eckert (bei kleinen Mach- 
Zahlen) zeigen, auch auf Profile beliebiger Dickenverteilung tibertragen. 

Verwendet man diese Werte im Zusammenhang mit der Korrektur nach Betz, so kann man die 
Ablenkung infolge der Profildicke als Anderung der Null-Auftriebsrichtung eines Einzelprofils 
im Gitterverband auffassen. Der Anstellwinkel bei einem gewiinschten Auftriebsbeiwert wird also 
eine entsprechend geiinderte Groéfe haben. 

Um eine Kontrolle der nach diesem Verfahren — Gitterkorrektur nach Betz, jedoch mit 
Bertcksichtigung des Dickeneinflusses — errechneten Werte zu haben, wurden Beispiele durch- 
gerechnet und mit dem Singularitiitenverfahren nach H. Schlichting® verglichen. 

Die Vergleichsrechnung fiir ein diinnes Kreisbogenprofil mit //t= 0,05, t/l=1,0, ca = 0,72 
und 2= —30° ergab eine vollige Ubereinstimmung der Umlenkwinkel. Die Differenz betrug 20’, 
ein Wert, der unter der Rechengenauigkeit liegt. 

Bei dicken Profilen werden die Umlenkunterschiede gréBer; so ergeben sich bei einem Profil 
G6 389, ca = 0,73, t//= 1,0 und 2= —30° folgende Verhiltnisse, wobei das gewahlte Profil nach 
Betz unter Beriicksichtigung der Zirkulationsverteilung korrigiert und das so erhaltene Gitter 
nach Schlichting fir drei Aufpunkte nachgerechnet wurde. 

*) Aus dem Institut fiir Stromungsmaschinen II der Technischen Hochschule Graz, 
Vorstand Prof. Dr. E. Niedermayer. 
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’ Nach Betz Nach Schlichting 
ohne Dickenkorrektur 75° 30’ a 
ate Be 56° 20 ee 
mit Dickenkorrektur - £, 75° 30’ 75° 30’ 
54° 36’ 51° 48’ 


. . 2 
Die Differenz der nach den beiden Verfahren ermittelten Austrittswinkel betrigt also 2° 48’. 
Dieses Ergebnis rechtfertigt mit Riicksicht auf die Vereinfachung, die die Betz-Korrektur darstellt 
die Anwendung dieses Verfahrens, wenn der DickeneinfluB beriicksichtigt wird. 


> 


Il. Vergleich der 2. und 3. Naherung ohne Beriicksichtigung des Dickeneinflusses 


Es wurde ein Turbinengitter — Profil Gé 389 — mit und ohne Beriicksichtigung der Zirkulations- 
verteilung untersucht. 

Als Vergleichswert wurde der Winkel fx 
in den Punkten P, und P, ermittelt. P, 
bzw. P, sind Punkte der Skelettlinie des 
Profils, die auf der Ein- bzw. Austrittskante 
liegen. fx ist der Winkel zwischen der 
theoretisch der Profilberechnung zugrunde 
gelegten mittleren Durchstrémgeschwindig- 
keit wo. und der tatsachlich in den gewahl- 
ten Punkten infolge des Gittereinflusses 
herrschenden Richtung der Grundstrémung. 

Die GroBe der mittleren Durchstrém- 
geschwindigkeit wurde, um die Rechnung 
nicht untibersichtlich zu machen, im ganzen 
untersuchten Bereich konstant, und zwar 
Woo=10m/s, gehalten. Der EinfluB der 
Langenkorrektur der Fligeltiefe infolge un- 
gleicher Geschwindigkeiten langs des Profils 
wurde wegen seiner geringen Auswirkungen 
auf die Winkelkorrektur nicht beriick- 
sichtigt. 

Die Abb. 2 und 3 zeigen die Ergebnisse 
dieser Rechnung. Darin ist die Winkel- 
korrektur fx tiber dem Gitterverhaltnis ¢/1 
aufgetragen. Als Parameter wurde der Staffelungswinkel 7 gewahlt. Winkelkorrekturen, die zu 

starkeren Kriimmungen der Skelettlinie fihren, 


Abb. 1. Skizze der Gitteranordnung 


a wurden positiv, Korrekturen, die eine Streckung 
Yi der Skelettlinie hervorrufen, negativ aufgetragen. 
£ 
4 be 
Wed 
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Abb. 2. Korrekturwinkel fx fiir P; bei ca = Abb. 3. Korrekturwinkel Br fir P, bei Ca = 
—0,73. —— 2. Naherung; --- 3. Naherung = (0,73. —— 2. Niherung; ——- 3. Naherung 
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Bei groBen Staffelungswinkeln / und mittleren Werten von ¢/l werden die Winkelkorrekturen 
negativ. Dieses anfangs iiberraschende Ergebnis stimmt mit den Angaben in der Literatur*»? gut 


iiberein. Wie die Abbildungen zeigen, ist die 


Ubereinstimmung am Austritt fir ‘alle Gitter- 10 
verhaltnisse gut, am Eintritt aber ergeben sich 
im Bereich von ¢/l=0,6 bis t/I= 1,2 groBere e 


Abweichungen. 


“30 35 40 US 50 S5ALT 60 - "90 35 40 5 50 S5-ALT 60 
Abb. 4. Bx far P, bei t/l= 0,8 Abb. 5. Bx fir P, bei t/l=0,8 


In den Abb. 4 bis 9 sind die Winkelkorrekturen tiber dem Staffelungswinkel aufgetragen. 
AuSerdem wurden noch die Zusammenhange bei zwei weiteren Auftriebsbeiwerten ¢a= 1,1 und 
Ca = 0,33 eingetragen. Wie die Abbildungen zeigen, besteht am Schaufelaustritt fir alle Verhialt- 
nisse eine gute Ubereinstimmung zwischen den 
Ergebnissen der 2. und 3. Naherung. ‘Am Ein- 
tritt trifft dies nicht zu, auch ist die Gro8e und 
Richtung der Abweichungen sehr stark von den 
einzelnen Parametern (t/l, ca und A) abhiangig. 


0 3S 40 YS 50 55 -ALT 60 30 BF 40 4g £0 S5-A[7/ 60 
Abb. 6. Br fir P, bei t/l=1,0 Abb. 7. Bx fiir P, bei t/1=1,0 


Da jedoch die Eintrittskorrektur erfahrungsgema8 keine so groBe Rolle spielt wie die Korrektur 
der Austrittskante, sind diese Umstinde von geringerer Bedeutung. 
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Die Ergebnisse des untersuchten Turbinengitters gelten qualitativ auch fir Gitter, die aus 


anderen Profilen aufgebaut sind, da durch die Variation der Auftriebsbei i anzli 
‘ x ? aria £ t bs f 4, € v 
verschiedener Druck- und daher auch Wirbel- bidnschilp ntact aan aieicianie 


verteilungen berticksichtigt wurde. Die Werte be i a 
der drei Parameter A, t/l und Zirkulationsver- fy 
teilung (hier durch cq ausgedriickt), fiir die die i= +— 
2. Naherung nicht mehr zulassig erscheint, kén- N am: 
nenaus den Abbildungen entnommen werden. 
4 SJ Say ee 
5 IS 
Sal | ve 
SS N 
= BSS 


30 LTP 4S 50 55 -ALT 60 30 IS 4) 48 50 S5-ALY 60 
Abb. 8. 6x fiir P, bei t/l=1,2 Abb. 9. Bx far P, bei t/l=1,2 


In den Abbildungen 4 bis 9 bedeuten: diinne Linie: 3. Naherung; dicke Linie: 2. Niherung 


Ca = 1s 
FS Ca OSes 
i ny oy ie Ca = 0,33 


Wie man aus diesen Betrachtungen erkennen kann, ist die Gitterkorrektur nach Betz auch 
nach dem 2. Naherungsverfahren geeignet, mit geringem Aufwand bei vielfach zufriedenstellender 
Genauigkeit die aus den Messungen an Einzelfliigeln bekannten Profilwerte auf Gitter anzuwenden. 
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Finfiihrung in die Stré6mungsmaschinen. Turbinen, Kreiselpumpen und Verdichter. Von M. Adolph. 
Mit 206 Abbildungen und 32 Berechnungsbeispielen, VII, 259 S. Berlin-Gottingen-Heidelberg: 
Springer-Verlag. 1959. DM 21.—, Ganzleinen DM 24.—. 

Das Buch ist vor allem zur Beniitzung durch Hérer von Staatsgewerbeschulen gedacht. Mit. 
dieser Feststellung ist bereits implizit eine gewisse Charakterisierung verbunden. Das Buch bringt 
zwar auf einem verhaltnismakig knappen Raum eine beachtliche Fille von Stoff, der Menge des 
Gebotenen miissen aber natiirlich andere Dinge geopfert werden. So wird bei den meisten Her- 
leitungen auf wissenschaftliche Strenge verzichtet; vieles Begriindbare wird ohne Begriindung hin- 
geschrieben. Andererseits wird mit einer recht grofen Zahl von Beispielen offenbar Wert auf die 
Erlernung eines routinemaiBigen Berechnungsganges gelegt. — Uber den Wert oder Unwert solcher 
Darstellungen ist schon viel diskutiert worden. Mancher Praktiker wird sie befirworten ; der wissen- 
schaftlich geschulte Ingenieur wird sie ablehnen. 

IT. Bednarczyk, Wien 


Elementare Mechanik vom héheren Standpunkt. Von L. Féppl. Mit 71 Bildern und 8 Bildtafeln, 
174 S. Minchen: R. Oldenbourg. 1959. DM 20,—. 


Das vorliegende Buch ist in sehr anregender Form geschrieben und zeigt die reiche Lehrerfahrung 
seines Autors. Da gerade erst in den neueren Lehrbiichern der Mechanik etwas mehr auf die ge- 
schichtliche Entwicklung und den Kreis der groSen Forscher dieses Gebietes eingegangen wird, 
so dient dieses Buch dem mit den Anfangsgriinden der Mechanik vertrauten Leser als ausgezeichneter 
Weg zu Verstindnis und Wiirdigung der Ergebnisse dieses Wissensgebietes. 

Wenn im Rahmen dieser vorliegenden ,,Elementaren Mechanik‘‘ auch eine kurze Einfihrung 
in die Ideen und Grundziige der Relativitaitstheorie eingebaut wurde, so ist dies sehr zu begriiBen. 
Dann steht jedoch zur Diskussion bis zu welchem Ausmaf die Grundziige der Elasto- und der 
Hydromechanik in diesen Rahmen einzufiigen waren. Da sich der Autor jedoch auf die Mechanik 
starrer Gebilde beschrinkt und auf die beiden genannten Gebiete nicht eingeht, hatte dies in 
einem etwas enger gefaBten Titel des Buches zum Ausdruck gebracht werden konnen. 

Die Gliederung und Darstellung des behandelten Stoffes erweist den groBen EinfluB H. Machs 
auf den Verfasser. Dies zeigt besonders der als ,,Newtonsche Mechanik* betitelte II. Abschnitt, 
in dessen Rahmen iiber die von Newton begriindete Massenpunktsmechanik weit hinausgegangen 
wird. Durch neuere Quellenforschungen (siehe z. B. C. Truesdell, ZaMM 1958) erscheint eindeutig 
bewiesen, da den grundlegenden Uberlegungen Newtons durch Euler, Daniel und Johann Bernoulli 
und Cauchy nicht nur Verfeinerungen oder Verbesserungen, sondern eindeutig fundamentale 
Prinzipien — im besonderen zur Behandlung der Bewegung starrer Kérper — hinzugefiigt wurden. 
Diese Tatsache wird nicht klar ausgesprochen. 

Der Berichter hatte es als vorteiJhaft empfunden, wenn das vektorielle Produkt zweier Vektoren 
nicht durch eckige Klammern, sondern in mehr tiblicher Art durch ,, x ‘“‘ gekennzeichnet worden 
ware. Weiters erscheint die Bezeichnung des Dralls als ,,Impulsvektor“ nicht gliicklich gewahlt, 
wenn in den anderen Abschnitten des Buches die VektorgréBe ,,mv‘‘ als Impuls in tiblicher Form 
bezeichnet wird. 

Von den angegebenen, leicht zu behebenden Mangeln abgesehen, stellt das Buch zweifellos 
eine sehr anregende und wertvolle Erginzung des Schrifttums der Mechanik dar. Fir jeden ernsthaft 
an der Mechanik interessierten Leser wird es ohne Zweifel eine Bereicherung des Verstehens und 
einen Weg zur Wiirdigung der Leistungen der groRen Forscher bedeuten. 

A. Slibar, Stuttgart 


Kinfiihrung in die piezoelektrische MeStechnik. Von W. Gohlke. Zweite, bearbeitete und erweiterte 
Auflage. (Technisch-Physikalische Monographien: Band 8.) Mit 219 Abb., X, 278 S. 
Leipzig: Akademische Verlagsgesellschaft Geest & Portig K.-G. 1959. Geb. DM 26.—. 

Das Buch vermittelt dem Leser einen guten Einblick in die Physik der Piezoelektrizitat sowie 
in die MeBtechnik mit piezoelektrischen Kristallen. 

Der Verfasser beschaftigt sich zuniichst mit der Druckmessung, stellt mathematische Betrach- 
tungen uber die verschiedenen Mefisysteme und deren Schwingungseigenschaften an und geht 
nachher auf die Physik der Piezoelektrizitit ein. Aber nicht nur das Problem der Druckmessung 
wird erdrtert, sondern auch die Probleme der Kraft- und Beschleunigungsmessungen finden ihre 
Behandlung, so da8 alle Verwendungsméglichkeiten fiir die piezoelektrischen ,,MeSwertwandler“ 
in erschopfender Weise gebracht werden. Sehr vorteilhaft sind die zahlreichen Schnittzeichnungen 
uber den Aufbau der Mefiwertwandler, die dem Leser die zweckmaBigste Konstruktion vor Augen 
fihren. Der Verfasser geht ferner auch auf die Verstiirker zur Erhéhung der MeSempfindlichkeit 


Ma 
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pegs beschreibt auBberdem noch die Priifung sowie Eichung der gesamten MeSanordnung. Ein 
umfangreiches Schrifttum verzeichnis im Anhang zeigt, welche Miihe fir die Zusammenstellung 
des Buches aufgewendet worden ist. Die Lektiire dieses Buches kann bestens empfohlen werden. 


F. Brandsiaetter, Wien 


Design of cireular cylindrical shells. Von J. Heiand. Mit zahlreichen Textabb., 25 


. . “a . 3 
re Vitenskapsakademi Series 2, No. 3, 1957.) Oslo: Oslo University Press. 195 
Tr. oe : 


Dieses Buch befaf8t sich ausfiihrlich mit der mathematischen Elastizititstheorie kreiszylindri- 
scher Schalen. Das Vorwort bringt einen kurzen Uberblick tiber die Entwicklung der allgemeinen 
Schalentheorie seit Beginn dieses Jahrhunderts mit der Aufzihlung der wesentlichsten erschienenen 
Bucher und Aufsdtze. Der Verfasser weist darauf hin, da® die Erfahrung, gewonnen an Schalen- 
konstruktionen, die im Kinklang zu den elastizitatstheoretischen Ergebnissen entworfen wurden, 
gezeigt hat, da®B die wirklichen Spannungen und Dehnungen ziemlich gut mit den durch die 
Theorie vorausgesagten tibereinstimmen. Dies gab Anlaf, die Biegetheorie der Schalen ausfihrlich 
auszubauen. So ist auch der Inhalt dieses Buches der Verbesserung der mathematischen Methoden 
gewidmet, die von den drei simultanen partiellen Fliiggeschen Differentialgleichungen ftir die 
hacen ara u, v, w der Schalenmittelflache zu den tatsichlichen Kraften und Verformungen 

thren. 

Das Buch ist in sechs Abschnitte eingeteilt. Zuerst werden die drei Fliiggeschen Differential- 
gleichungen in tblicher Weise hergeleitet. Dann wird dargestellt, wie die Verschiebungen und 
hieraus alle statischen Groen durch eine Spannungsfunktion ausgedriickt werden kénnen. Die 
nachsten Abschnitte haben den Zweck, die Wirkung von Randbelastungen zu studieren. Die 
Spannungsfunktion wird in der gewohnten Weise durch passende Fourier-Reihen ausgedriickt. 
Die darauffolgenden mathematischen Entwicklungen werden weitgehender, ausfihrlicher und im 
groferen Umfang als bisher bekannt, bis zu den Ausdriicken der statischen GréBen durchgefihrt. 
Die gewonnene neue Naherungslosung wird mit den Ergebnissen anderer Autoren verglichen. 
Der Vergleich ergibt in den statischen GréBen eine unter 1% liegende Abweichung der neuen 
N4aherung des Verfassers gegentiber Pliigge, demgegentiber die Donnell-Niherung eine bei Extrem- 
fallen bis zu 8%ige Abweichung zeigt. 

An diese Behandlung des homogenen Teiles der Eingangsdifferentialgleichungen schlieBt sich 
der Abschnitt, in dem die verteilte Belastung in Betracht gezogen wird, und gezeigt wird, wie 
Partikularintegrale zu den Differentialgleichungen gefunden werden, wenn die Lasten in doppelte 
trigonometrische Reihen entwickelt werden. Der Fall der Belastung durch Eigengewicht wird im 
einzelnen behandelt und die Ergebnisse werden mit der Membrantheorie und mit den Partikular- 
integralen, gefunden durch Jenkins, verglichen. Kin weiterer Abschnitt befaBt sich mit dem 
Zusammenwirken von Schale und Randversteifungen an den geraden Rindern. Daran schlieBen 
sich dann im letzten Abschnitt des Buches Rechenbeispiele, die den Gebrauch der entwickelten 
Methoden in praktischen Fallen zeigen. E. Tschech, Graz 
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An Introduction to Plasticity. Von W. Prager. VIII, 148 8. Reading Massachusetts: Addison- 
Wesley Publishing Company, Inc., 1959. $ 6.50, 49s. 
Beim vorliegenden Buch handelt es sich um die englische Ausgabe des 1955 bei Birkhauser, 
Basel, erschienenen deutschen Originals. Es darf deshalb auf die in dieser Zeitschrift (Band 11, 
Seite 163, 1957) gegebene Besprechung verwiesen werden. EEDA TRATES 


Mechanische Schwingungslehre. Teil 1: Einfache Schwinger. Von Max Schuier. Zweite, bearbeitete 
Auflage. Mit 123 Abbildungen, VIII, 158 S. Leipzig: Akad. Verlagsges. Gees u. Portig. 1958. 
Geb. DM 15.—. 

Die zweite Auflage dieses bekannten Werkes stellt, bis auf einzelne Verbesserungen, im wesent- 
lichen einen Neudruck der ersten, schon seit einigen Jahren vergriffenen Auflage dar. Uber den 
Inhalt eines iiblichen Lehrbuches iaiber Schwingungslehre geht dieses Werk insofern hinaus, 
als viele Fragen iiber den Antrieb und die Ganggenauigkeit von Uhren sowie das Ingangsetzen 
eines gediimpften Schwingers erértert werden, Gebiete, in denen der Verfasser als Fachmann 
bekannt ist durch die von ihm erfundene photoelektrisch gesteuerte Uhr hoher Ganggenauigkeit 
und das nach ihm benannte Pendel. 

An mathematischen Kenntnissen wird nicht mehr vorausgesetzt, als jeder Student in der 
Mitte seines Studiums und jeder praktisch tatige Ingenieur haben soll, so da® dieses Buch sowohl 


als Lehrbuch als auch als Nachschlagewerk dienen kann. uaDeeover Bviea 
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Finfiihrung in die Technische Mechanik. Von I. Szabé. Vierte, verbesserte und erweiterte Auflage. Mit 
525 Abb., XII, 4348S. Berlin-Gottingen-Heidelberg: Springer-Verlag. 1959. Ganzleinen DM 22.50. 


An der vorliegenden vierten Auflage dieses so rasch beliebt gewordenen Lehrbuches wurden 
nur mehr kleinere Anderungen und Erginzungen vorgenommen. Die wichtigste besteht in der 
Aufnahme eines Anhanges, der eine Einfiihrung in die ,,Energieprinzipe’ der Mechanik enthalt. 
Fir Beweise wird auf die vom gleichen Verfasser stammende ,,Hohere Technische Mechanik‘‘ 
verwiesen. 


Das Buch kann dem Studierenden wirmstens empfohlen werden. 
H. Parkus, Wien 


Hiitte. Mathematische Formeln und Tafeln. Von J. Szabé. Herausgegeben vom Akademischen 
Verein Hiitte. Mit 142 Abb., 34 Tafeln, XV, 287 S. Berlin: Verlag Wilhelm Ernst & Sohn. 1959. 
Geb. DM 19.80. 

Die vorliegende, mit Tafeln versehene Formelsammlung soll dem in Forschung und Praxis 
stehenden Ingenieur, dessen Bedarf an mathematischen Kenntnissen immer noch anwachst, ein 
handliches Nachschlagewerk sein, das einerseits nicht zu umfangreich und unhandlich ist, anderer- 
seits aber doch ftir alle wesentlichen Bediirfnisse ausreicht. Formelsammlungen der hédheren 
Mathematik sind an sich schon diinn gesit. Die vorhandenen sind meist in erster Linie auf die 
Bediirfnisse des Studenten zugeschnitten, wodurch die Auswahl der Formeln und Tabellen eine 
ganz andere ist, als sie den Bediirfnissen des Ingenieurs entsprechen wiirde. Grundlage zu dem 
vorliegenden Band bildete der an sich ausgezeichnet zusammengestellte Abschnitt ,,Mathematik* 
in ,,Hiitte. Des Ingenieurs Taschenbuch‘, Bd. I, einem vielbewahrten Werk. Dieser Abschnitt 
wurde nunmehr zu einer selbstindigen Publikation verarbeitet und zu diesem Zweck in vielerlei 
Hinsicht erginzt und modernisiert. Insbesondere wurde der Abschnitt tiber Statistik neu ge- 
schrieben. 

Der Formelteil gliedert sich in folgende Abschnitte: Arithmetik, Kreis- und Hyperbelfunktionen, 
Differential- und Integralrechnung, Lineare Vektoralgebra (einschlieBlich Determinanten, Matrizen 
und Tensoren), Vektoranalysis, Analytische Geometrie (einschlieBlich Differentialgeometrie), 
Funktionen einer komplexen Verdanderlichen, Differentialgleichungen (einschlieBlich Randwert- 
problemen, Variationsrechnung, Integralgleichungen), Wahrscheinlichkeitsrechnung und mathe- 
‘matische Statistik, Praktische Mathematik, Inhalte von Flachen und Korpern. Von der ebenfalls 
sehr reichlich gehaltenen Sammlung von Tafeln, die natiirlich Potenzen, Wurzeln, Logarithmen 
und verschiedene wichtige Zahlwerte enthalt, seien insbesondere noch erwahnt: Elliptische Integrale, 
Loésungen wichtiger transzendenter Gleichungen, Nullstellen der Bessel-Funktionen, Bessel- 
Funktionen, Kugelfunktionen und Gammafunktion. Insgesamt enthalt der Tafelteil 26 Tafeln, 
deren jede von groBem Nutzen sein kann. 

Aus den kurzen Andeutungen wird man schon entnehmen konnen, da’ der Verfasser wirklich 
fir fast alle in der Praxis vorkommenden Falle das notwendige Riistzeug bereitgestellt hat. Eine 
besondere Leistung, wenn man bedenkt, da&® dies alles auf verhaltnismafig kleinen Raum kom- 
primiert werden mufte! Die Beniitzung der Formelsammlung wird wesentlich durch das Stich- 
wortverzeichnis erleichtert. Fiir Fille, in denen weiterfiithrende Literatur benétigt wird, ist auch 
noch ein verhaltnismaifig reichhaltiges Literaturverzeichnis beigegeben. Ein Hinweis auf alte 
geloste und ungeloéste mathematische Probleme kénnte, sofern er beachtet wird, so manchem, 
der sich immer noch mit Problemen wie der Quadratur des Kreises beschiaftigt, unniitze Arbeit 
ersparen ! 

Man kann wohl diesem Buch eine weite Verbreitung voraussagen. Es wird vielen Ingenieuren 
ein stets zuverlissiger und treuer Begleiter sein! 

E. Bukovies, Wien 


Allgemeine Schalentheorie und ihre Anwendung in der Technik. Von W. S. Wlassow. 661 S., 193 Bil- 
der. Berlin: Akademie-Verlay. 1958. Geb. DM 75.—. 


Das vorliegende Buch ist aus dem Russischen tibersetzt (wissenschaftliche Redaktion: Prof. 
Dr.-Ing. A. Kromm, Technische Hochschule Dresden/Graz) und bringt eine ausfihrliche Dar- 
stellung der Statik der Schalen. Die Stabilitatstheorie dimnwandiger Konstruktionen wird in diesem 
Buch nur kurz dargelegt. Den Inhalt des Buches bilden fiinf groBe Teilabschnitte. Sie sind das 
Ergebnis langjahriger Forschungstitigkeit des bekannten Verfassers. 

Der I. Teil ist der Membrantheorie der Schalen gewidmet. Er beginnt mit flichentheoretischen 
Grundlagen und der Aufstellung der auf die Hauptkriimmungslinien bezogenen Gleichgewichts- 
bedingungen der allgemeinen Membranschalen. Fiir die Rotationsschalen wird dieses Gleichungs- 
system im Falle des Fehlens der auferen Belastung auf eine Differentialgleichung zweiter Ordnung 
fir eine Spannungsfunktion zuriickgefiihrt. Daran schlieBen sich Betrachtungen, die den Anwen- 
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dungsbereich der Membrantheorie festlegen. Sodann wird gezeigt, da® die Grundgleichungen im 
Falle des homogenen Problems fiir die durch Flichen zweiter Ordnung mit positivem GauBschem 
Krimmungsmafi K gebildeten Schalen auf die Cauchy-Riemannschen Differentialgleichungen 
zuruckgefiihrt werden kénnen. Dies gestattet, die Schnittkrafte solecher Schalen mit Hilfe der 
funktionentheoretischen Methoden zu bestimmen. Fur eine Reihe neuer Aufgaben werden Lésungen 
in geschlossener analytischer Form angegeben, so fiir die elliptische Rotationsschale mit Einzel- 
kraften und -momenten, die in den Polen angreifen, mit Einzelkriften, die durch die Drehachse 
gehen und auf ihr senkrecht stehen, fiir die Torsion einer elliptischen Rotationsschale, fiir beliebig 
belastete elliptische und sphirische Kuppeln usw. 

Im II. Teil wird die allgemeine Theorie der biegesteifen Schalen behandelt. Ausgangspunkt 
bildet die Ableitung der Grundgleichungen der Elastizitaétstheorie in drei Dimensionen in orthogo- 
nalen krummlinigen Koordinaten. Diese werden zur Herleitung des Grunddifferentialgleichungs- 
systems einer beliebigen elastischen Schale herangezogen, wobei der Ubergang vom dreidimensio- 
nalen auf das zweidimensionale Problem unter Zugrundelegung der Kirchhoffschen Hypothese 
durchgefiihrt wird. Aus diesen allgemeinen Schalengleichungen heraus werden die Ausgangs- 
differentialgleichungen der Kreiszylinderschale und der Kugelschale angeschrieben und jedes dieser 
Systeme fiir den Fall der Belastung nur in Richtung der Schalennormalen auf eine Differential- 
gleichung achter Ordnung fir eine skalare Funktion zurickgefiihrt. Den Abschluf8 dieses Teiles 
bildet die Erstellung der allgemeinen Vertraglichkeitsbedingungen der Schalen. ; 

Im III. Teil entwickelt der Verfasser sodann die allgemeine technische Theorie derjenigen Scha- 
len, die entweder K = 0 haben (zylindrische und konische Schalen) oder fiir die X sehr klein ist, 
wie dies bei den schwachgeneigten Schalen der Fall ist. Als erste Anwendung wird die technische 
Biegetheorie der Kreiszylinderschalen in ausfithrlichster Weise behandelt. Unter verschiedenen 
Randbedingungen werden eine Reihe von Problemen der geschlossenen Kreiszylinderschale und 
soleche aus der Theorie der Tonnendicher behandelt. Daran schlieBt sich die allgemeine Theorie 
der schwachgeneigten Kugelschale mit Lésungen fiir verschiedene Sonderfille, wie: die Kugelschale 
unter einer im Scheitelpunkt angreifenden Einzellast, die Kugelschale unter gleichmabig verteilter 
Normallast, die schwachgeneigte, elastisch gebettete Kugelschale usw. AnschlieBend daran werden 
kompliziertere Aufgaben der Theorie schwachgeneigter Schalen behandelt. Die Losungen fithren 
zur praktischen Berechnung beliebiger schwachgeneigter, gelenkig gelagerter Schalen mit recht- 
eckigem Grundri®B. Die erzielten Ergebnisse dieser theoretischen Forschungen werden bei der 
Projektierung von diinnwandigen Deckenkonstruktionen angewandt. Der III. Teil schlieBt ab mit 
einigen Betrachtungen tiber die Stabilitatsbedingungen solcher Schalen. Sige 

Der IV. Teil enthalt die Theorie der orthotropen Kreiszylinderschalen, wobei die Biege- und Drill- 
momente in Lingsrichtung vernachlissigt werden, die Schale also in Lingsrichtung nur durch 
Membrankrafte beansprucht gedacht wird. Nach Herleitung der Grundgleichungen werden zuerst 
Probleme der Krafteinleitung in geschlossene unendlich lange Zylinderschalen behandelt. Dann 
folgt die Berechnung der geschlossenen Zylinderschale mittlerer Lange unter der Einwirkung 
radialer Belastung und bei verschiedenen Randbedingungen. Insbesondere werden die Spannungen 
i i i ine i 1 nkt angreie et 
in einer solchen Schale, die durch eine in einem beliebigen Punkt angreif nde radial gerichtete 
Einzelkraft belastet ist, ermittelt. Mit Hilfe der dargestellten Methode lassen sich die Spannungen 

i i li ) Angriffspunkt der Einzellast sehr genau 
und Deformationen in geniigend groBer Entfernung vom Angriffspu u ; 
bestimmen. In der Umgebung des Lastangriffspunktes mub man jedoch den Rano 
nach der genaueren Schalentheorie berechnen, die auch die Torsions- und die bee eee ene e 
beriicksichtigt. Das letzte Kapitel dieses Teiles beschaftigt sich mit der Darstellun g ae eae 
analytischen Methode zur Lésung von Randwertaufgaben fir offene Zylinderschalen. ie a = e 
ist sehr allgemein; sie gestattet, den Spannungszustand komplizierter sone paeey ac ef 
systeme zu untersuchen, die aus mehreren Zylinderschalen bestehen, fiir die oe en nae ae 
Langsrindern beliebig gegebene Randbedingungen gelten, die an den Mn Pan ee 
barten Schalen verstirkt sind und die sich unter der Emwirkung einer willkirlich gege 

ung befinden. ; 

ee v und letzte Teil des Buches enthalt eine Darstellung der allgemeinen Theorie op es 
In dieser Theorie werden mit Hilfe geeignet id fers a La seer Ie PU ye oS uae 
der zylindrischen und prismatischen orthotropen 5¢ halen au Ae ie fee et eee 
Differentialgleichungen zuriickgefihrt. Pees ee ae fe dee Randbedingungen an den 
prismatischer Schalen mit freien Liangsrandern und dann ne ae Paar cose ee 
Lingsrandern. In einem weiteren Kapitel wird pee en : Ero Ieee eee 
methode prismatischer Schalen dargestellt, die die Mog ichkeit g : eehiuinl nm Crt 
mathematischen Hilfsmitteln eine Reihe ziemlich komplizierter und pra 2 a eee toe 
aus der Theorie zylindrischer und Prem AE CIE: Pe eae ne ie ee 
aot olasniny Cane seas ear iete ae d indern der Schale beliebig vorge- 
ebenso kénnen die Randbedingungen an den Langs- und Querra g 2 
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geben werden. Die bisher angefiihrten Berechnungsmethoden prismatischer Schalen wurden bei 
Vernachlissigung der Schubverzerrungen durchgefiihrt. Im weiteren werden nun die Schubver- 
zerrungen beriicksichtigt. Die diinnwandige Konstruktion von der Form eines zweistéckigen Rah- 
mens, diinnwandige Kastentriger mit mehreren Zwischenwanden und der Kastentrager mit 
deformierbarer Querschnittsform erfahren eine ausfithrliche Behandlung. 

Das Buch steht auf hervorragendem wissenschaftlichen Niveau. Es vermittelt eine groBe Fille 
von wichtigen Erkenntnissen im Schalenbau. Durch die ausfihrliche und vorzigliche Darstellung 
der theoretischen Grundlagen und durch die groBe Anzahl von Losungen wird dem Schalenbauer 
die Méglichkeit gegeben. kompliziertere Fille selbst zu rechnen. Das Buch kann allen Fachkollegen 


warmstens empfohlen werden. 
E. Tschech, Graz 


Absolute Analysis. Von F’. und R. Nevanlinna. Grundlehren der mathematischen Wissenschaften : 
Bd. 102. Mit 4 Abb., VIII, 259. S. Berlin-Gottingen-Heidelberg: Springer-Verlag. 1959. 
DM 36.—, Ganzleinen DM 39.—. 

Es ist ein sehr dankenswertes Unternehmen der Verfasser, daf sie den von ihnen mehrfach vor- 
getragenen, systematisch und vielfach eigenstandig entwickelten ,,koordinatenlosen** Aufbau der 
Analysis, wie er sich in den letzten Jahrzehnten, etwa seit den Arbeiten von Hilbert, Banach, 
Riesz, Frechet und speziell J. Neumann angebahnt hat, nunmehr in einem ,,Lehrbuch* dargestellt 
haben. 

Die Darstellung beschrankt sich auf das Endlichdimensionale, wobei aber die Erweiterung auf 
einen Hilbert- oder Banach- Raum stets deutlich wird. Nach einer Einfithrung in die lineare Algebra 
wird zunachst die Differentialrechnung nach dem Frechetschen Ansatz entwickelt und deren 
Aufbau (Taylorsche Formel, partielle Ableitungen, Vertauschungssatz, implizite Funktionen) in 
diesem allgemeinen Rahmen gegeben. Die Integralrechnung schlieBt an das Integral eines alter- 
nierenden Differentials nach Cartan an, sodann werden die Integralsitze von Stokes und GauB 
mit einigen Anwendungen gebracht. 

Ein nachstes Kapitel gilt den Differentialgleichungen. Erst wird fiir gewohnliche Differential- 
gleichungen (Normalsysteme) der Eindeutigkeitssatz von Osgood formuliert und ein Existenzbeweis 
gegeben. Dann wird fur lineare Differentialgleichungen erster Ordnung (mit mehreren unabhangigen 
Variablen) ein Existenzsatz gezeigt, der — in Verwendung der Goursatschen Idee — nur die 
Differenzierbarkeit des linearen Operators verwendet. 

Im letzten Kapitel wird die Differentialgeometrie nach der allgemeinen Methode behandelt: 
zunachst die Kurven, dann die allgemeinen Flachen des Rn mit Verwendung der Tensorrechnung 
(Theorema egregium, Parallelverschiebung, Satz von Gauf-Bonnet). ; 

Die im Buch gegebenen Ausschnitte aus verschiedenen Gebieten geben ein eindrucksvolles 
Bild von der Kraft der koordinatenlosen oder besser absoluten Methode, mit der verschieden- 
artige Dinge einheitlich erfa8t, ein Formalismus gebildet und das Wesentliche in den Vordergrund 
gestellt wird. 

Moge das Buch dazu beitragen, diesen Methoden auch im Hochschulunterricht die ihnen zu- 
kommende Verbreitung und Anwendung moglichst bald zu sichern. 

H. Hornich, Wien 


Higentiimer, Herausgeber und Verleger: Springer-Verlag in Wien I, Mélkerbastei 5. — Fiir den Inhalt verant- 
wortlich: Prof. Dr. H. Parkus, Wien IV, Technische Hochschule, Karlsplatz 13. — Druck: Steyrermiihl, Wien VI, 
Gumpendorferstragwe 40—44, 
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Die in relativ kurzer Zeitspanne erforderlich gewordene sechste Auflage wird in einer umfas- 
_- genden Uberarbeitung vorgelegt ; richtungsweisend fir die getroffenen Ergadnzungen des in Statiker- 
__ kreisen des In- und Auslandes weit verbreiteten Standardwerkes waren die nachgelassenen Auf- 
__zeichnungen des Verfassers. ‘ 
___ Die im ersten Teil behandelten theoretischen Grundlagen wurden vor allem in bezug auf 
# _ Tragwerke mit gelenkigen Stabanschlissen ergainzt, deren Berechnung sich gegeniiber friheren 
: Auflagen erheblich vereinfacht. Zum besseren Verstandnis ist ein Kapitel iber das Wesen unver- 
packaeblioticn und verschieblicher Tragwerke mit zahlreichen Abbildungen eingefiigt; weiterhin 

- wurden fir Studium und Praxis tibersichtliche Tafeln mit einer Zusammenstellung der wichtigsten 

 Bedingungsgleichungen aufgenommen. is 

- _Der zweite Teil des Buches, der ausschlieBlich der praktischen Anwendung der im ersten Teil 
behandelten Rechnungsgrundlagen vorbehalten bleibt, ist hauptsdchlich durch solche Zahlen- 
beispiele erweitert worden, die bereits in dem Spezialwerk des Verfassers tiber die Cross-Methode 

- enthalten sind und somit einen Vergleich dieses Berechnungsverfahrens mit dem ,,Drehwinkel- 

_ verfahren‘‘ zulassen. 

. Die von allen Praktikern sehr geschatzte Zusammenstellung bequem benutzbarer Hilfstafeln 
im dritten Teil des Werkes, die sowohl fiir das Drehwinkelverfahren als auch fiir die Cross- 
Methode und zum Teil auch fiir die Methode der Festpunkte unmittelbar verwendet werden kann, 
wurde im wesentlichen durch Zahlen- und Kurventafeln sowie EinfluBlinientafeln fiir ,,Gelenk‘‘- 

__-baw. ,,Symmetriestibe“ bereichert. Damit wird eine rasche Ermittlung aller statischen Grundwerte 
= fir die Berechnung von Rahmentragwerken und Durehlauftragern mit und ohne Vouten ermoglicht. 


Hauptabschnitte 


_ Erster Teil: Rahmentragwerke ohne Vouten. — Rahmentragwerke mit beliebig verdnder- 
lichen Stabquerschnitten. — EinfluBlinien fir statisch unbestimmte Tragwerke. — Die Wirkung 
von Temperaturinderungen bei statisch unbestimmten Tragwerken. — Der Durchlauftrager mit 

- yerinderlichen Stabquerschnitten unter Beriicksichtigung aller Sonderfille. — Zweckmabige 
- Auflésungsverfahren fir lineare Gleichungssysteme. — Vereinfachte Berechnung hochgradig 
_gtatisch unbestimmter Tragwerke. — Verschiedene Methoden und Naherungsverfahren zur Be- 


7 - rechnung von Rahmentragwerken. 
F Zweiter Teil: Zahlenbeispiele. 


sé ritter Teil: Hilfstafeln. 

4 Zu beziehen durch Ihre Buchhandlung 
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Einleitung A 


Geschichtliches. Anwendungsgebiete des Dieselmotors. Kraftstoffe des Dieselmotors. Schmierstoffe des Dieselmotors. 
A 


bgrenzung des Bandes 12. Konstruktionsgrundsitze. Bezeichnung der Hinzelteile. © 
Allgemeine Richtlinien fiir den Aufbau von Dieselmotoren kleiner und mittlerer Leistuns 


Gesichtspunkte fiir die Wahl von Viertakt- oder Zweitaktverfahren. Grundlagen fiir den Bau von dtertaktanothrde e < 
ohne Kreuzkopf. (Wahl der Kennzahlen, Wahl der Zylinderzahl und Zylinderanordnung, konstruktiver Aufbau, 


normale Aufladung bei Viertakt-Tauchkolbenmotoren, Abgasturboaufladung, mechanischer Antrieb der Auflade- 


gruppen, fremdangetriebenes Aufladegeblise, Hochaufladung bei Viertaktmotoren.) Grundlagen fiir den Bau von 3 
Zweitaktmotoren ohne Kreuzkopf. (Wahl der Kennzahlen, Spiilung, Gleichstromspiilung, Umlenkspiilung, Spier, © 
Wahl der Zylinderzahl und Zylinderanordnung, konstruktiver Aufbau, Aufladung von Zweitakt- Tauchkolbenmotoren.) = 


Allgemeine Richtlinien fiir den Aufbau von Dieselmotoren groBer Leistung 


Gesichtspunkte fiir die Wahl von Viertakt-, einfachwirkenden Zweitakt- oder doppeltwirkenden Zweitaktmotoren. — 
Grundlagen fir den Bau von Viertaktmotoren in Kreuzkopfbauart. Grundlagen fiir den Bau von einfachwirkenden — 


Zweitaktmotoren mit Kreuzkopf. (Wahl der Kennzahlen, Sptilpumpen und Geblise, Aufladung beim einfachwirkenden 
Zweitaktmotor mit Kreuzkopf, konstruktiver Aufbau.) Grundlagen ftir den Bau von doppeltwirkenden Zweitakt- 
motoren. (Wahl der Kennzahlen, konstruktiver Aufbau.) Bezeichnung der Zylinder und Motoren. 

Einzelteile — 
Bauteile fir Tauchkolben-Vier- und Zweitaktmotoren. (Kurbelwelle, Kurbelwellenlager, Pleuelstange, Kolben mit 


- Kolbenringen und Bolzen, Grundplatte, Kurbelgehause, Zylinderblock, Zylinderbuchse, Zylinderdeckel, Zuganker, 


Steuerwellenantrieb, Einspritzpumpen, EHinspritzventil, Hinsaug- und Auspuffventile und ihr Antrieb, AnlaBein- 
richtung, Regler, Umsteuerung, Ol- und Wasserpumpen und ihr Antrieb, Geblase fiir Zweitaktmotoren, weitere 
Bauteile.) Bauteile fiir Motoren mit Kreuzkopf. (Kurbelwellen, Pleuelstange, Kreuzkopf und Gleitschuh, Kolben- 
stange, Kolbenstangenstopfbiichse, Kolben, Kolbenkiihlung, Grundplatte, Grundlager, Stander, Gleitbahn, "Zylinder- 
block mit Laufbuchse, Zylinderdeckel, Einspritzpumpe, Einspritzventil, Kolbenspiilpumpe, weitere Bauteile am 
cine Pyben ee: (AnlaBgefaBe, Olkiihler und SuSwasserriickkihler, polis tae Ome une Tilger, Drehvor- 
ric g. 
Beschreibung ausgefiihrter Motoren 


Kinfachwirkende Tauchkolben-Viertaktmotoren. (Motoren der Firmen: Maschinenfabrik Pi nach eee Nt inher AG., 
Werk Augsburg; Daimler-Benz AG., Untertiirkheim ; Kléckner-Humboldt-Deutz AG., K6In-Deutz; Maybach-Motoren- 
bau GmbH., Friedrichshafen am Bodensee; Motorenwerke Mannheim; Maschinenbau Kiel AG., Rurmeister & Wain, 
Kopenhagen; Gebr. Sulzer AG., Winterthir. ) Hinfachwirkende Tauchkolben -Zweitaktmotoren. (Motoren der Firmen: 
Maschinenfabrik Augsburg-N iirnberg AG., Werk Augsburg; Kléckner-Humboldt-Deutz AG., KéIn-Deutz; Maschinen- 
bau Kiel AG.; Motorenfabrik Darmstadt ‘AG: VEB Dieselmotorenwerk Rostock (DDR); J enbacher-Werke, Jenbach; 
Burmeister & Wain, Kopenhagen; Fiat Stabilimento Grandi Motori, Turin; Gebr. Sulzer AG., Winterthur.) Hoch: 
aufgeladener, einfachwirkender Viertaktmotor mit Kreuzkopf der MAN. "Binfachwirkende Zweitaktmotoren mit 
Kreuzkopf. (Motoren der Firmen: Maschinenfabrik Augsburg-Nirnberg AG., Werk Augsburg; Burmeister & Wain, 
Kopenhagen; William Doxford & Sons, Ltd. Sunderland; Fiat Stabilimento Grandi Motori, Turin; Gebr. Sulzer AG., 
Winterthur.) Doppeltwirkende Zweitaktmotoren. 


Sonderkonstruktionen 
Dieselnotstromanlagen. Dieselgasverfahren. — 


Schwerdlbetriech 


Viskositat. Ziindwilligkeit. Gehalt an Unverbrennlichem, an Schwefel und Asphalt. Verkokungsriickstand. Gehalt an 


organischen oder Mineralséuren. Beispiele von mit Schwer6l betriebenen poterchagecn. 


Beispiele von Dieselanlagen 


Ortsfeste Anlagen. Schiffsanlagen. (Direkter Antrieb des Propellers, Propellerantrieb liber Gettiere, Dieselelektrischer 
Schiffsantrieb.) 


Formelsammlung und Umrechnungstafeln 
Formelsammlung. (Leistungsformeln, Formeln fir Kolbengeschwindigkeiten.) Umrechnungstafeln. 
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